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Loo NEW 
Du Pont “SSS” Static- 


To meet the demand and need of 
seismic crews for an even better, safer 
electric blasting cap, Du Pont explo- 
sives technicians—through their pro- 
gram of continuing research—have 
now improved the popular, widely 
used Du Pont “SSS” Electric Blast- 
ing Cap. These caps have long been 
tops in the field . . . now they’re even 
better. Note these three new and im- 
portant advantages: 


































1. Much higher static-resistance. 


2. Dependable performance in shal- 
low-hole pattern shooting. 


3. No risk of premature firing 
through passage of high amperage 
stray current between cap shell and 
leg wires. 


Higher Static-Resistance. The key 
to high static-resistance in the ‘“‘SSS”’ 
cap is its conductive rubber plug. 
Means have been found to im- 
prove the efficiency of this plug 
and, as a result, the Du Pont 
“SSS” cap is now much more 
static-resistant than ever 
before. It cannot be exploded 
by voltages far higher than any likely 
to be encountered in the field. 
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ADVANTAGES in 
Resistant (SR) Electric Blasting Cans 


Shallow-hole Pattern Shooting. 
When this typeof work was carried out in 
highly conductive ground, current leak- 
age frequently caused cap misfires. To 
eliminate these, and at the same time 
maintain high static-resistance, was 
quite a problem. It has been solved, how- 
ever, and the new Du Pont ‘“‘SSS”’ static- 
resistant cap will give dependable per- 
formance in shallow-hole pattern shoot- 
ing. 

Premature Firing. There was a remote 
possibility of premature firing in the old 
“SSS” cap through the passage of a high- 
amperage current between the leg wires 
and the cap shell. This also has been 
eliminated in the new Du Pont ‘‘SSS” 
cap. 





Other Dependable 
Seismic Products by Du Pont 


“NITRAMON” S*—the safest blasting 
agent available. Cannot be detonated 
with ordinary blasting caps, open flame, 
or by sudden shock or friction. Readily 
detonated, however, with a “‘Nitramon’”’ 
S Primer . . . itself relatively insensitive. 
Packed in water-tight, threaded metal 


*Reg. trade-mark for Du Pont 
nitrocarbonitrate blasting ag 





cans that simplify making up strong, 
rigid charges. 


“‘HI-VELOCITY” Gelatin 60%—espe- 
cially recommended for deep-hole work 
and sleeper charges. Withstands high 
water pressure, prolonged immersion 
and adverse storage. It is available with 
Du Pont ‘‘Fast-Coupler.”’ 


SEISMOGEL—an economical gelatin for 
use under less severe water conditions. 
Also available with ‘‘Fast-Coupler.”’ 


SEISMOGRAPH BOOSTER— This is a new 
Du Pont product designed for use with 
gelatins and either ‘‘SSS” Electric Blast- 
ing Capsor Primacord. The Booster gives 
additional insurance that the charge will 
detonate under exceptionally adverse 
conditions of water-depth or exposure 
period. It will often give a completely 
successful shot, or at least fire enough to 
clean out the hole, where ordinary prim- 
ing will give nothing. 

Ask the Du Pont Explosives represen- 
tative serving your district for further 
information about these and any other 
Du Pont seismic products. E. I. du Pont 
de Nemours & Co. (Inc.), Explosives 
Dept., Wilmington 98, Delaware. 
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in a contest to obtain the design of an official crest symbolizing the science of EXPLORATION 
GEOPHYSICS. First prize will be a paid life membership, including subscription to Geophysics, 
in the society. 
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1. Entries, to qualify, must be received at the business office of the Society of Exploration Geo- 
physicists, Box 7248, Tulsa 18, Oklahoma, on or before March 1, 1952. Only active, associate or 
student members in good standing of the Society of Exploration Geophysicists shall be eligible 
to enter. Current officers and members of the Standing Committee on Honors and Awards are not 


eligible to receive prizes. 


2. Line drawings must be made with India ink on plain white paper or on tracing cloth. Lettering, 
if used, must be of sufficient size to be legible after reduction. Designs shall be judged on simplicity 
and clarity of meaning. They may be of any geometric form considered practical for use in the 
embellishment of award certificates, plaques, corporate seal or society letterhead. As a guide to 
the general nature of suitable designs, your attention is invited to the official crests of other societies, 
such as Geological Society of America, American Geophysical Union, American Institute of 
Metallurgical Engineers, Houston Geological Society etc., or the designs used by some of the 
geophysical companies whose advertisements appear in Geophysics. 
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ELASTIC WAVE VELOCITIES IN ROCKS AT HIGH PRESSURES 
AND TEMPERATURES* 


D. S. HUGHES} anp J. H. CROSST 


ABSTRACT 


The variation in dilatational and rotational velocities in rock samples with pressure and tempera- 
ture has been studied. The measurements were taken at pressures from 1 bar (atmospheric) to 
5,000 bars, and temperatures of -25°C. to 200°C. for all samples and up to 300°C. for two samples. 
We can thus reach a pressure equivalent to some 18 km or 60,000 feet of burial but our highest 
available temperature probably at most corresponds to no more than 8-9 km. 

The measured rotational velocities check very well with independent measurements. Laboratory 
measurements of the dilatational velocity are not available for comparison. From the measured 
dilatational and rotational velocities, the elastic moduli and Poisson’s ratio may be computed. 
Values of Poisson’s ratio are tabulated for all samples. In general, highly quartzitic rocks have low 
values 0.13-0.20 whereas the majority of rocks have values in the range 0.26-0.33. 

The effect of interstitial water has been investigated in one sandstone and one limestone. The 
sandstone shows an increase in velocity at low pressure and a decrease at high pressure whereas the 
limestone shows the opposite effect. 


INTRODUCTION 


A method of measuring both the dilatational and rotational wave velocities 
- that is applicable to small laboratory samples has been developed and has been 
in use for some time (Hughes, Pondrom, and Mims, 1948). It is particularly 
applicable to rocks and to measurements in pressure chambers (Hughes and 
Jones, 1950) since the velocities are obtained from measurements of the trans- 
mission time of pulses through the sample. It is thus possible to measure high 
loss materials to which resonance methods would either be inapplicable or 
inaccurate. With a few minor exceptions to be discussed later, the boundary 
conditions of the sample are unimportant. 


APPARATUS 


A block diagram of the pressure apparatus is presented in Figure 1. It consists 
of two systems, including a low pressure system consisting of a hand pump and 


* This work was supported by the Office of Naval Research under Contract N6onr-266, Task 
Order VIII, and by research grants from the Shell Oil Company and the Humble Oil and Refining 
Company. Manuscript received by the editor May 21, 1951. 

t Department of Physics, the University of Texas, Austin, Texas. 
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dead weight gauge that are used over the range 1.0-1,000 bars. At the latter 
pressure the high pressure valve is closed and an auxiliary pumping system is 
used to raise the piston in the pressure intensifier. The volume displacement of 
the intensifier is sufficient to raise the pressure in the sample and gauge cylinders 
to about 6,000 bars. The sample cylinder is 1.40 inches inside diameter and 7.0 
inchesfinside height. Two insulated electrical leads are brought through the 
chamber head. An electric furnace can be placed around the sample cylinder and 
the temperature raised to 300°C. 
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Fic. 1. Block Diagram Pressure Apparatus. 


A block diagram of the electronic equipment is presented in Figure 2. This 
has been modified somewhat from that described in the references. The essential 
method is the same but the number of components has been reduced by sub- 
stantially rebuilding the TS/100 AP range oscilloscope into a double oscilloscope 
capable of presenting A, J, and R displays. The oscilloscope delivers a trigger 
pulse at intervals of 2 to 3 milliseconds. These pulses are delayed about 5.0 uv 
sec by the delay line, amplified by the pulse amplifier, and impressed on the 
driving crystal. The signal from the receiving crystal is amplified and returned 
to the oscilloscope. The time between the trigger pulse and the received pulses 
can be read to about 0.02 yu sec on the oscilloscope. 
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The samples are cylindrical cores about 1.0 inch in diameter by 1.5-3.25 
inches long. The longer cores are preferred but good values of Vp can be obtained 
with short cores. These are cut from well cores or block samples with a diamond 
drill. The cores are usually jacketed by placing them in a close fitting copper 
tube, with wall o.o12 inch, and cutting the tube off even with the rock. The 
ends of the core are ground to provide a good surface and copper caps 0.010 
inch thickness placed over the ends. The caps are spun down on the tube and 
overlap about 0.3-0.4 inch. The seal between the caps and the tube has been 
made in several different ways. If the core is to be measured.to 300°C. a special 
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Fic. 2. Block Diagram Electronic Apparatus. 


lead-silver solder is used. This has the severe disadvantage that the core must be 
dry in order to make the seal. Dr. Birch has used rubber tubes for jackets and the 
authors have been very successful with neoprene tubing, but these cannot be 
heated much above 200°C. in the pressure chamber. Neoprene tubing has been 
used, however, for all cores with included water. Dr. D. H. Reno of the Stanolind 
Oil Company has suggested the use of Teflon tubing. This material will stand 
the temperature and pressure but is difficult to seal at low pressures. 


THEORY 


The appearance of the transmitted pulse is shown in Figure 3. This record 
was taken with a sample of -Solenhofen limestone in the pressure chamber. The 
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first arrival is the direct dilatational wave through the sample. The length of 
this sample is 8.252 cm. An arrival time of 13.56 u sec gives 6,086 meters/sec 
for Vp. At the higher pressures it will be noted that following the first arrival 
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Fic. 3. Pulse Arrivals Solenhofen Limestone. 


there are others at intervals of about 6.60 uw sec. These Jater pulses have been 
shown to arise from the incident dilatational pulse being transformed at the 
boundary and crossing the rod as a rotational wave and being transformed into 
a dilatational wave at the opposite boundary (Mason and McSkinin, 1947). 

In general we may write (Hughes, Pondrom, and Mims, 1948), 


tmn = (2n + 1)L/Vp + mD(V p? — Vr?)!!2/V Ve (4.4) 


as the times of arrival of the pulses from a single excitation of the driving crystal. 
Where L=length of sample 
D= diameter 
Vp=dilatational velocity 
aVr=rotational velocity 


m and n are integers subject to the condition 


(2n + 1)L = mD tan 0 (4.b) 
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where 
6 = sin“ Vr/V p. (4.c) 


There is no well defined limit to the values that m and m may take subject to 
the above condition. With metal rods of good acoustic properties, values of x 
as high as 8 or g and values of m of about 7 have been observed. With rocks, 
Figure 3 shows about the best results. The first arrival corresponds to n=o, 
m=o. The second arrival corresponds to n=0, m=1. The arrival with n=o, 
m= 2 is visible but not readable even at the highest pressures. At the very end 
of the traces, the arrival n=1, m=o is visible. This, of course, is the pulse that 
has travelled the length of the rod three times. 

Evidently if the arrival times of the 0,o and 0,1 pulses can be read, both Vp 
and Vr can be computed. Since these velocities are computed from pulse arrival 
time, they are actually the velocity corresponding to the band of frequencies 
necessary to form the break. In the present work this is approximately 1.0-5.0 
mc. 

RESULTS 


In spite of the high equivalent frequency, our velocities check very well with 
other measurements at lower frequencies, with one possible exception to be 
discussed later. Table I presents a comparison of our measurements of Vr with 
Birch (1942). 














TABLE I 
Vr (meters/sec) 
Pressure Solenhofen Limestone 
(bars) Birch Authors 
I 2,750 "2,854 
500 2,910 3,008 
4,000 3,080 3,066 








In addition to the Solenhofen limestone, samples of the Quincy granite, 
Twin Sisters dunite, Cheshire quartzite, and Danby marble have been measured 
by Dr. Birch and ourselves. Our value on the Quincy granite is about 10 percent 
low and apparently must be ascribed to a leak. Dr. Birch’s values for the other 
rocks have not as yet been published. Our value for the quartzite checks well 
with other quartzites studied and the one value for the dunite is close to that of 
Birch but this rock shows some anisotropy.* 

Dr. Birch’s measurements were made at frequencies of several kilocycles by 
a resonance method. Because of the damping under pressure by the transmitting 
fluid, Vp cannot be measured by this method. We have not been able to get a 
direct check in our measurements of Vp in rocks (Hughes and Jones 1950). 


* Personal communication. 
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The values measured so far have been very reasonable. A better check is probably 
the values of Poisson’s ratio computed from Vp and Ver. Our computed values 
are regarded as very reasonable and since a small change in the velocities will 
make a large change in this quantity, it constitutes an excellent check. Table XI 
shows computed values for Poisson’s ratio for several rocks under different 
conditions. It will be noted that quartzitic rocks tend to show a low value in 
general. 

Table II shows the Vp and Vez for dunite from Twin Sisters Mountain in 
Washington. This dunite, from thin section study, is go percent olivine and 
10 percent augite. The olivine crystals form an interlocking mosaic with a scat- 


TABLE IT 


Core #123, Sample 28 
Density 3.160 


Description: Dunite 
Origin: Twin Sisters Mtn., Washington 























T—C. 24.5 100 200 300 24 
Fig (bars) Vo Vr Vo Ve Vo Ve Vo Ver Vo Ver 
70 8,603 4,185 | 8,367 4,224 8,725 4,390 
138 8,694 4,244 | 8,434 4,289 | 8,155 4,133 | 7,701 4,016 | 8,788 4,425 
345 8,746 4,406 | 8,635 4,303 | 8,448 4,286 | 8,033 8,874 4,456 
690 8,820 4,435 | 8,748 4,434 | 8,691 4,442 | 8,404 4,320 | 8,907 4,467 
1,035 8,873 4,456 | 8,791 4,439 | 8,710 4,512 | 8,572 4,368 | 8,928 4,490 
1,725 8,905 4,515 | 8,821 4,470 | 8,783 4,467 | 8,652 4,444 | 8,973 4,510 
2,415 8,926 4,545 | 8,884 4,480 | 8,793 4,440 | 8,703 4,307 | 8,983 4,526 
3,105 8,936 4,553 | 8,862 4,502 | 8,834 4,472 | 8,743 4,422 | 9,017 4,530 
4,140 8,956 4,533 | 8,892 4,526 | 8,886 4,492 | 8,784 4,433 | 9,062 4,549 
5,170 8,977 4,529 | 8,087 4,559 | 8,895 4,514 | 8,792 4,459 | 9,006 4,539. 























tering of large augite crystals. It has apparently been metamorphosed. This 
tabulation shows several features that are characteristic of rocks in general. 
First, the increase in velocity with pressure is very rapid at the lower pressures 
and increases more slowly thereafter. Part of this is undoubtedly due to the 
closing up of pore space, thus forcing individual crystals into better contact. The 
effect of temperature is also much larger at the low pressures. This rock shows a 
fairly regular variation with temperature but most rocks, as will be shown, tend 
to change somewhat erratically with temperature. The room temperature run 
after the high-temperature run also shows a somewhat higher velocity. For this 
rock, the change is only about } of one percent but in many cases it is substantial. 
This effect is interpreted to indicate that subjection to pressure and temperature 
has compacted the rock. 

This interpretation cannot always be applied, however. The corresponding 
data for a sample of granite from Barrifield, Ontario, is shown in Table III. This 
rock from thin section study has a granitic texture and is composed mostly of 
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quartz and microcline with a rose color probably due to iron stain during altera- 
tion. Here a room temperature run after the 100°C, and 200°C, runs shows a 
slight “improvement” but after the 300°C. run it is evident that the elastic 
properties have been drastically altered. In addition to a substantially lower 
velocity, the initial steep rise has disappeared as though all pore space had been 
filled and the material altered to a low velocity form. No apparent change in the 
sample could be detected. 

This same effect has been observed in samples whose jackets failed and 
permitted the pressure fluid to penetrate the rock. In this sample, however, the 
gain in weight was only 42 mg. Thus the added fluid was at most only 4 parts in 
10,000. When the jacket was removed no trace of kerosene could be detected. 


. TABLE IIT 
Core #150, Sample 26 


Description: Granite 
Origin: Barriefield, Ontario 






































T—°C. 28 100 200 28 300 28 
P (bars) Vp Ve Vp Ver Vp Vr Vp Ve Ve Vr Vp Ve 
20 5,640 2,866 
5° 5,824 2,930 | 5,966 2,996 | 5,062 2,768 | 5,774 2,036 5.753 2,784 
100 5,875 2,961 | 6,o1r 3,045 | 5,100 2,865 | 5,850 2,077 5,749 2,797 
250 6,080 3,046 | 6,130 3,107 | 5,520 2,980 | 6,066 3,061 | 4,488 5,815 2,785 
500 6,215 3,109 | 6,243 3,159 | 5,061 3,083 | 6,235 3,120 | 4,838 5,831 2,803 
750 6,288 3,163 | 6,282 3,105 | 6,138 3,160 | 6,302 3,172 | 5,101 5,865 2,703 
1,000 6,336 3,155 | 6,325 3,107 | 6,237 3,166 | 6,347 3,170 | 5,223 5,898 2,813 
1,500 6,370 3,169 | 6,354 3,193 | 6,304 3,181 | 6,380 3,174 | 5,401 5,941 2,800 
2,000 6,384 3,209 | 6,372 3,202 | 6,317 3,103 | 6,415 3,105 | 5,526 5,907 2,814 
2,500 6,403 3,233 | 6,301 3,211 | 6,336 3,181 | 6,423 3,232 | 5,618 6,022 2,825 
3,000 6,412 3,215 | 6,405 3,226 | 6,354 3,186 | 6,438 3,244 | 5,676 6,053 2,850 
4,000 6,429 3,231 | 6,428 3,243 | 6,387 3,196 | 6,481 3,256 | 5,706 6,105 2,869 
5,000 5,447 3,228 | 6,446 3,244 | 6,420 3,224 | 6,505 3,210 | 5,874 6,191 2,910 





Table IV shows results in a sample of Solenhofen limestone. This sample 
again shows a slight “improvement” after the temperature runs. This sample 
did develop a very slight leak. After the runs the surface was slightly discolored 
and a section showed that a thin surface layer about 0.05 cm thick had been 
altered. This is also shown by the data. No o1 arrival could be detected at 300°C, 
and on the return to room temperature. We have observed that a layer of ma- 
terial around the cylindrical boundary of the core with a velocity somewhat less 
than that of the rock tends to destroy the transformation necessary to give or. 
A copper jacket o.o12 inch thick has no apparent effect. A neoprene jacket 
over the length of the core reduces the transformed wave below detectibility. 

Figure 4 shows to a still greater degree some of the effects discussed above. 
This is a sample of Pierre shale from the vicinity of Limon, Colorado. The 
original core was taken from a subsurface depth of 111 feet. The sample was 
jacketed and placed in the pressure chamber. Initially no energy was transmitted. 
This is a common occurrence and usually is due to the end caps of the jacket not 
being seated. The pressure was raised slowly and at 750 bars a fairly strong 
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arrival appeared and remained visible as the pressure was lowered to 50 bars. 
The velocity at this pressure was then 2,675 meters/sec as compared to 1,875 
meters/sec as determined by The Carter Oil Company from seismic field work. 
This difference is undoubtedly due partly to the fact that the sample was carried 
up to a substantial pressure before the velocity was measured; however it is 
highly probable that this is also partly due to dispersion. As will be shown, this 
sample compacted substantially under pressure and it is probable that in its 
original state frequencies in the megacycle band were very strongly absorbed. 

Starting from the 50 bar point, the pressure was increased step by step and 
the resulting velocities measured. At 1,000 bars, the velocity was 3,325 m/sec. 


TABLE IV 
Core #134, Sample 24 


Density: 2.656 
Description: Solenhofen Limestone. 
Origin: Bavaria, Germany. 









































T—°C 25 100 200 300 25 
P (bars) Vp Vr Vo Vr Vo Ver Vp Ve Vp Vr 
° 5,974 2,879 6,028 

100 6,000 2,952 | 5,892 2,049 | 5,756 2,890 | 5,547 6 ,046 
250 6,009 2,977 | 5,909 2,057 | 5,764 2,906 | 5,629 6,059 
500 6,025 2,991 | 5,925 2,962 | 5,784 2,928 | 5,691 6,080 
750° 6,038 2,909 | 5,932 2,972 | 5,799 2,930 | 5,721 6,088 
I ,000 6,051 3,006 | 5,049 2,070 | 5,814 2,037 | 5,744 6,105 
1,500 6,063 3,013 | 5,973 2,074 | 5,838 2,942 | 5,771 6,127 
2,000 6,079 3,018 | 5,989 2,986 | 5,853 2,944 | 5,802 6,143 
2,500 6,100 3,023 | 6,005 2,004 | 5,880 2,047 | 5,816 6,165 
3,000 6,108 3,026 | 6,017 2,998 | 5,892 2,048 | 5,844 . 6,177 
4,000 6,123 3,046 | 6,023 3,004 | 5,910 2,961 | 5,866 Q;, 212 
5 ,000 6,134 3,040 | 6,037 3,004 | 5,928 2,974 | 5,892 6,218 





On lowering the pressure to 750 bars, the velocity continued to increase reaching 
3,385 m/sec. This almost certainly indicates a time lag, i.e., the 1,000 bar pressure 
was not held long enough to establish equilibrium. The pressure was then re- 
turned to 50 bars and the velocity measured as 3,000 m/sec. Runs II and III 
were then made from 50 bars to 1,000 bars and return. On these runs the sample 
showed little hysteresis, i.e., the velocity at a given pressure was nearly the same 
going up or down. The velocities for these runs are very nearly the same. This 
would indicate that the rock was stabilized for this pressure range. Run IV was 
then made carrying the pressure up to 3,000 bars. The velocities duplicate runs 
II and III very accurately up to 1,000 bars. Beyond 1,000 bars the velocity 
increases rapidly to a value of 3,970 m/sec at 2,000 bars and then slowly to 
4,120 m/sec at 3,000 bars. On the run down the sample shows some hysteresis 

















ELASTIC WAVE VELOCITIES IN ROCKS 585 


down to 2,000 bars but beyond this point continues along almost the same 
straight line to a value of 3,915 m/sec at 750 bars. At this pressure one of the 
electrical connectors broke and the sample had to be removed. Later, run V 
was made. Here it became evident that a true first arrival was not readable below 
500 bars and hence data is plotted beginning at this point. This run gives every 
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Fic. 4. Vp vs. Pressure Pierre Shale. 


evidence of duplicating the original. At 3,000 bars, the connection was again 
broken and the sample had to be removed. It is evident, however that the sample 
has been fairly well stabilized for the range o—3,000 bars. The original length 
of the sample was 7.394 cm..After the first run to 3,000 bars it was 7.175 cm. 
and after the second run it was 7.158 cm. The diameter changed proportionally. 
The density increased from 2.20 to 2.42. 

Table V shows data taken with a sample of sandstone from the Caplen 
Dome, Galveston County, Texas, subsurface depth 6,500—7,000 feet. The porosity 
was measured as 5.1 percent. The core was dried and the standard runs taken. 
As would be expected in a material high in quartz, the temperature effect is 
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small. Vp curves rise rapidly at first and finally very slowly with pressure. Vr 
changes very slowly. Poisson’s ratio is low but increases substantially with pres- 


TABLE V 


Core #154, Sample 29 
Density 2.543° 
Porosity 5.1% 





Description: Sandstone. 
Origin: Caplen Dome, Galveston County, Sun Oil Company Well, Depth, 6,500-7,000 ft. 


Remarks: Dry. 



































T—C 27 100 200 30 
P (bars) Vo Ve Vo Vr Vp Ver Vo Vr 

° 3,672 2,963 

50 3,936 .2,483 3,601 3,977 3322 

100 4,042 2,514 3,819 31295 3,586 
250 4,367 2,504 | 4,089 2,472 3,801 4,030 2,583 
500 4,577 2,700 4,382 4,144 2,600 4,428 2,644 
75° 4,758 2,807 4,587 2,688 | 4,411 2,637 4,670 = 2,705 
1,000 4,870 2,847 4,737 2,751 45575 2,675 4,812 2,764 
2,500 5,020 2,884 4,929 2,833 4,852 2,797 5,047 2,868 
2,000 5,110 = 2,925 5,061 2,876 4,992 2,846 5,193 2,932 
2,500 5,167 2,044 | 5,151 2,922 5,102 2,879 5,278 2,970 
3,000 5,127 2,956 5,208 2,966 55157 2,922 5,341 2,992 
4,000 5,279 2,975 5,295 2,969 5,225 2,044 5,419 3,004 
5 ,000 5,363 2,971 51345 2,973 55312 2,975 5,461 = 3,025 

TABLE VI 


Core #155, Sample 29 
Density 2.606° 
Porosity 5.1% 


Description: Sandstone. 





Origin: Caplen Dome, Galveston County, Sun Oil Company Well, Depth, 6,500-7,000 ft. 
Remarks: Water saturated. 
































T—C. 27 100 200 28 
P (bars) Vo Ve Vo Ve Vo Vr Vo Vr 
° 4,453 2,319 4,487 
5° 4,548 2,432 4,213 4,507 
100 4,462 2,687 4,196 3,789 4,536 
250 4,469 2,761 4,217 3,802 4,577 2,404 
500 4,505 2,896 4,233 3,873 4,659 2,509 
75° 4,537 3,008 4,268 31932 4,739 2,603 
1,000 45553 3,104 4,307 3,960 4,788 2,688 
1,500 4,608 3,263 4,324 4,074 4,903 2,506 
2,000 4,651 35377 4,373 4,180 5,007 2,669 
2,500 4,689 3,468 4,411 4,295 5,070 2,730 
3,000 4,706 3,598 4,442 4,397 5,160 2,968 
4,000 4,782 3,651 4,535 4,600 5,252 3,105 
5,000 4,887 3,945 4,660 4,796 55323 3,176 
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sure. The return to room temperature shows a small deviation from the original. 

Table VI shows some data taken with a water saturated sample of the same 
material. The first room temperature run does not show the initial large rise in 
velocity. The reason for the 200°C. velocities rising above the 100°C. velocities 
is not evident. This together with the return to room temperature run would 
indicate that the sample had lost water. However, the sample lost no. weight 
during the run so this explanation is not possible. 
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Fic. 5. Comparison Vp vs. Pressure Caplen Dome Sandstone, Saturated and Dry. 


Figure 5 shows the initial room temperature runs for both the dry and 
saturated samples plotted together. In this material the water has increased 
substantially the velocity at low pressure and reduced largely the rate of increase 
with pressure. Qualitatively this checks former observations (Hughes and Jones, 
1950) on igneous rocks in which oil was permitted to penetrate the rock. The 
curves may be explained by assuming that in the dry core the pulse is travelling 
alternately through rock and air. With increase in pressure the air space is 
closed up and the velocity is essentially that through solid rock. In the saturated 
core, the pulse is travelling alternately through rock and water. Water has a 
much higher velocity than air hence the initial velocity is higher. An increase in 
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pressure now, however, cannot close up the pore space and hence the initial 
velocity increases slowly with pressure. 

The large change in room temperature velocity after the 200°C. run shown 
in Table VI can be ascribed to the water combining chemically with some material 
in the rock but it is almost impossible to adequately check this hypothesis. That 
the above explanation cannot always hold is shown by Figure 6. These data were 
taken on Solenhofen limestone. Core 133 was run saturated then dried and re-run. 
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Fic. 6. Comparison Vp vs. Pressure Solenhofen Limestone, Saturated and Dry. 


The re-run checks very closely with the data taken on Core 134, showing that 
the water resulted in no permanent change. Here the saturated velocity is smaller 
at low pressure but increases more rapidly with pressure. This is effectively the 
opposite of the behavior of the Caplen Dome sandstone. 

The porosity of the Solenhofen limestone is about 4 percent and thus is very 
close to that of the Caplen Dome sandstone which is 5 percent. There is a great 
difference in permeability however. In our method of saturating the rock, it is 
placed in a glass tube and evacuated until warming it to about 100°C. gives no 
rise in pressure above 0.002 mm of Hg. It is then sealed off and the tube immersed 
in distilled water that has been boiled and pumped to remove gas. The tip of the 
tube is broken permitting water to enter the tube and sample. The tube is then 
placed in the pressure cylinder and the pressure slowly raised to 5,000 bars and 
held for about one hour. The pressure is then slowly released. If any water is 
forced into the sample by pressure it will be revealed by an equal amount of 
kerosene in the tube after removal from the pressure cylinder. The Caplen Dome 
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TABLE VII 
Core #141, Sample 32 
Density 2.514 
Description: Travis Peak Sandstone. 
Origin: Lassater Field, Marion County, Texas, 7,571 ft. subsurface, Magnolia Oil Company Well. 





























T—C. 27 100 27 
P (bars) Vo Vr Vo Ve Vo Vr 
° 3,726 3,926 
50 3,951 4,222 4,099 
100 4,113 4,367 4,224 
250 4,581 4,767 4,724 
500 5,086 5,048 5,183 
75° 5,310 5,310 55342 
1,000 55373 55354 5,382 
1,500 5427 5,402 5,403 
2,000 55457 5433 5433 
2,500 5,483 5,408 55459 
3,000 5,500 5,484 55495 
4,000 51533 5,518 5,514 
5 000 55553 55557 5,552 





sample was effectively saturated by atmospheric pressure, that is, no kerosene 
entered the tube. With the Solenhofen limestone, however, almost all the water 
was placed in the sample by the pressure. 

It appears to the authors that two separate mechanisms operate when a 
fluid is placed in the rock. First, the replacement of air or gas in the pore space 
leads to behavior similar to that shown by Figure 5. This explanation is checked 
by our former work (Hughes and Jones, 1950). Here, however, the fluid was No. 


TABLE VIII 
Core #122, Sample 27 


Description: Danby Marble 
Origin: Rutland, Vermont. 
































T—°C. 26 100 200 

P (bars) Vo Ve Vo Ve Vo Ver 
7° 5,866 2,823 4,530 2,512 4,296 2,336 
140 6,060 3,071 4,944 2,673 4,653 2,550 
345 6,250 3,104 5,819 2,947 55326 2,750 
690 6,501 35147 6,290 3,061 5,938 2,965 
I ,035 6,551 3,170 6,409 3,104 6,235 3,055 
1,725 6,621 3,195 6,542 3,125 6,395 3,090 
2,415 6,651 3,118 6,576 3,138 6,436 3,101 
3,100 6,665 3,209 6,585 35145 6,460 3,106 
4,140 6,673 3,209 6,602 3,146 6,468 3,106 
5,170 6,660 3254 6,621 35133 6,474 3,096 



























































vez'e = =6g1‘9 £61‘€ zvo‘g were - 1282 °0 Sgz‘€ = zO1‘9g SZr°¢. . Lo0%9 SZ1‘f 6609 £Ee*e. SLr‘9 ooo S 
fgz‘E = Lgr‘g Ser‘t.. goo’o || gee‘. o60°9 |. a9e5f Zorto gl1°f = gg6‘S 1g1‘f zfo'g | eee‘t . SEx‘9g 000‘ ¥ 
Soc‘e . Sex‘g obi‘ 196‘S oz‘ LSo‘g Sgz‘€ Sz1‘9g rere .—yso's Sy1*€ Geotg | Ore‘sE gor‘g 000‘ 
Seo Sra". | Ser"  O266S 1rz‘€ = Elofg | ESzE zir‘g | grrse bz6‘S | gge‘E g66‘S | Erz‘E ggo‘g ooS ‘z 
oe ote'€ _ ror‘g | ovr‘ gggfS gg1‘f z10‘9 Lye*e  ecOofg | girst . veg‘S weese z6Z‘S | ggrfE 6So0'9g 000'z 
wa , 
) tice: -rlo'9 |. pert. . ghg’s Lyr*€ olg*s ziz‘E = rLofg | 660'F  gfg‘S g6r'e: . 186°S .|-.9$3°t - €So%q oo$ ‘1 
1S) ol1‘€ zH0‘9 SLo‘€ ogl*S ver'S: '-026*S o61‘f zfo‘g SSote. ..2£2°S o6g°S | vE1SS o96‘S 000‘ I 
as Gerce... Lente |S9e°s SSL*S | Err’e (S6e°S | Gers oat .| -Leort- fen °s Leg*sS | fer1°t  gt6*s oSl 
en ali‘t. -696°S gzo'S = fg9‘S Lgoff = gg ‘S SZrte- -99g6"s goo'e r129‘S ogl ‘Ss gz1‘f = g6g‘S oo$ 
s 
Ss oS1‘f 616‘S of6‘z SgS‘S gSofe adgl‘S Ber’ <Sobts gS6‘z 39 fgh‘S Gol‘S | tors | ¥19°S oSz 
= 1f1‘E =. 1gQ‘S g6giz = olb‘S tho‘f obl‘S 6goE = Leg‘S ob6‘z 6gb‘S Lyg‘S Igo‘ zbl‘s OoI 
sy bz1‘'f gSgis veo'S 9th *s SLof€ = 61g*S go6‘z = vEh‘S 1z9‘S | fvo'S o1L'S of 
Ma} Llo’e ~ r189°S ° 
= 
= 4A S4 4A 4 4A aA <4 cA aA TA 4A aA bad, | a4 (s1eq) J 
z Qz 00z ool Lz 00z OOoI lz "I.—L 
wn 
S 





"12M Auedurog [IO [24S “sexay, ‘AyuN0D uojydy ‘Burppaq 0} se[noIpuadiag *3} £06‘6-006‘6 *1-y suey Joureg :UIsZUC 
(ODFD %04) auojsaurry snosseyjisz1y : uondussag 





6€L°2 Ayisueg 
1S1# 9109 


XI TIVE, 


ob sna a NE Ee 








ELASTIC WAVE VELOCITIES IN ROCKS 591 


10 motor oil, the pressure was limited to 1,000 bars, and the temperatures were 
200°C. or less. In addition, however, there is certainly the possibility of chemical 
change in the material due to interaction of the fluid and rock. In almost every 
instance of a jacket leaking the leak has occurred on the 300°C. run at high 
pressure. In every instance this penetration of hot kerosene under high pressure 
has very substantially lowered the velocity. With a core of Danby marble the 
rock became very soft and could be crumbled easily. It would no longer transmit 
any energy in our range. We have also observed two very porous sandstones of 
g percent and 28 percent porosity that would transmit no signal at all when 


TABLE X 
Core #152 
Density 2.731 


Description: Argillaceous Limestone (70% CaCOs) Parallel to Bedding. 
Origin: Barnet Farms A-1, 9,900-9,903 feet, Upton County, Texas. Shell Oil Company Well. 



































T—°C. 26 100 200 30 
P (bars) Vo Ve Vp Vr Vo Ve Vo Vr 
° 6,107 3,191 
50 6,028 3,030 5,987 3,166 5,926 3,131 6,114 3,210 
100 6,046 3,134 6,020 3,177 5,931 3,145 6,119 = 3,214 
250 6,088 3,178 6,056. 3,189 5,954 3,162 6,156 3,226 
500 6,131 3, 205 6,104 3,209 5,995 3,168 6,187 3,248 
75° 6,155 3,226 6,134 3,220 6,036 3,201 6,224 3,254 
I ,000 6,173 3,239 6,164 3,225 6,047 3,217 6,236 3,262 
1,500 6,241 3,262 6,182 3,452 6,088 35233 6,248 3,282 
2,000 6,278 3,282 6,212 3,259 6,111 3,232 6,265 3,284 
2,500 6,296 3,201 6,223 3,256 6,140 3,229 6,290 = 3, 283 
3,000 6,321 3,290 6,234 3,265 6,151 3,233 6,328 3,205 
4,000 6,344 3,296 6,257 3,274 6,160 3,244 | 6,351 3,204 
5,000 6,368 3,291 6,293 3,282 6,207 3,253 6,368 3,205 





saturated. The 9 percent core gave good data when dried. The 28 percent core 
would still transmit no energy when dry. 

Table VII shows dilatational velocities for the Travis Peak sandstone. The 
or arrival could not be identified and hence Vz could not be computed. This 
sample showed a porosity of 9 percent. Saturated with water it would transmit 
no detectible energy. The data of Table VIII were taken after drying. The 
velocities are very close to those measured on the Caplen Dome sandstone despite 
the great difference in origin. 

Several runs have been made on the Danby marble most of which were 
invalidated by leaks. This material, which is apparently a very pure marble, 
shows great changes with only slight leaks. Table VIII shows the best data that 
have been obtained. 

Tables IX and X show data on two cores cut from an ictal limestone 
well core. These cores were cut one perpendicular and one parallel to the bedding 
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from very closely neighboring positions. The data shows a definite anisotropy, 
with the higher velocity parallel to the bedding. This checks qualitative earlier | 
observations, but this data is particularly good because of the close spacing of 
the cores. 

CONCLUSIONS 


A résumé of data on several rocks is presented in Table XI. The Vp, Vr, and 
Poisson’s ratio o have been computed for various depths using a temperature 
increase with depth of 30°C./km. This gradient was chosen as an approximate 
average. The increase of Vp with depth down to 9,000 feet is very much less 
than is usually observed in the field. Below 9,000 feet the velocities increase 
very slowly or even decrease slightly. 

- A reasonable explanation would appear to be that in our measurements 
insufficient time has elapsed to attain equilibrium. This time could, of course, 
be extremely long. Observation thus far shows no detectable change in a few 
hours with the exception of the one point of the Pierre shale (Figure 4). 

The values of Poisson’s ratio are very reasonable and in general change but 
little with depth. The ratio for the Cheshire quartzite (147) is slightly irregular. 
The values for the quartzite, sandstone (154), and shale (116) are substantially 
lower than the other rocks. 

Samples 151 and 152 were cut from a single well core (Shell Oil Company, 
Barnet Farm A-1, 9,900 feet) perpendicular and parallel to the bedding. There is 
a substantial difference in the velocities. Our methods of computing Vr and o 
are not valid for anisotropic materials and thus the values entered for these 
cores are for comparison only. 
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A SEISMIC WAVE GUIDE PHENOMENON* 


K. E. BURG,t MAURICE EWING,{ FRANK PRESS,} E. J. STULKEN{ 


ABSTRACT 


On one particular prospect in shallow water repetitive patterns appeared on short spread seismo- 
grams in such prevalence as to jeopardize identification of desired reflections. It is demonstrated that 
under favorable conditions, less restrictive than thought necessary heretofore, a layer of water com- 
prises an effective wave guide for seismic energy propagation. Reinforcement fronts formed by multi- 
ple reflection of sound in water can develop into a set of waves completely overshadowing other seis- 
mic arrivals. With but minor modifications conventional wave guide theory applies. 

Examples from the prospect are presented to illustrate various reinforcement patterns. Observed 
frequency characteristics, group velocity, and phase velocity magnitudes are investigated for normal 
modes of propagation. 


INTRODUCTION 


Oil prospecting over water has revealed sound transmission characteristics 
either absent or less obvious on land. Among these can be mentioned repeated 
bursts of energy traceable to oscillations of the explosion bubble, simple multiple 
reflections between bottom and surface of the water, offside reflections from 
objects near the surface, and wave guide propagation which usually involves 
dispersion. 

This last characteristic, wave guide propagation, stems from reflections be- 
tween top and bottom of the water layer. The high frequency of first energy 
arrivals associated with long range direct waves in shallow water comprises one 
manifestation. Seismic observations reveal that this initial high frequency wave 
marks the beginning of a band of waves whose frequency rapidly decreases. 
These waves of continuously lower frequency persist until terminated by the 
Airy phase. The theory for the case of a liquid bottom has been given by Pekeris 
and for a solid bottom by Press and Ewing. In both these cases the propagation 
is due to total reflection within the liquid layer involving rays whose angle of 
incidence at the bottom lies between the critical angle and the grazing angle. 

In this paper we deal with a closely related phenomenon involving more nearly 
vertical incidence and very gradual dispersion. No doubt in most areas low 
reflection coefficients in effect at the bottom of water layers prevent the develop- 
ment of this wave guide action. The example here described came to light 
through two special circumstances which favor real or apparent sustained rein- 
forcement at practically constant amplitude level. First, over the seismic prospect 
which yielded the observations, the bottom of the water layer is smooth and 


* Presented at the St. Louis meeting of the Society April 25, 1951. Manuscript received by the 
editor July 16, 1951. ; 
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composed of rock, producing a high reflection coefficient even for normal inci- 
dence; second, automatic volume control on the amplifiers increases the gain for 
waves whose energy decreases due to leakage during a large number of reflections. 
This gives the impression of a long train of waves of almost constant amplitude. 

Previous theoretical work can readily be extended to include the case of 
only partial reflection at the bottom. It will be shown that in this case there is a 
limiting frequency for which the group velocity approaches zero and the phase 
velocity approaches infinity. Normal mode propagation in simply the water 
layer produces seismograms showing waves which could easily be mistaken for a 
long series of excellent reflections from deep horizons. 


OBSERVATIONS 


Characteristic Seismograms 


Since our introduction to the problem came through observation of seismo- 
grams on which regular reflection arrivals were obscured by water-transmitted 
energy, perhaps this discussion can be best introduced in a similar manner. 
Consider first then the seismogram illustrated in Figure 1, which may at first 
give the erroneous impression that it is a very good regular reflection record. 
The seismogram illustrates, among other things, the unusually long time through- 
out which wave amplitudes persist at levels that, because of automatic volume 
expansion, appear to be practically constant. Certain characteristics of the 
seismogram appear to repeat. A check on time intervals between them establishes 
periodicity. Further examination shows that almost every detail of trace move- 
ment is repeated over period after period. Essentially there is very little on the 
record beyond the features contained in a single interval such as the one indicated 
by marks. 

Figure 2 shows a seismogram which, from the standpoint of seismic pros- 
pecting, is more obviously hopeless than the first and therefore less treacherous. 
This second seismogram is simpler, consisting mostly of a long train of nearly 
sinusoidal waves. It develops that Figures 1 and 2 involve the same pattern of 
propagation; they differ only in the particular normal modes present. 


Empirical Formulas 


Examination of a number of seismograms over a particular offshore prospect 
indicates frequency of this wave motion is closely related to depth of water. Trial 
shows the frequency formula 


f = 3V/4H (1) 
describes the oscillations in Figure 2 well within the limit of accuracy with which 


water depth H is determined by sounding line or fathograms. Here V is velocity 
of sound transmission in water. In general, on many other records, 


f = (2n — 1)V/4H ; (2) 
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where n=1, 2, 3..., fits the observations well. These are essentially the fre- 
quencies of resonance in vibrating bars clamped at one end or in organ pipes 
closed at one end. 


Examples 


Figure 3 shows a resonant record of relatively low frequency—about 22} 
cycles per second. Applying the formula for the second normal mode, n=2 in 





Fic. 3. 16 to 40 cycle per sec seismogram showing second odd harmonic or second normal mode 
of vibration. Depths: Water, 165 ft; seismometers, 8 ft; dynamite, 5 ft. Distances from shot: Trace 1, 
2,555 {t; Trace 12, 2,275 ft; Trace 24, 2,480 ft. Recording spread: 2,000 ft. 


Equation (2), and using 5,000 feet per second for propagation velocity gives 
H=166 feet. The fathometer reading was 165 to 167 feet. Figure 4 shows reso- 
nance in about 50 feet of water. Frequency here is 26.8 cycles per second. Com- 
puting depth, assuming the fundamental mode of vibration, gives 47 feet. Evi- 
dently the fourth normal mode appears in Figure 5; the 51 cycle frequency 
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Fic. 4. 16 to 60 cycle per sec seismogram showing fundamental or first normal mode of vibration. 
Depths: Water, 53 ft; seismometers, 8 ft; dynamite, 5 ft. Distances from shot: Trace 1, 1,015 ft; Trace 
12, 420 ft; Trace 24, 1,165 ft.-Recording spread: 2,000 ft. 




















K. E. BURG, MAURICE EWING, FRANK PRESS, E. J. STULKEN 


Ay thy valey 
iy 
Hl Wit ait qheaVaghge 
ti) nt hy Al} iia 
} i li = ; f iM) de pe 
Vl mn NNN 
: | 0 


Fic. 5. 20 to 60 cycle per sec seismogram showing the fourth normal mode of vibration. 
Depths: Water, 164 ft; seismometers, 10 ft; dynamite, 5 ft. 





applied in the appropriate formula gives 172 feet for H. The sounding records 
‘ indicate a water depth of about 165 feet. In Figure 6 one finds the fifth normal 
j mode. In deep water even extremely high frequencies recorded with wide pass 
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Fic. 6.16 to 60 cycle per sec seismogram showing the fifth normal mode of vibration. Depths: 
Water, 214 ft; seismometers, 30 ft; dynamite, 214 ft. Distances from shot: Trace 1, 1,180 ft; Trace 12, 
500 ft; Trace 24, 1,070 ft. Recording spread: 2,000 ft. 








bands fit the resonance formula surprisingly well. In Figure 7 are shown sustained 
high frequency vibrations of many modes or harmonics. Recorded depth was 
106 feet. Tentatively identifying an observed periodicity of 4.3 milliseconds 
with the eleventh normal mode and computing depth yields 113 feet. 

As would be expected, the frequency response curves of our instruments 
determine roughly which modes appear. At least the range of frequencies that 
may be expected on a given seismogram narrows down to the limits of the pass 
band used. Although generally harmonics lying within the pass band are repre- 
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sented on a seismogram, it is sometimes difficult to account for the predominance 
of some and the absence of others. A typical frequency response curve appears 
in Figure 8. The arrow tips under the curve indicate a sequence of modes cor- 
responding to the water depth above which the seismogram shown at the right 
was shot. On this record one can identify the third and fourth normal modes, 
which lie in the frequency range receiving optimum amplification. 





Fic. 7. Unfiltered seismogram showing various higher harmonics. Depths: Water, 106 ft; seis- 
mometers, 8 ft; dynamite, 5 ft. Distances from shot: Trace 1,2530 ft; Trace 12, 2,160 ft; Trace 24, 
2,065 ft. Recording spread: 2,000 ft. 


Figure 9 shows a portion of a seismogram obtained in 248 feet of water with 
a very flat uniform bottom. The cutout fitted to the top traces is a mathe- 
matically computed curve in which third, fourth and fifth normal modes were 
computed from the frequency formula and simply added. One can easily see that 
the curve applies almost equally well to any appropriate part of the seismogram 
segment shown. 

THEORY 

General Remarks 


The unusually high multiples of reflections within the water layer offer evi- 
dence that reflection coefficients are close enough to unity so that even after 
many reflections their attenuating effect is substantially offset by an automatic 
increase in amplifier gain. Within the range of observations we have no specific 
data on amplitude level—only the illusion that oscillations continue practically 
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Fic. 9. The cutout at the top represents a mathematically computed curve of a combination of © 
third, fourth, and fifth normal modes of vibration in water 248 ft deep. Note correlation with the 
actual seismogram below it. 


undiminished. In first consideration of theory on an elementary level, then, the 
situation may as well be considered without reference to the effect of reflection 
coefficients, wave divergence, diffraction and other factors affecting attenuation. 


Repetitive Arrivals of a Pulse 


In this discussion we shall consider the shot a point source at (0, d) in an 
appropriately placed rectilinear coordinate system (x, z) with z=o repre- 
senting the surface of the water and z= ZH representing a level bottom. 


0,0 x 





sf 





(o,d) 


(0,2) + (x, 2) 








Y 





Fic. 10. Rays from (a, d) to (0, Z) multiple reflected at top and bottom of water layer and a particular 
ray path from (0, d) to (x, z) involving two reflections off the free surface. 
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Along «=o, a vertical line through the shot point, normal incidence applies 
upon reflection of the pulse at the boundaries of the water layer. It follows from 
the law of reflections that through successive reverberations the ray path of a. 
pulse is here retraced repeatedly. A particular point (0, Z) in water H units deep 
therefore receives the pulse alternately from above and below until the energy 
dies down. See Figure to. 

Other points, (x, z) in general, are also bombarded by a succession of pulse 
arrivals from a single explosion, but the mechanics are somewhat different. A 
particular ray path applies at most only once, at a single instant of time. To 
present the matter in more detail, four types of propagation from (0, d) to (x, z) 
are shown in Figure 11. The classification is based upon the direction of the 
vertical component of wave front propagation initially at (0, d@) and upon arriving 
at (x, z). Thus, Pa designates a ray which starts upward and passes (x, 2) 











Fic. 11. Four types of ray paths from (0, d) to (x, s): Pa, initially upward, terminally downward 
propagation; P,, initially upward, terminally upward; Qa, initially downward, terminally downward; 
Qu, initially downward, terminally upward propagation. 


downward; P,,, a ray which starts upward and passes (x, z) upward; Qu and Q, 
represent rays which start downward and pass (x, z) downward and upward 
respectively. Thus, for a single pulse at (0, d) four not necessarily separate energy 
arrivals are indicated at (x, z) for any given number of reflections off the top. 
Furthermore, each of the four types of paths provides multiple arrivals in direct 
relation to the different number of reflections off the free surface. 

Successive ray paths to (x, 2) are indicated in Figure 12 for Qa paths. If m 
denotes the number of reflections a given ray experiences off the top surface, it is 
apparent that the angle the ray makes with the vertical can be associated with 
the discrete number m and designated as 09,,,. Assuming a constant velocity of 
propagation V, path length is V¢ where / is travel time, so this incident angle 
can be expressed in the two following ways: 

x x 


= sin“! —. 
amH —d+z Vi (3) 
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(x,2) 








Fic. 12. The Qg wave is schematically illustrated by ray diagrams for one through five reflections off 
the free surface between (0, d) and (x, z)—that is for m=1, 2, 3, 4, and 5 respectively. 


Figure 12, as well as these expressions, shows convergence of successive rays 
through (x, z) toward the vertical as a limiting direction. 0n4:—0, thus becomes 
increasingly small. Ray length corresponding to 69,,, is also (2mH—d+z) 


COS DQ am +% SIN BQ: 
Since similar statements apply for the Qum, Pam, and Pum rays the situation 


may be summarized as follows: 


Angle of Incidence Travel Time Path Length 


tp, =*/V sin dp, Vir,,= [2(m— 1)H+d+z] cos Op, +x sin 6p, 





x 
Op, =tan 


—] _ 
& 2(m—1)H+d+z 











x ‘ “ 
6p, ,, = tan 1 ee Ee tp, =x/V sin Op, Virin= [2mH+d—s] cos Op, +x sin OP. 
x . , . 
994m = tan aie dae tom =*/V sin 0Q,,, Vloam= [2mH —d+s] cos 901m + Sin 89,,, 
x 
09, =tan™ tg =x/V sin@ Vig. =|2(m+1)H—d—z| cos 0 x sin 0 
= mba) Hadas on =2/V sinde,,, Vbo,.,™ l2(m-+1) ] COS 6 0p: +2 S12 60 


For any one of these four types of ray paths some convenient pulse charac- 
teristic, such as the extremes of pressure, plotted for successive arrivals at 
(x, z) with respect to time would appear somewhat as shown in Figure 13. A 
180° phase shift with each reflection off the free surface requires that alternate 
arrivals be plotted below the axis. 


Pulse Characteristics 

In the case of the artificial seisms which produced the records previously 
shown, the origin of each disturbance was a dynamite explosion. The work of 
Cole and others indicates an underwater explosion of this type is propagated 
outward in the form of a pulse with roughly the shape illustrated in Figure 14, 
where pressure is indicated qualitatively with respect to time. 


(4) 

















OBA PRI RSA EAD MRE TR 








K. E. BURG, MAURICE EWING, FRANK PRESS, E. J. STULKEN 








Fic. 13. Energy arrivals at (x, z) for a particular type ray for successively higher orders of 
multiple reflections. Pressure is plotted schematically with respect to travel time. 


Experience and theoretical work by other writers indicates that this pulse 
is essentially a complex wave form representable by a continuous frequency 
spectrum of sine waves. In this discussion it will simply be assumed that the 
arrivals at any point (x, z) have the same general shape as the original pulse 


r 














Fic. 14. Generalized pressure-time curve of an explosion. 


for a few reflections—for small values of m. As m increases these arrivals are 
manifested by whatever frequency components are favored by interference 
effects between the multiply reflected waves and the frequency response of the 
recording system. 


Wave Development 


In line with this view on pulse characteristics one would then expect the 
multiple arrivals portrayed in Figure 13 to appear in the form of those wave 
components which fit the pulse arrival pattern. The lowest frequency range 
expected is shown fitted to these arrivals in Figure 15. Development of even 
harmonics of this wave will evidently be prevented by the alternations in phase. 
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Fic. 15. Energy arrivals of Figure 13 fitted with an appropriate wave form. 


—_ 








Odd harmonics may, however, be fitted as well as the fundamental. The com- 
bination of eligible wave components may therefore be a quite complex wave 
form. Thus, Figure 9 shows the form developing out of the combination of three 
particular odd harmonics or normal modes. To investigate this wave development 
further let us express the four types of travel in the form of velocity potentials 
of the expected waves at any point (x, z) at a time ¢ after m reflections from the 


free surface.! 

Pam = A(—1)™e""! exp iar{ x sin @ + [2(m — 1)H + 2+ d] cos@ — Vi}/L, 
Pum = A(—1)™e" exp i2r{« sin 6+ [2mH — z+ d] cos 6 — Vt}/L, 

Qam = A(—1)™e™ exp i2r{ x sin 0+ [2mH + z — d] cos0 — Vt} /L, (5) 
Qum = A(—1)™e"*! exp i2r{ x sin + [2(m+ 1)H — z — d| cosé@ — Vit}/L, 
where 06=sin“ «/Vt, L is the wave length measured along the oblique ray path, 
A is a variable amplitude factor which is left unevaluated, and ¢ is the bottom 


reflection coefficient depending on 6 and the elastic properties of the two layers 
involved. Phase shifts during m reflections off the free surface are accounted for 


by the factor (—1)”. 
Let 


6 = exp i2m(2H cos 6)/L = exp ir(4H cos @)/L. (6) 


Then 
Pam = A(—1)™e™—15™— exp i2m[x sin 0+ (z + d) cos 0 — Vt|/L, 
Pum = A(—1)™e"5™ exp izr[x sin 6 — (z — d) cos @ — Vt]/L, 
Quam = A(—1)™e™5™ exp i2r|x sin 6+ (z — d) cos 6 — Vt]/L, 
Qum = A(—1) ™e™*1§™+1 exp iar|x sin 6 — (z+ d) cos @ — Vt]/L. 


(7) 


1 These wave expressions could, with minor modifications, be applied to direct analysis of other 
characteristics besides velocity potential, such as pressure, displacement, or vertical component of 


particle velocity. 
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The disturbance at any point (x, z) is due to the superposition of the various 
P and Q rays. Let us combine these as follows: 


Qa = Pacmyiy + Qam 
Qy == Pun + Our) 


Q may be thought of as the velocity potential due to the doublet source of the 
shot point and the image of its first reflection in the top surface. Carrying out 
the additions we obtain: 


Qa = — A(—1)™e"5™i sin (2rd cos 6/L) exp izr[x sin 6 + z cos @ — Vt]/L, 
Q,, = A(—1)™e"5™i sin (2d cos 6/L) exp i2r[x sin 0 — z cos 0 — Vt]/L. 


(8) 


(9) 


Constant factors have been included in the unevaluated factor A. Q, and Q, will 
both contribute to the disturbance at any one point. Their combined effect is 


o= Qa+t g.., or 
@ = A(—65)™ sin (2rd cos 6/L) sin (24z cos 6/L) exp i2r[x sin @ — Vt|/L, (10) 
where constant factors have again been included in A. 

The waves given by ¢ for m reflections interfere with those of @ for m+1 


reflections. Reinforcement occurs if ¢mn-++-¢m41 is a maximum or when the quantity 
1—ed=1—e€ exp i272H cos 6/L is a maximum. This condition is fulfilled for 


those wave lengths which satisfy the equation 

2r2H cos 0/L = (2n — 1)r (11) 
where m=1, 2, 3... The wave length factor is then 

L = (4H cos @)/(2n — 1). (12) 


Here the integer m determines the harmonic or the normal mode of propagation 
and the dependence of 6 upon / establishes L as a function of t. 

Since 6= —1 under the condition for optimum wave development we may 
write 


ar or 
o, = Ae™ sin (= 3 COS s) sin (= d cos s) exp i2m(x sin 0 — Vt)/L. (13) 


With the aid of Equation 12 velocity potential can be written 


on = Ae™sin [(2n — 1)rz/2H | sin [(2n — 1)rd/2H | exp i2r(x sin @ — Vt)/L.(14) 


This expression represents a progressive wave traveling in the water layer in 
the x direction with dispersive characteristics, since L=f(t), and with a standing 
wave pattern between the surface and bottom which is independent of 6. Nodal 
planes are consequently parallel to the surface. The summation of the velocity 
potentials ¢, for all modes, from »=1 to ©, and over as large a range of mas 
applies, gives the over-all disturbance at any point or instant of time. The 
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factor e” shows the effect of partial refraction into the bottom at each reflection. 

A complete description of the resulting wave form would, of course, require 
a careful analysis of amplitude relationships embodied in Ae”. No record was 
kept of automatic volume gain expansion. A typical figure on volume control 
capacity is 10°. Thus, while high amplitude particle motion arrives shortly after 
a shot, amplifier gain is throttled to the extent that trace motion on the seismo- 
gram may have an amplitude comparable with that produced by only 1o~ times 
as much particle motion with controls completely released. Since true amplitude 
ratios can be identified on observations only within limits as broad as this, and 
since elastic constants of the bottom were not measured, it was not considered 
worthwhile at this stage to identify expected amplitudes by theoretical expres- 
sions. 


Dispersion Characteristics 


The rate of progress of a particular phase is identifiable by the relation 


2m 
d — (x sin 8 — V2) 
OL 





= 0, 15) 
a (15 

Expressing L and @ in terms of x and Vi and differentiating (v?/? —x*)'/? leads to 
the phase velocity expression 
c = V/sin @. (16) 


Group velocity measures the rate of progress of the wave group associated with 
each wave length. From simple geometry in this case, or from the general 
formula 


U=c + Kdc/dK (17) 
where K =(27/L) sin 6, one can identify group velocity as 
U = V sin @. (18) ° 
Since sin 0=x/Vt, we also have 
U = aft. (19) 
For each n or each normal mode of propagation frequency 
f = V/L = (2n — 1)V/4H cos 8. (20) 


For reasons explained earlier we limit our discussion to values of @ in the range 
6.2020 where 8, is the critical angle of reflection at the bottom, the angle with 
the vertical whose sine is the ratio of water velocity to the velocity of underlying 
material. The corresponding range of frequency is f,=f2fo where /, is the cutoff 
frequency beyond which the earlier theoretical discussions apply (see Pekeris or 
Press and Ewing, 1950) and fo is the cutoff frequency given by (2n—1)V/4H 
for 0=o0. See Equation 2. 
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INTERPRETATION OF THEORY 


An impulsive point source of sound in a water layer generates spherical waves 
which undergo numerous multiple reflections from the surface and a level bot- 


























Fic. 16. Each graph shows wave frequency plotted with respect to arrival time for a particular 
record. Each plotted point represents an average of several actual observations. The smooth curves 
fitted to collections of average points represent the frequencies predicted by theory. 


tom. The disturbance at any distant point is obtained by the superposition of 
all the waves arriving at this point by the different ray paths indicated.” The 


SSRN IR RIE REG aE pe 


2 In fact a spherical wave can be represented by a sum of component plane waves with appropriate 
amplitude factors (Press and Ewing, 1950; Officer, 1951). Our elementary discussion can in a’ sense 
be considered an analysis of these plane waves leaving out the amplitude factor. 
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result of the combination is a single disturbance which travels outward from the 
source and vertically has a standing wave pattern of motion from the surface to 
the bottom. 

The initial disturbance generated by an underwater explosion can be ana- 
lyzed into a broad continuous frequency spectrum. However, only particular 
frequencies are reinforced according to normal mode propagation theory, and 
these are identified by Equation 20. Significantly, it is members of this set of 
frequencies that are observed. See Figure 16. Because these frequencies are a 
function of the incident angle 6, and in turn time /, for any particular x, they can 
be assigned characteristic velocities. The concepts of phase and group velocities 
were then introduced to distinguish between the velocity with which a constant 
phase proceeds and the velocity with which a given frequency proceeds forward. 

_A description of the system of waves associated with values of @ greater than 
the critical angle has been given elsewhere (Press and Ewing, 1950), where it 
was shown that the last waves of this group arrive at a time corresponding to 
propagation with a velocity of about 0.6 V. Following this group an infinitely 
long train of waves arrive associated with normal mode propagation for values 
of @ smaller than the critical angle. If one focuses attention on any given mode 
(by means of filters for example) the frequency of the waves decreases slowly with 
time, the rate of decrease becoming more and more gradual as time increases 
and as the cutoff frequency for 0=o is approached. Phase and group velocity, 
however, are quite sensitive, approaching infinity and zero respectively as 0 
approaches o. 

These waves are attenuated due to partial refraction at each reflection from 
the bottom. For a given shot detector distance, the number of reflections m in- 
creases as 6 decreases, hence attenuation due to reflection increases as travel time 
increases. For a given travel time and distance frequencies of the various modes 
have the ratios 1:3:5:7: ---. To a large extent the resultant wave pattern 
actually recorded is determined by interference between the various modes ad- 
mitted within the filter pass bands used. 

The pressure or velocity field generated by an underwater explosion can be 
discussed most simply in terms of the pressure fields of the component normal 
modes. Expressions for pressure or vertical particle velocity w can be derived 
from Equation (14) and the relationships 
a ge ae 
? . dt 0z 
Thus the vertical variation of pressure and particle velocity are respectively 
given by the factors 


w(CS)q] ~ «(Cd 


nm = 1, 2,3,°°*. (21) 
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The following conclusions are evident from these equations: 

(1) The surface =o is a nodal plane in pressure or an antinodal plane in 
vertical particle velocity for all frequencies and modes. 

(2) The bottom z= is an antinodal plane in pressure and a nodal plane in 
vertical particle velocity for all frequencies in all modes. 

(3) There are n—1 nodal planes in pressure for all frequencies and modes 
beneath the surface located at depths given by 

z 2 4 6 2% — 2 


Stee ’ ’ greets lett s 
H 2n — 1 2n— 1 2an— I 2n— I 











n= I, 2,°°° 


(4) There are m nodal planes in vertical particle velocity located beneath 
the surface at depths given by ~ 
Z I 3 5 an — 1 


<a ’ ’ ss Bet 
H 2n — I 2n— 1 2n — 1 an— 1 








a he 


The above equations would indicate that the response to these waves of a 
hydrophone sensitive to pressure changes and a geophone sensitive to vertical 
particle velocity varies with depth because of the vertical standing wave pattern. 
Since the ideal location of a detector for elimination of these waves is at a nodal 
plane, the preceding criteria may be used to determine theoretically the best 
depth for locating detectors. A similar discussion may be given for charge loca- 
tion. 


MISCELLANEOUS ASPECTS 


It must be obvious that the preceding remarks about the observed multiple 
reflection phenomenon are essentially introductory. Vital aspects of the problem 
have not been adequately treated. Of these, amplitude has already been men- 
tioned. 

A study of reflections from plates of finite thickness might disclose the 
possibility of a more detailed description of bottom conditions than given. 
A. B. Wood in A Textbook of Sound indicates thickness of the reflector is quite 
an important factor. Certainly a complete discussion of the problem would in- 
clude a consideration of possible reverberation in shallow solid or semi-solid layers 
below the water. 

Effect of corrugations along the bottom or top of the water layer has been 
entirely neglected. Nevertheless, a complete understanding of reinforcement 
phenomena would require more attention to this point. 

In localities where the bottom is rugged no periodicity was observed. Even 
gently sloping bottoms must radically modify the described dispersion pattern 
or prevent the development of resonance altogether. Poulter’s interesting analysis 
of slope determination from multiple reflection behavior appears to offer a promis- 
ing approach to a study of the effect of bottom relief. The fact that parallel 
plates comprise only one example of electromagnetic wave guides or cavity 
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resonators suggests that seismic counterparts may not be restricted to the simple 
example here described. 

The assumption that velocity of propagation is constant appears to be 
justified in our example by the close fit between observations and theoretical 
predictions. If dimensions of the region considered were radically altered, varia- 
tions in propagation velocity might be pertinent factors. 

Finally, this most excellent example of a seismic ae guide suggests that 
there may be others less perfectly developed. Numerous instances of periodicity 
on reflection records have been experienced in land shooting. Some may bear an 
interesting relationship to the reinforcement example herein discussed. 

Establishment of a sound wave guide requires, in general, that 

(1) A source of sound energy provide essentially true continuous wave motion 

or pseudo-continuous wave motion through pulse reverberation. 

(2) Boundary geometry, transmission characteristics, and reflection charac- 

teristics be favorable to the development of particular reinforcement 
waves progressing away from the shot. 


SUMMARY 


Among the different types of waves observed in seismic measurements over 
water we can include the normal mode propagation produced by constructive 
interference as described in various papers by Ewing, Pekeris, and Press. A 
particular example of this type of wave develops relatively close to the source 
of energy provided bottom conditions are favorable. Seismograms from a par- 
ticular prospect show this phenomenon so pronounced that all normal reflections 
from underlying beds are completely lost. 

At the shot point the resonant condition which develops is analogous to that 
which develops in an organ pipe closed at one end. The geometry of spherical 
waves multiply reflected at the boundaries of the water layer provides for a 
spreading of the resonant zone out from the shot point. This progressive growth 
is accompanied by dispersion which takes the form of gradual changes in domi- 
nant frequencies. 

Geometry and terminology of electromagnetic wave guides applies very 
nicely to this particular instance of normal mode propagation. The body of shal- 
low water herein described might therefore be termed a leaking wave guide. 
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SEISMIC RECORD SECTIONS* 


H. R. PRESCOTTT 


ABSTRACT 


A method is described whereby seismic records may be secured which have sufficient dimensional 
stability that a composite reprint of several records, from a continuous profile, when placed in proper 
sequence, yields a cross section of seismic arrivals. A method is described for entering corrected time 
scales on the seismic record section itself. Double recording is provided, with traces in transposed 
positions but alike in phase relation, which makes possible the construction of closed traverse sec- 
tions. A few examples of the record section are shown and some uses suggested. 


INTRODUCTION 


There are times when the seismic interpreter strongly desires to see a con- 
tinuous traverse of seismic records in composite form for critical study. More 
than this, the routine use of a scheme which produces a composite seismic arrival 
section can be made both desirable and convenient for the interpreter. One 
method of producing the record sections will be given. 


APPARATUS FOR SECURING THE FIELD RECORD 


Essentially, the seismic recording equipment is conventional in design with 
a few exceptions. The recording camera produces records of very uniform dimen- 
sion along the time axis. This is accomplished by use of a synchronous motor 
drive in the camera, electronically tied to the timing source for stability. To 
allow high-quality reprints, transparent film is used. To provide a final composite 
section having a minimum width for convenience of routine study and final re- 
view, the width of the seismic record itself is held to less than that of conven- 
tional records. A width of two inches is currently in use. It is desirable that the 
quality of the optical traces be finer and’ sharper than is necessary on conven- 
tional records for the reason that the width of the record is restricted. 

Although not an actual necessity for production of record sections, the use of 
double recording accomplishes some additional objectives. These objectives will 
be described later. 

Figure 1 shows a field record produced by recording equipment intended for 
record section use including double recording. As currently practiced, the 
original double field record is four inches wide and the separate two inch wide © 
records are later secured by splitting the original record. The traces on the 
halves of Figure 1 are alike in phase relation but are transposed in position, as 
diagrammatically indicated. For example, a specific channel indicated by the 


* Presented at the St. Louis Meeting of the Society, April 23, 1951. Manuscript received by the 
Editor May 14, 1951. 
t Continental Oil Company, Houston, Texas. 
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Fic. 1. Double recording where the traces on the halves are alike in phase relation 
but are transposed in position. 
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numeral 1a produces traces 1b and tc of like polarity; but on one of the two 
halves it appears as the bottom trace, whereas on the other half it appears as 
the upper trace. This feature of transposed position of traces is maintained by 
choice as will be explained later. 

The timing lines on the original field record are restricted to short ‘tips’ 
along the edges of the two inch records as shown by Figure 1. This leaves most 
of the record clear for later application of corrected time scales. 


APPARATUS FOR PROCESSING THE FIELD RECORD 


Although the camera which secures the field record may be quite precise in 
speed, the final record itself will show some dimensional variations between 
equal time intervals. The magnitude of the dimensional instability will depend 
on the material used, on the environmental humidity, and on the age of the 
record following photographic processing. It has been found that cellulose acetate 
film is generally sufficiently stable dimensionally for the purpose. The length 
will vary about one part per thousand for each ten percent change of relative 
humidity. After prolonged aging, the length may reduce nearly one percent. To 
obviate this instability, the following process has been adopted. 

In the field office, the original record is positioned in a stretching frame 
and gently stretched to match carefully in length a set of timing lines etched on 
a zinc plate. Several zinc line etchings are available to provide a range of dimen- 
sional choices from areas of low relative humidity to areas of high humidity in 
order that the record need not be stretched unduly. 

After the record is stretched to accurately match the zinc plate in use, 
cellulose acetate ink is applied by a brayer to the zinc line etching. Next, the 
record is superimposed to position properly the zinc scales for corrected ‘times 
and the film record is pressed against the zinc line etching to receive an ink 
impression. The ink is applied to the back of the film in order that the timing 
scale marks may be removed by a solvent without disturbing the emulsion 
should it become desirable to erase them. Silver or orange-colored ink is used 
to contrast more readily with the black timing line ‘tips’ on the original film. 

Figure 2 shows the process somewhat more in detail. Each of the numerals 
23 and 24 indicates one half of a double record where the ground was occupied 
from one shotpoint to another and where a shot from each shotpoint was re- 
corded. The numerals 21, 22, and 25 indicate in diagrammatic outline the position 
of the zinc line scales. The numeral 28 indicates the position of the instant of 
blast. On record 23 the correction of the computing trace 29 for weathering and 
elevation is 0.025 second. The zinc scale 21 was moved along the time axis until 
the zero mark coincided with a value of 0.025 second after the instant of blast. 
The times read along this scale are corrected values. 

On record 23 the computing trace 30 bears a correction of 0.022 second. The 
zinc scale 25 was moved along the time axis until the zero mark came at 0.022 
second. This scale then produces corrected times for the computing trace 30. 
The uncorrected time scales on the common edges of records 23 and 24, for ex- 
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Fic. 2. The application of corrected time lines is illustrated on this pair of records obtained by 
occupying the surface between two shotpoints where a shot from each shotpoint was recorded. 
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ample, mismatch physically by the difference of the corrections applied, 
(0.025 —0.022 = 0.003 second). 
Similar treatment of the record 24 with zinc scales produces the result shown 
in Figure 2. 
THE RECORD SECTIONS AND THEIR USE 


After the original field record has received the ink impression of the corrected 
timing scale and the original double record has been separated into its halves, 
the records may be assembled in proper sequence to produce the record section 
as illustrated by Figure 3. Generally one day’s work is assembled on a printing 
table and reprints secured, each of which is a single sheet. More than one day’s 
work may be readily joined by pasting separate sheets together. 

By choice, no selection of reflections is made on the original film, but selections 
are marked on the reprints. This permits fresh, unmarked reprints to be secured 
at will for new review efforts. 

As interpretation of an assignment progresses, the reviewer may quickly 
obtain record section lines for study of the significant structural and stratigraphic 
features. He may readily prepare separately, in a single record section piece, 
closed loops for critical study of correlations and misclosures as illustrated by 
Figure 4. The transposed position of traces on the two halves of the original 
record allows assembly of continuous single pieces regardless of the orientation 
of recording equipment on the ground. The double recording illustrated by Figure 
1 cannot be obtained by reverse printing of the original record or by a mirror 
image. To do so would transpose the position of the traces as desired but would 
invert the phase relation. 

Although the double record anticipates the advantageous use of different 
equipment adjustments by choice, it has been interesting to observe that oc- 
casionally an interpreter can advantageously use a transposed record, on which 
there is no difference except transposition of the traces, for evaluation of a critical 
reflection selection. Although the records in this case would carry precisely the 
same traces, a reflection sometimes has a different appearance on the transposed 
record, presumably for the reason that the wavelet is not symmetrical. 

Beyond the simplified type of record section illustrated by Figure 3 and 
further illustrated by Figures 4, 5, 6, 7, 8, and 9, other refinements can be made. 
For example, the arrivals have stepout times which tend to detract from cross- 
section use. Adjacent timing scales are staggered by the difference in correction, 
although this is not usually troublesome. The recording equipment can be 
designed, if desired, to produce the oscillograph type of record in the field and 
variable density traces (or variable area) at the same time. Reprocessing of the 
variable density record can be done by photo-electrical cells viewing variable 
density traces through a slit. The slits can be staggered in position to correct for 
weathering, elevation, and a variable dimension with time to compensate for 
normal stepout time of a selected velocity function. Such a reprocessed record 
in oscillograph trace form would then show the seismic arrivals in true relation- 
ship for all traces. 
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Fic. 3. A continuous straight line record section profile taken in southern Oklahoma 
showing prominent thickening toward the basin within the Pennsylvanian section. 
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Fic. 4. A closed loop record section taken in Gaines County, western Texas. 
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Fic. 5. A continuous straight line record section profile taken in Oklahoma 
showing a prominent anticline. 
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Fic. 6. A continuous straight line record section profile taken in Central Oklahoma showing a 
prominent fault in the Ordovician section not detectable in the Pennsylvanian section. 
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Fic. 7. A continuous straight line record section profile taken in southwestern Texas showing a 
prominent unconformity between the Frio and Vicksburg sections. 
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Fic. 8. A continuous straight line record section profile taken in southern Texas showing 
a fault in the Claiborne and Wilcox sections. 





H. R. PRESCOTT 


Fic. 9. A continuous straight line record section profile taken in about 100’ of water offshore 
Texas across a prominent structure shown in the Miocene section. 
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CONCLUSIONS 


Probably the greatest single objection to record sections will come from 
those accustomed to use of records of conventional dimensions. Experience has 
generally shown that this objection disappears as an individual uses the sections. 
Where this objection persists, the double recording feature could be utilized to 
yield one conventional record and one suitable for record section assembly. The 
design and construction of a suitable constant-speed camera and inking devices 
are considered essential to the method and may be considered by some to be 
disadvantageous. Where records lack character and lateral continuity, the 
method may not clearly demonstrate advantages. 

The record section, however, allows an interpreter to review a relatively 
large amount of data and to screen quickly seismic work to identify the critical 
places. It is a rather interesting experience to start with a record section and, py 
the simple procedure of marking a reflection arrival; have a value of time suitable 
for entry on maps. The record section allows both the interpreter and the ad- 
ministrative staff to secure quick comprehensive impressions of the features 
mapped. The designer of seismic equipment may rapidly screen a great many 
records for clues which are important in suggesting improvements. The feld 
observer can more readily recognize the need of apparatus adjustment of various 
sorts by a glance at the record section than by seeing the separate records as 
they are processed. The repetition from record to record, on a composite section, 
of faulty apparatus adjustment becomes quite evident. 

The confined dimensions of the record section allow use of lines at selected 
positions to be filed with final maps to augment the knowledge of the history of 
the structural development pictured on the separate maps where the maps at 
several levels may fail to give the complete story. Experience with their use for 
both current interpretation and subsequent review has shown the seismic record 
section to be a convenient and useful tool. 
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NEUTRON LOGGING OF DRILL HOLES: THE 
NEUTRON-NEUTRON METHOD* 


C. W. TITTLE,f HENRY FAUL,{ ann CLARK GOODMAN$ 


ABSTRACT 


Experiments were performed to determine the distribution of thermal neutrons and of indium 
resonance neutrons in continuous hydrogenous media and in pipes passing through hydrogenous 
media. Included in the study were water, brine, mixtures of sand ad water, and mixtures of sand 
and_brine. Experiments in a continuous typical barite drilling mud showed that the neutron distribu- 
tions were essentially the same as in water. Also, from the point of view of these experiments, oil 
and fresh water are nearly identical. 

These experiments show that well fluid (and, by inference, cement) imposes serious limitations 
on the sensitivity and accuracy of the neutron-neutron logging method. The indium resonance 
neutron response (or, in general, the intensity of epithermal neutrons) is a more reliable indicator of 
hydrogen content of the formation than is the thermal ncutron response. The neutron-neutron 
method of chlorine determination was found to be not sensitive enough to be useful with brines of 
the concentrations ordinarily found in reservoirs. 


INTRODUCTION 


The present paper is concerned with the neutron-neutron or u-n method, 
in which the neutrons themselves are detected, after slowing down and diffusing 
in the formation and, in some instances, returning to the hole to be picked up 
by a neutron detector. The field use of this method has only recently been re- 
ported in the literature (Deegan, 1950), but no detailed information on the pos- 
sible advantages of the method has been available. The response is similar to 
that of the neutron-gamma (or u-y) method (Faul and Tittle, 1951). By virtue 
of the greater sensitivity of the m-x method to the ability of the formation to 
absorb slow neutrons, it is possible under certain favorable circumstances to 
distinguish brine from oil or fresh water. Both the n-n and n-y methods detect 
hydrogen and hence can locate porous zones. Since porous formations are nearly 
always fluid laden, and since the fluids encountered are always hydrogenous 
(oil, water, brine, gas), the hydrogen log can be interpreted in terms of porosity 
(Bush and Mardock, 1950). Gypsum, which contains hydrogen in solid combina- 
tion, might erroneously be interpreted as a porous formation. Gas might also 
be misinterpreted as petroleum or water in a less porous formation, if hydrogen 
content alone is employed as a measure of porosity. The methods discussed in 
the present paper in conjunction with the density log (Faul and Tittle, 1951), 
may lead to more quantitative determinations of porosity. | 
There is a serious limit to the ability of the m-n method to detect brine. We 
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find that casing and well fluid are important factors. In most instances they ob- 
scure the picture to an extent that makes the method impracticable. Attention 
has been directed toward these fundamental limitations. Laboratory experiments 
designed to show what can and what cannot be done have been performed under 
controlled conditions. An attempt has been made to develop a quantitative 
method of interpretation. 

Cement outside the casing, since it is hydrogenous (typically 55 percent 
water by volume), also imposes limitations and hinders interpretation, but a 
given thickness of cement would be equivalent to a somewhat smaller thickness 
of well fluid (typically 80 percent water or equivalent by volume). We have 
performed no experiments with cement, but the experiments on well fluid can 
be applied to problems involving cement by arguing in terms of equivalent 
thicknesses. 


After Vn 
impact: 
Before impact: ym 


Op as Pn (se 


m Vem ~ Vem 


Vr 


Fic. 1. Classical model of elastic neutron scattering; m=neutron mass, A=mass of struck 
nucleus, Vo=initial velocity of neutron, V,=final velocity of neutron, V,=recoil velocity of struck 
nucleus, ¢=angle of scattering, V-m=velocity of the center of mass. 


GENERAL CONSIDERATIONS 


Let us consider an essentially point source of neutrons surrounded by a mass 
of matter so large that only a negligible fraction of the neutrons can escape; i.e., 
in an effectively infinite medium. Let us further assume, for the time being, that 
the medium consists of a single chemical element of atomic mass A and atomic 
number Z. 

Every neutron emitted by the source will then collide with the nuclei (col- 
lisions between neutrons and electrons can be disregarded) of the medium, and 
the collisions will be of three distinct types: 

A. Elastic Scattering (Figure 1)—The neutron may strike the nucleus much 
like a marble hitting a billiard ball and rebound in such a way that the total 
momentum of the neutron and the nucleus is the same before and after the 
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collision. Some of the kinetic energy of the neutron will be lost to the nucleus 
(recoil energy), and the fractional amount will depend only upon the relative 
masses. The velocity and direction of the center of mass of the system remain 
constant throughout, and the ‘entire process can be described in terms of ele- 
mentary mechanics of elastic particles (Marshak et al., 1947 and 1949). 

B. Inelastic Scattering —Instead of a simple mechanical collision, the neutron 
may impart a greater portion of its kinetic energy to the nucleus, leaving it in an 
excited state. The emitted neutron has considerably less energy than the incident 
neutron, and the excited nucleus returns to ground state usually by the emission 
of one or more gamma rays. 

C. Capture.—In some cases, the impinging neutron forms a compound nucleus 
which either: (a) disintegrates by the emission of charged particles (protons, 
alphas, or fission fragments), (b) settles to a stable ground state by gamma 
ray emission (radiative capture), or (c) decays by 8 radioactivity (with or with- 
out accompanying gamma rays). 

The probability that an atom will undergo any of these processes is generally 
expressed in units of the nuclear area (1 barn=10~* cm?) or cross section o for a 
particular reaction. The cross sections are different for different isotopes of the 
same element and show a marked dependence on the energy of the impinging 
neutron. Some isotopes have large capture cross sections at certain neutron 
energies, known as resonance energies (Figure 2). Most of the pronounced capture 
resonance peaks occur at neutron energies between 1 and 50 ev (Goldsmith 
et al., 1947; Adair, 1950). 

Let us assume that the cross sections for inelastic scattering and capture are 
small in our medium, and for the moment, consider only elastic collisions. The 
source emits neutrons with a velocity vo. The velocity of the center of mass of 
the system, consisting of a neutron and the target nucleus, is then %m=mv0 
/(m—A), assuming the nucleus to be stationary with respect to the source. But 
m is practically unity in atomic mass units, and we may write v¢m=v/(1—A). 
If we place the origin of our co-ordinate system at the center of mass of the two 
particles, the velocities will be vo/(1—A) for the nucleus, and Avo/(1—A) for 
the neutron. In the center of mass system, the magnitudes of these velocities 
will remain the same before and after the collision, although the directions will 
change. (The following discussion is largely based on the lectures of E. Fermi, 
1949.) 

We may then compute the final velocity », of the neutron with respect to the 
source by adding the velocity vector of the center of mass to the velocity vector 
of the neutron in the center of mass system. It follows that 


E = 0,2/2 = v9?(A? — 2A cos @ — 1)/2(1 — A)? (1) 
where E is the final energy of the neutron, @ is the angle of scattering in the 


center of mass system, and the neutron mass m=1. The corresponding angle 
¢ in the laboratory system (with respect to the source) is given by 


cos ¢ = (A cos 6 + 1)/(A? + 2A cos 6 + 1)!/2 (2) 
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from which it can be shown that the average value 
cos = 2/3A (3) 


showing that the neutrons tend to be scattered in the forward direction. This 
forward scattering is greatest for collisions with the lightest nuclei (hydrogen, 


A=1). 
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Fic. 2. Total neutron cross sections of indium and cadmium (after Goldsmith, Ibser, and Feld) 


The neutron energy loss per collision can be computed from 
E,/E = (A + 1)?/(A? + 2A cos @+ 1) (4) 
which has the mean logarithmic value 
(A — nn A+ r 
(5) 


Ey 
¢-[mZ| =I- 
E av 2A 7 ne 


| For hydrogen, the geometric mean energy E (Bethe, 1937) of the scattered 








neutron is given by 
(E/Eo)x = 1/e 











630 C. W. TITTLE, HENRY FAUL AND CLARK GOODMAN 





where e is the base of the Napierian logarithms (2.718 -- -). Hence, it is seen 
that on the geometric average a neutron loses more than half its energy in a 
collision with a proton, but it loses only a very small fraction of its energy in an 
elastic collision with a nucleus of a heavy atom. Consequently, for a given 
number of collisions, hydrogenous material is by far the most effective as a slowing- 
down medium for neutrons. In addition, the elastic scattering cross section of 
hydrogen is fairly large: 21 barns for neutron energies between 1 ev and 10 Kev. 
Only about twenty collisions with hydrogen are required to slow the fastest 
neutron from a Ra-a-Be source (approximately 13 Mev) down to thermal 
energies (roughly 0.025 ev). 

As the neutrons approach thermal energies, the average energy loss per 
collision becomes smaller and smaller. A thermal neutron in equilibrium with its 
surroundings neither gains nor loses energy in an average collision, due to the 
thermal motion of the atoms in the material. On the whole, thermal neutrons in 
a medium behave somewhat like the molecules of one gas in another, and the 
diffusion of neutrons can be calculated by using much the same reasoning as 
applies to the diffusion of gases. 

While the neutrons are being slowed down, we can generally neglect capture, 
especially if the substance is hydrogenous, for the capture cross sections of most 
nuclei for fast and intermediate energy neutrons are small. As soon as the neu- 
trons approach thermal energies, however, capture becomes important. The 
capture cross sections of common terrestrial elements (hydrogen, oxygen, carbon, 
silicon, aluminum, etc.) are relatively small, with the notable exception of 
chlorine which has a relatively high thermal capture cross section (¢,= 31.5 
barns). This difference is of importance in neutron logging. (Boron and cadmium 
have much higher thermal capture cross sections (710 and 2,410 barns, respec- 
tively), but they are relatively rare in common rocks. These and other capture 
cross section data are taken from Ross and Story (1949), except that of chlorine, 
which is from Hibdon and Muehlhause (1950). 

So far, we have assumed that our medium consists of only one element, but 
we have shown that the slowing-down power of hydrogen is far greater than 
that of any other common element. Consequently, if we are dealing with a mixed 
medium containing an appreciable concentration of hydrogen, the heavier nuclei 
contribute very little to the actual slowing-down process although they do 
materially affect the spatial distribution of the neutrons through their scattering 
properties. 

In the experiments described below, neutron distributions in continuous 
media were investigated, and the effect of holes on the neutron distributions was 
studied. 


THE EXPERIMENTS 


The experiments in the present work on the n-n log have employed a Ra-Be 
source. The source contained 1.03 grams of radium and was constructed by the 
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modern pressed pellet process. It was handled by a special electromagnet on a 
10-foot pole. 

The experiments have been of two general types, which might be called ex- 
periments in continuous media and pipe experiments, respectively. Since the 
scattered neutrons are detected in the n-x method, the spatial distributions of 
neutrons in continuous media and in the pipes running through continuous media 
were measured. The distributions of slow neutrons (thermal and epithermal), 
resulting from the slowing down and diffusion of the fast primary neutrons from 
the source, are generally of the greatest interest. 

Experiments in the following continuous hydrogenous media were performed 
to aid the interpretation of the pipe experiments: 


a. water 

b. saturated brine 

c. sand-water mixtures 

d. a typical barite drilling mud. 


The neutron distributions in the drilling mud were practically the same as in 
water, hence to avoid complications water was used in the pipe experiments to 
simulate well fluid. 

Various sizes of iron pipe were used as casing. Aluminum pipe was used to 
simulate uncased holes, since aluminum has low neutron scattering and capture 
cross sections (1.43 and 0.22 barns respectively, for thermal neutrons). Specific 
items which were investigated include: 


a. effect of more than one string of casing 

b. masking effect of well fluid and casing 

c. eccentricity of logging instrument position 

d. effect of an empty hole, or of a long essentially “empty” logging instru- 
ment, in distending the neutron distribution axially. 


The Slow Neutron Spatial Distributions: Thermal and Epithermal Neutrons.— 
In the experiments, slow neutrons of two different energy ranges have been 
detected and their spatial distributions determined: (a) thermal neutrons, or 
neutrons substantially in thermal equilibrium with the medium and characterized 
by a Maxwell-Boltzmann number-speed distribution law with a most probable 
speed corresponding to an energy kT, where k& is Boltzmann’s constant and T 
the absolute temperature of the medium, and (b) epithermal neutrons, or 
neutrons just above thermal energies and generally ranging from about 0.2 ev 
upward. The particular epithermal energy employed in the present work was 
1.44 ev, the energy of the principal neutron resonance in indium. 

The thermal neutron distribution depends on the temperature of the medium. 
Since subsurface formations commonly have elevated temperatures of the order 
of 150 or 200°F., temperature needs to be considered in applying the results of 
laboratory experiments, ordinarily performed at about 20°C. (68°F.), to sub- 
surface conditions. 
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EXPERIMENTAL RESULTS 


Details of experimental arrangements and correction of data appear in 
Appendices A and B, respectively. The cadmium correction factor for indium 
resonance neutrons and the Bothe correction factor for thermal neutrons have 
been applied to the data. 

Experiments in continuous media.—Thermal and epithermal spatial neutron 
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Fic. 3. Distribution of indium resonance and thermal neutrons around the Ra-Be source in 
water and in drilling mud (61.5 percent water, 35 percent barite, and 3.5 percent bentonite by weight, 
gravity 1.38). Note: Throughout this paper, neutron density is reported as the saturation counting 
rate of the front plus the back of a thin indium foil (see Appendix for details). The corresponding 
ordinates are marked “counts/minute.” The distances from the source, plotted on the abscissae, 
are marked Z when measured longitudinally along a pipe and r when measured radially outward 
through the material under study. 
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distributions in water, a typical barite drilling mud, a.concentrated brine (NaCl 
solution), a mixture of sand and water, and a mixture of sand, water, and air 
are presented graphically in Figures 3-5. Details are given in the captions. 
Figtre 3 shows the distribution of indium resonance neutrons and thermal 
neutrons around the Ra-Be source in water and establishes that the neutron 
distributions in drilling mud and in water are nearly the same, i.e., water can 
be used to simulate well fluid. Figure 4 indicates that the indium resonance 
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Fic. 4. Neutron distributions in concentrated brine (98 percent saturated NaCl solution). 
The distributions in water are shown for comparison (dashed lines). 


neutron distributions in water and concentrated brine (98 percent saturated) 
are nearly the same, but there is a pronounced difference between the thermal 
distributions, principally due to the large thermal capture cross section of chlo- 
rine. As expected, Figure 5 shows that if the hydrogen concentration is reduced 
by the addition of sand to water, the neutron distributions fall off less rapidly, 
since hydrogen dominates in the slowing-down and diffusion processes. The 
effect is emphasized when most of the water is drained from the sand. This 
mixture represents the only instance in which the medium under study was not 
large enough to be considered effectively infinite in extent. The finite size tends 
to make the distributions fall off somewhat too rapidly at large distances from 
the source. 

Experiments with pipes through water—A number of experiments were made 
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Fic. 5. Neutron distributions in sand of 39 percent porosity, saturated with water (A), and 
containing 9.5 percent water and 29.5 percent air by volume (B). 


TABLE I 
DIMENSIONS OF IRON AND ALUMINUM PIPES USED IN THE EXPERIMENTS 











Nominal Inside Wal Outside 

size (IPS) diameter thickness diameter 
4-inch iron 4.03 in. 0.24 in. 4.50 in. 
s5-inch iron 5.05 0.26 5.50 
— iron 5.50 0.25 6.00 
6-inch iron 6.07 0.28 6.68 
8-inch iron 7.98 0.32 8.63 
4-inch aluminum Same as 4-inch iron (used as instrument) 





6-inch aluminum Same as 6-inch iron 














NEUTRON LOGGING OF DRILL HOLES 635 


in pipes (iron and aluminum) of various dimensions summarized in Table I. 

Figures 6-9 show the thermal and resonance neutron distributions along the 
axes of these pipes traversing water. Figure 6 indicates that the iron in the 4- 
inch pipe depresses the thermal neutron density as compared with the “‘uncased” 
hole, which is simulated by an aluminum pipe of the same size. Also shown is 
the effect of an additional 5-inch string of iron casing. The resonance distributions 
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Fic. 6. Neutron distributions in cased and uncased holes through water. The uncased hole is 
simulated by a 4-inch aluminum pipe, which is practically transparent to neutrons. 


were not affected very much by the presence of the iron. The space between the 
two pipes was filled with water. Figure 7 shows how the neutron distributions 
inside a 4-inch aluminum pipe, representing a logging instrument (outside 
diameter 4.5 inches), depend on the size of the casing. The space between casing 
and instrument was filled with water in each instance. Since the medium beyond 
the pipe was also water, it was anticipated that the effect of the iron would be the 
least for the largest pipe. Figure 8 shows how the neutron distributions in 
aluminum pipes through water depend on the pipe size. The neutron distributions 
in the larger hole fall off more slowly with distance. Figure 9 shows the same 
effect with iron pipes. The holes contained only the neutron source and detecting 
foils with their holders. 

Experiments with pipes through concentrated brine. 





Figures 10-14 present 













































































ee 
ie ~~ he, 
<3. Mo 
jor PSS 
= SSAC es aN! ile 
z “= fiass st | B 
E aa 
> SS —_ 
r= On 
310° DNS 
" oa dk 
. 
-— A uncased °so SS 
9 S 
B. cased, 8-inch rey 
Pas 
C cased, 6-inch A a, : 
lo2+— 1D cased, 5-inch S~ 
| | | 
0 10 “20 30 40 50 60 70 
L,cm— 


Fic. 7. Neutron distributions in cased holes of various diameters in water. Distributions in 
4-inch uncased hole are shown by dashed lines for comparison. 
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Fic. 8. Neutron distributions in uncased holes of two different diameters in water. 
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Fic. 9. Neutron distributions in essentially empty iron pipes of four different diameters in water. 
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Fic. 10. Neutron distributions in cased and uncased holes through 98 percent saturated brine. The 
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distribution in an uncased hole through water is shown for comparison (dashed line). 
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Fic. 11. Neutron distributions in 6-inch cased holes through water and saturated: brine with 
logging instrument in coaxial position (dashed line) and extreme eccentric position (solid line). Well 
fluid is water. 
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Fic. 12. Neutron distributions in 8-inch cased holes through water and saturated brine with logging 
instrument in coaxial position (dashed line) and extreme eccentric position (solid line). 
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results of various experiments with pipes in a 98 percent saturated sodium chlo- 
ride solution. Figure 10 shows the neutron distributions along the axes of the 
pipes. The 4-inch aluminum pipe served as the logging instrument in all instances 
except the one marked “Cased, two strings (4 and 5 inches),” in which both 
pipes were iron. The space between pipes was filled with water, serving as well 
fluid. The masking effect of well fluid is clearly shown. The larger the layer of 
well fluid, the more nearly is the distribution like that of the uncased hole in 
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Fic. 13. Neutron distributions in saturated brine radially outward from cased holes. 
The neutron source is on the axis of the cased hole. 


water. With the 13-inch layer of well fluid in the 8-inch cased hole, the well fluid 
alone is being logged. 

Eccentricity of the logging instrument.—Figures 11 and 12 illustrate the effect 
of eccentricity of the logging instrument in 6-inch and 8-inch cased holes. Water 
served as well fluid between probe and casing; the medium beyond the casing 
was water or concentrated brine. The difference in thermal distributions in the 
instrument for water and concentrated brine beyond the casing was increased 
slightly by the eccentricity but not enough to be decisive. The masking effect of 
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well fluid is still severe. It is seen that variations in radial position of the instru- 
ment in the well will produce fluctuations in response, so that it is desirable to 
use a positioning device to hold the probe in either a coaxial or an eccentric 
position. 

Radial neutron distributions around a well.—Figure 13 shows the neutron 
distributions radially outward into concentrated brine for 5-inch and 8-inch 
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Fic. 14. Contours of equal thermal neutron flux in a cased 8-inch well through 
concentrated brine (idealized). 


cased holes. The aluminum instrument was essentially empty, as in all the 
experiments involving it, and the space between instrument and iron pipe was 
again filled with water. There is little difference between the two conditions. 
Figure 14 is a sketch of contours of equal thermal neutron flux in and around the 
cased 8-inch well. The effect of the hollow inner pipe in stretching the neutron 
distribution is clearly shown. Since measurements were made only axially and 
radially, positions of the curves at intermediate points are only suggestive. 
Many more experiments of this type, in which the neutron distributions are 
measured around wells in various media of geological interest, will have to be 
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performed before a real understanding of the neutron logging processes is 
attained. 

Experiments with pipes through brines of different concentrations.—Figure 15 
shows the neutron distributions in uncased holes (lined with aluminum pipe 
4.5 inches in outside diameter) through brines of various concentrations. Figure 
16 reveals similar variations for 5-inch cased holes, in which water filled the 
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Fic. 15. Neutron distributions in an uncased hole through brines of different concentrations. 


annular space between the 4.s-inch aluminum instrument and the 5-inch iron 
pipe. 

Experiments with pipes through a mixture of sand and water—The silica sand 
used in these experiments was obtained from the Ottawa Silica Co., Ottawa, Ill. 
The mixture obtained by pouring the sand into the experimental tank while 
spraying it with water had a pore space which was 39 percent of the total volume, 
considerably more porous than typical oil reservoirs (10-30 percent). Figure 17 
shows the neutron distributions along the axes of empty 4-, 5-, and 8-inch iron 
pipes through the sand-water mixture. One of these experiments showed that 
introduction of the aluminum logging instrument pipe had a negligible effect on 
the neutron distributions. Figure 18 shows the axial neutron distributions in 
the 4.5-inch O.D. logging instrument with no casing and with 5-inch and 8-inch 
casing. In the last two instances water filled the space between instrument and 


casing. 





Fic, 16. 


Fic. 17. Neutron distributions in empty 4-, 5-, and 8-inch pipes through a silica sand of 
39 percent porosity saturated with water. 
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Fic. 18. Neutron distributions in cased and uncased holes in the saturated sand-water mixture. 


Experiments with pipes through a mixlure of sand, waler and air—The mixture 
was obtained merely by draining most of the water from the mixture of sand 
and water just described. The porosity was unchanged, .i.e., 39 percent. The 
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Fic. 19. Neutron distributions in empty 4-, 5-, and 8-inch iron pipes through a silica sand of 39 per- 
cent porosity containing 9.5 percent water and 29.5 percent air by volume (drained sand). 











644 C. W. TITTLE, HENRY FAUL AND CLARK GOODMAN 


water and air contents were 9.5 and 29.5 percent by volume, respectively. Figure 
19 shows the neutron distributions along the axes of empty 4-, 5-, and 8-inch 
iron pipes through the sand-water-air mixture. Figure 20 shows the axial neutron 
distributions in the 4.5-inch O.D. aluminum logging instrument with no casing 
and with s-inch and 8-inch casing, water filling the space between instrument 
and casing in the last two instances. 
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Fic. 20. Neutron distributions in cased and uncased holes through the 
sand-water-air mixture (drained sand). 


Experiments with pipes through sand-brine mixtures of different chloride 
concentrations—The brines used in the sand-brine experiments were dilute, 
containing 1.18, 2.32, 3.69 and 5.11 percent Cl by weight. The most concentrated 
brine was 31.5 percent saturated. The porosity was the same as in the sand-water 
mixture, 39 percent. The pore space was completely filled with brine. Figure 21 
shows the neutron distributions along the axis of a 4-inch iron pipe passing 
through the sand-brine mixtures. 


INTERPRETATION OF RESULTS . 


The quantitative relations between u-n log response and hydrogen content 
and the possibility of detecting brine by the thermal neutron n-n method are 
covered mostly by the data of Figures 15-21. 

Neutron-neutron log response vs. hydrogen content of formation.—A strictly 
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functional relation between the n-n log response and the hydrogen content of the 
formation is not to be expected because of the influence of the non-hydrogenous 
components on the neutron distributions. This is especially true of the thermal 
neutron distribution, in which the macroscopic thermal capture cross section 
plays a strong part as illustrated by the data obtained with brine. Properties 
related to the rate at which the distributions (thermal or epithermal) decrease 
with distance along the pipe reflect the hydrogen content more accurately than 
the thermal neutron response itself, because the latter depends on the thermal 
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Fic. 21. Neutron distributions in empty 4-inch iron pipe through sand of 39 percent porosity 
saturated with brines of different concentrations. The brines contained 1.18, 2.32, 3.69, and 5.11 
percent Cl by weight. 


neutron density which in turn varies with capture cross section. The epithermal 
neutron response, however, is a fairly good indication of hydrogen content, since 
the epithermal neutron density depends on the neutron slowing-down properties 
of the medium and hydrogen tends to dominate the slowing-down process. Hydro- 
gen content will be expressed as percent water by volume in the discussion to 
follow. Capture will in general be negligible at epithermal energies in cases of 
geological interest. 

Equivalence of cement and well fluid——A typical cement contains 55 percent 
water by volume. Some cements contain boron in the form of borax or boric 
acid which serves as a retarding agent, and due to the high capture cross section 
of boron (710 barns for thermal neutrons), such cements will be more effective 
in masking the formation. Roughly, 1.5 inches of boron-free cement are equiva- 
lent to 1.0 inch of well fluid. 
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Terminology.—Since aluminum is essentially transparent to neutrons, the 
holes produced by aluminum pipes will be referred to as uncased, with a diameter 
equal to the outer diameter of the pipe, whereas the iron pipes will always be 
designated by their imner diameter. In the following discussion, the terms ‘“un- 
cased,” ‘“‘cased,”’ and “iron pipe” have the meanings described below. 

Uncased: The aluminum pipe serves as logging probe with the formation in’ 


direct contact. 
Cased: The 43-inch aluminum logging probe is located coaxially in an iron 
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Fic. 22. Thermal decay length Zg (th) as a function of iron pipe size and 
percent water in formation (no well fluid). 


casing of specified inner diameter (such as “‘5-inch cased”’). In all instances where 
the term “‘cased”’ is used in the discussion, it is implied that well fluid is present 
in the space between probe and casing. Thicknesses of well fluid are as follows: 


5-inch cased: 4 inch 
8-inch cased: 12 inch. 


Iron Pipe: The term “‘iron pipe” will be used to designate holes lined with 
iron but containing no well fluid. These holes were essentially empty. Hole size 
designation means inner diameter of pipe. 


Thermal decay length as a measure of hydrogen content.—The decay length 
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La(th) of the thermal neutron flux along the axis of the pipe, determining the 
rate at which the flux falls off in the Ce~“/“¢ representation of the distribution, 
is shown in Figure 22 as a function of iron pipe diameter and percent water in 
the formation. The decay length increases with increasing hole diameter and 
with decreasing hydrogen content of the formation. In the region of interest in 
o'l production—1o0 to 30 percent water content, or the practical equivalent of 
10 to 30 percent porosity whether filled with oil or water—L,(th) is a fairly sensi- 
tive indicator of hydrogen content for a given hole size, but it should be noted 
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Fic. 23. Thermal decay length Za (th) as a function of percent water in formation 
showing the effect of well fluid between probe and casing. 


that conditions are comparatively ideal because of the absence of well fluid or 
cement. Figure 23 shows the effect of well fluid in reducing the sensitivity of 
the method. Even }-inch of well fluid in the 5-inch cased hole seriously affects 
the sensitivity, and 13 inches of well fluid masks the effect of the formation behind 
the casing. The desirability of displacing as much of the well fluid as possible is 
thus demonstrated, but in a cased hole with a fairly thick layer of cement beyond 
the casing—typically 1 to 2 inches—the method would be impractical for quanti- 
tative work. 

Resonance (epithermal) decay length as a measure of hydrogen content.— Figures 
24 and 25 show that there is little difference between the decay lengths of the 
thermal and resonance neutron distributions as measures of hydrogen content. 
The masking effect of well fluid is much the same. Another parameter which 
has the same significance as the decay length is the ratio of two detector readings 
at different distances from the source. 
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Fic. 24. Resonance decay length Lg (res) as a function of pipe size and 
percent water in formation (no well fluid). 
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Fic. 25. Resonance decay length Lg (res) as a function of percent water in formation, 
showing effect of well fluid between probe and casing. 
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Resonance response as a measure of hydrogen content——The response Nres of a 
long detector sensitive to the indium resonance neutrons and extending along 
the axis of the pipe at points 30 to 60 cm from the source was obtained by inte- 
grating the measured distribution over this range. The length and position of 
this equivalent detector were chosen arbitrarily for the purpose of discussion. 
In most instances the detector lies within the exponential part of the distribution. 
Figure 26 shows the resonance response of the long detector as a function of iron 
pipe diameter and percent water in the formation. The flattening of the curve 
for the 8-inch pipe is a result of the particular detector position chosen; at larger 
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Fic. 26. Resonance response of a long detector extending from 30 to 60 cm from source, 
as a function of iron pipe size and percent wate~. 


distances the effect would not be so pronounced. Figure 27 demonstrates the 
effect of well fluid on the sensitivity of the method. The 3-inch layer of well 
fluid has little effect, as shown by the parallel curves for the 43-inch uncased and 
5-inch cased holes. In the 8-inch cased hole the responses differ by 20 percent for 
10 and 25 percent water content, a considerable improvement over the sensi- 
tivity of the decay-length method. 

The effect of moving the resonance detector is illustrated by Figures 28 and 
29, which represent the response of a point detector at 60 cm from the source. 
Sensitivity, as indicated by the slope of the curves, is increased, and the flattening 
of the curve for the 8-inch iron pipe is almost eliminated. Thus with proper de- 
tector positioning, the response of the resonance detector would be a fairly 
reliable indication of hydrogen content of the formation. 

Ratio of responses of thermal and resonance detectors——Figures 30 and 31 
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Fic. 27. Resonance response of a long detector extending from 30 to 60 cm from the source, 


showing the effect of well fluid between probe and casing. 
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FiG. 23. Resonance activities at 60 cm from the source, as a function of pipe size and percent water. 
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show the ratio R of the response of a long thermal detector and that of a long 
resonance detector at the same location, extending 30 to 60 cm from the source. 
Since the resonance response is a measure of the hydrogen content and the 
number of thermal neutrons in the hole depends inversely on the ability of the 
formation to capture thermal neutrons, R is a measure of capture probability. 
The curves of Figures 30 and 31 are for chlorine-free formations. Any medium 
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Fic. 29. Resonance activities at 60 cm from the source, as a function of percent water, 
showing the effect of well fluid between probe and casing. 


in which a strongly absorbing element (such as chlorine) occurs will be charac- 
terized by a smaller thermal to resonance response ratio, and R can therefore 
be used as a parameter that indicates the chlorine content of the formation, on 
the assumption that other strong absorbers are absent. Figure 31 indicates a 
pronounced dependence of the parameter R on well fluid which will seriously 
limit the practicality of the method, as discussed below. 

Chlorine determination by the n-n method —Figure 32 shows how the param- 
eter R depends on chlorine content of the brine in a sand formation of 39 per- 
cent porosity. The hole was lined with iron but contained no well fluid. R is a 
fairly sensitive function of chlorine content. However, in this case the porosity 
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Fic. 30. Ratio of responses of long thermal and resonance detectors (30 to 60 cm from source) 
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as a function of pipe size and percent water in formation. 
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Fic. 31. Ratio of responses of long thermal and resonance detectors (30 to 60 cm from source) 
as a function of percent water, showing the effect of well fluid between probe and casing. 
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is higher than normal and the absence of well fluid or cement makes the prospects 
somewhat optimistic. Figure 33 is a compilation of all the data obtained in the 
present series of experiments concerning the possibility of detecting brine. The 
ratio of parameters R(water)/R(brine) is plotted as a function of percent Cl in 
the brine for various conditions. Enough data were taken to allow estimation of 
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Fic. 32. Ratio of thermal to resonance response of long detector, extending from 30 to 60 cm 
from the source, in an air-filled 4-inch iron pipe surrounded by a sand-brine mixture, as a function of 
percent chloride by weight in the brine. Formation porosity is 39 percent of total volume. Pore space 
is filled with brine. The most concentrated brine shown is 31.5 percent saturated NaC] solution. 


the curve for a 5-inch cased hole through a sand and brine mixture of 39 percent 
porosity. For a 5 percent brine, R(water)/R(brine) =1.3. Since this still repre- 
sents rather idealistic conditions—only 3} inch of well fluid, no cement, brine 
considerably more concentrated than average, porosity twice as high as typical— 
it can be concluded that under field conditions the m-n method will not be a 
satisfactory procedure for determining chlorine content of the formation. Even 
in uncased holes, the layer of well fluid and the uncertainty of the hole size in 
most instances would impose severe restrictions on the method. 


FIELD APPLICATIONS 


Thermal and epithermal neutron detectors—The method of thermal and 
epithermal neutron detection used in these experiments is obviously not suitable 
for field application. Equivalent results can be obtained by use of a proportional 
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counter containing boron, for example the gas BF;, preferably with boron 
enriched in the isotope B!°. Normally such a counter would be sensitive to neu- 
trons of all energies down to and including thermal energies but its sensitivity 
would be greatest at thermal energies, and ordinarily the density of the thermal 
neutrons is greater than that of the higher energy neutrons. Hence the response 
of the counter would primarily depend on the thermal neutron distribution. If 
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Fic. 33. Effect of casing and well fluid on ability. of the n-n log to detect brine. Ratio of 
Nin/Nres ratios as a function of percent chloride in brine in cased and uncased holes. 


the counter is shielded by a layer of cadmium at least 0.02 inch thick, the thermal 
neutrons will be absorbed and will not reach its sensitive volume. The response 
will then be largely a function of the epithermal neutron density. Sensitivity will 
be rather low, but this can be compensated by use of a more intense neutron 
source. These considerations are intended to show only that practical detectors 
of the type needed in the n-n method can be devised; many arrangements are 
possible. 

Field limitations—It has been shown that the thermal neutron response 
depends critically on well fluid (and on cement). It would therefore be a poor 
indicator of hydrogen content of the formation. Furthermore, absorbing elements 
such as chlorine would tend to produce complications. The effect of well fluid, 
casing, cement, and neutron absorption can be reduced by detecting epithermal 
neutrons instead, and the epithermal neutron-neutron logging method should 









































NEUTRON LOGGING OF DRILL HOLES 655 


therefore be a satisfactory indicator of hydrogen content. By virtue of the 
adverse influence of well fluid, casing, cement, and variations in hole diameter, 
the thermal-epithermal method of chlorine detection is faced with too many 
limitations to be practical in the field. 

Comparison of n-n and n-y logging methods——A previous paper of this series 
(Faul and Tittle, 1951) is concerned with the neutron-gamma logging method. 
It is appropriate to say here by way of summary that the methods are of com- 
parable sensitivity if epithermal neutrons are detected in the n-n method, but 
more work will have to be done to determine which method is the more accurate 
indicator of hydrogen content. 


APPENDIX A: EXPERIMENTAL ARRANGEMENTS 


The cadmium difference technique with indium metal foils was employed to 
determine the thermal and indium resonance neutron densities as a function of 
position. After an exposure of known duration, the indium was removed and its 
B-radioactivity ‘(half-life of 54 minutes) measured by a thin-window Geiger- 
Mueller counter. 

Figure 34 illustrates a typical setup of Ra-Be neutron source, foil holders, 
pipes, and media in the experiments. It represents a cased hole with well fluid 
between the aluminum logging probe and the iron pipe. The holder had a funnel 
at the top for ease in inserting the source. The indium detecting foils were placed 
in aluminum or cadmium holders for the thermal + resonance or resonance 
determinations, respectively. The holders were mounted on an aluminum channel 
with slots at fixed intervals relative to the source. 

The pipes were about 4 feet long. The arrangement was centered in a cubical 
4-foot tank containing the material representing the formation under study. In 
the experiments with low hydrogen content the media were not extensive enough 
to be considered as effectively infinite. In such cases, and especially in the case 
of nonhydrogenous media, the amount of material required in neutron diffusion 
experiments is vast. 

The indium foils were circular discs 3 cm in diameter and approximately 
0.003 inch thick. The average surface density was about 58 mg/cm*. The foils 
were intercalibrated for both thermal and resonance neutron sensitivity. A set 
of some 20 foils is required for work of this type. The foils were first placed in the 
holders and the channel assembly inserted in the inner pipe. The source was then 
inserted for a timed interval—usually 54 minutes. The exposure was terminated 
by removing the source. After a delay of 3 minutes to allow for handling and to 
let the 13-second indium activity decay to negligible proportions, the most 
distant foil was counted on an end-window Geiger-Mueller counter shielded with 
lead. The other foils were counted in sequence. Each foil was counted both 
with its front and its back facing the counter window, and after correction for 
background count and for decay of the activity while waiting and counting, the 
activities of front and back were added to obtain a measure of the neutron 
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density (Amaldi and Fermi, 1937). An automatic continuous timing device was 
used to establish a rigid counting schedule. It was found desirable to store the 
indium foils in a cadmium box surrounded by paraffin to prevent stray activation. 
For the experiments in continuous media the same source and foil holders were 
used. 
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Fic. 34. Arrangement of pipes, foil, and source in some of the experiments (not to scale). 


APPENDIX B: CORRECTIONS APPLIED TO INDIUM ACTIVATION DATA 


For accurate work it is necessary to apply two corrections to the data: (a) 
cadmium is not perfectly transparent to neutrons of indium resonance energy 
(1.44 ev), so that for foils of the thickness used in these experiments (58 mg/cm’), 
and for a cadmium thickness of 0.040 inch, it is necessary to multiply the measured 
resonance activity by a factor 1.13 (Kunstadter, 1950); and (b) the introduction 
of a thermal neutron detector depresses the thermal neutron density because it 
is necessarily a neutron absorber and acts as a sink for the diffusing neutrons 
(Bothe, 1943). The essential criterion is the probability that the same neutron 
may cross the foil more than once in its zigzag path through the medium. The 
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foil may thus ‘“‘cast a shadow” upon itself by absorbing neutrons which other- 
wise would have returned to it. These corrections are discussed in detail by 
Tittle (1951). 

The thermal neutron depression (‘‘self-shading”’) correction factor F was 
nearly 1.0 for foils placed along the axis of the pipes in experiments, because the 
innermost pipe was always empty and had a large radius compared to the radius 
of the foil. However, the correction is necessary and was applied where the 
neutron intensity was measured radially outward into the medium and the foils 
were surrounded by the material under study. The different correction factors 
are shown in Table II. The effect is negligible for indium resonance neutrons in 
hydrogenous media. 


TABLE IT 


SELF-SHADING (Bothe’s Selbstbeschattung) THERMAL NEUTRON DEPRESSION 
CORRECTION FAcTORS FOR VARIOUS MEDIA 














, Correction 

Medium Factor F 
Water DE 
Drilling mud, gravity 1.38 (61.5 percent H,O by wt) E.10 
Brine, saturated 1.06 
Sand, water-saturated (39 percent porosity) 1.05 
Sand, drained (39 percent porosity, 9.5 percent water by vol) 1.02 
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THE INVESTIGATION OF EARTH RESISTIVITIES IN THE 
VICINITY OF A DIAMOND DRILL HOLE* 


A. R. CLARK! anp D. J. SALT? 


ABSTRACT 


The method of electrical images has been used to derive an expression for the difference in po- 
tential between two points close to a good conducting sphere in a poor conducting medium. The 
results indicate that the presence of a conducting ore body of commercial size near a drill hole can be 
detected by suitable electrical measurements made in the hole. 

The results of a test of the method in a drill hole near a small conducting sulphide ore body are 


given. 


INTRODUCTION 


Geoelectrical investigations, in which the detection of buried conductors 
depends on the distortion by the conductor of a given potential configuration 
applied at the surface, must be confined to relatively shallow depths. A diamond 
drill hole in an area offers an opportunity to greatly increase the depth which 
may be reached (1). 

If a current electrode C; is placed in a drill hole in a homogeneous isotropic 
medium and four other current electrodes C: are placed on the surface sym- 
metrically spaced around the collar of the hole at large distances from the hole, 
and the resistance to earth of each current electrode made the same, then C; 
will be a point of electrical symmetry as well as space symmetry. The vertical 
line through the point source C, will also be a line of symmetry with respect to 
the surface electrodes. If the material is not homogeneous the vertical line is no 
longer a line of symmetry and when current is sent from the central point elec- 
trode to the surface electrodes in turn, the difference in the drop in potential 
between any two points in the original line of symmetry provides a measure of 
the comparative conductivities of the materials in any of the four directions 
from the bottom of the hole. Thus a lowered potential measurement when the 
current flow is in a given direction indicates material of greater conductivity in 
that direction. 


DISTORTION OF EQUIPOTENTIAL SURFACES BY A CONDUCTING SPHERE 


A rigorous mathematical treatment of the problem encounters severe diffi- 
culties in attempting to apply the boundary conditions which must be satisfied 
at the rock-ore and rock-air interface even for bodies of very simple shapes. 
However, when the ore bodies are at large distances from the surface, the current 


* Manuscript received by the Editor March 5s, 1951. 
' Associate Professor of Physics, University of British Columbia, Vancouver. 
? Geophysicist-Mining Geophysics Corp., Ltd., 11 Jordan Street, Toronto, Ontario. 
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flow near the body will not be appreciably affected by the boundary conditions 
at the air interface, and an approximation of the magnitude of the potential 
differences may be obtained by assuming the body to be a spherical conductor 
embedded in an infinite homogeneous medium, by regarding the current elec- 
trodes as small spheres, and by using the method of electrical images. In such 
a method the electrodes C, and C2 (Figure 1) are replaced by charges of magni- 
tude+p:J/4m where p; is the resistivity of the medium in which the electrodes 









SPHERE 
RESISTIVITY (, 
RADIUS a 


HOMOGENEOUS MATERIAL 
RESISTIVITY (, 

















Fic. 1. Electrical images in a conducting sphere. 


are situated and J is the current flow between C; and C2. The charges at C; and 
C2 induce charges on the conducting sphere whose effect at points external to 
the sphere is equivalent to point charges or “images” at A; and Az, the distance 
of A; and A: from the center of the sphere along lines joining O and Ci, C2, where 
a=radius of sphere, being a2/OC; and a2/OC2 respectively (2). By the use of the 
condition that the current flow across the boundary between sphere and sur- 
rounding medium is a continuous function, the magnitude of images A; and A» 
may be shown to be 
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where pz is the resistivity of the sphere. For a good conducting sphere in a poorly 
conducting medium pp: is negligible compared to p,; and the magnitude of the 
images becomes 


pl a pl a 








— an F . 
4n OC, 4r OC2 


The potential at any point P outside the sphere will be the sum of the potentials 
due to the two objects at C;C: and the images A; and Az. If, as an added simplifi- 
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Fic. 2. Electrical set-up at Drill Hole 57. 


cation, the center of the sphere and all electrodes be in the same plane the dif- 
ference in potential between any two such points P; and P2 will be given by 
AV where 


pil I a/l I aft I I I 
GS-93 GA 
4r b c fitd g fo h Sj k 
where f:=OCj, fo=OC2 and the other letters represent distances as shown in 
Figure 1. 
If CZ is used as a current electrode, then the image Az must be replaced by 


image Az, and the potential, AV, between P,P: will be greater than AV since A’ 
is nearer the potential electrodes than A». When the size of the sphere and its 
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Fic. 3. Vertical N-S geological cross section through Drill Hole 57. 
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position relative to the electrodes is known a value of the ratio AV/AV may be 
obtained. For a conducting sphere of radius 100 feet, whose center is 150 feet 
both horizontally and vertically from C,, and where C; and Cz are 500 feet above 
and 1,000 feet horizontally from C,, and where P,P: are respectively 100 and 200 
feet above Ci, the ratio is 1.04. Since a sulphide ore body of this radius would 
contain approximately one half million tons the effects produced by bodies of 
commercial size would be small but easily detected. The magnitude of the 
calculated results further emphasizes the necessity of making the resistance to 
earth of the surface electrodes exactly the same. 


FIELD TEST OF THE METHOD 


During the summer of 1949 the authors had the opportunity of testing the 
method on the property of the Lake Dufault Mines Limited, where the size and 
shape of a small sulphide ore body were accurately known from extensive diamond 
drilling. The relation of the test hole and the ore body is shown in plan and 
section in Figures 2 and 3. | 











TABLE I 

Distance Pp» AV (North) + AV (West) + 

to collar AV (South) AV (East) 
1,250 feet 1.03 1.02 
1,300 feet 1.02 1.10 
1,350 feet 1.05 5.37 
1,400 feet 1.21 1.05 
1,450 feet 1.20 4.93 
1,500 feet 1.04 1.02 
1,550 feet 1.08 0.79 
1,600 feet 1.09 1.02 
1,650 feet 1.01 1.01 
1,700 feet 1.04 1.04 
1,750 feet 1.14 1.10 
1,800 feet 1.02 1.02 
1,850 feet 0.04 1.00 








Current electrodes were placed on surface 2,000 feet north, south, east, and 
west of the collar of the hole and their resistance to earth made exactly the same. 
One current electrode C, and the two potential electrodes P;P2 were placed in 
the hole. The electrodes were moved along the hole at a constant electrode spac- 
ing, the potential between P,P: and the current through C; being measured at 
suitable intervals. All surface electrodes were in parallel. 

Current electrode C, was next left near the bottom of the hole and P,P: 
moved as a unit along the hole. The potential between P,P: and the current 
through C,; were measured at suitable intervals, using in turn each surface 
electrode as the return current electrode and finally all surface electrodes in 
parallel for the return current. The values AV of the potential between P,P: at 
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each interval when the current flow was to the north or to the west were taken 
as the references and the potential when the main current flow was in the other 
two directions were shown as fractions of these potentials. See Table 1. 
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Fic. 4. Resistivity profiles from Drill Hole 57. 


CALCULATION OF RESISTIVITIES 


Resistivities were calculated on the assumption that the medium was homo- 
geneous from the expression 
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where a and b were the distances C; P; and C,P2 respectively. The results thus 
calculated gave only apparent average resistivities, but the relative magnitude 
of these values afforded valuable information concerning changes in the con- 
ductivity of materials close to the drill hole. The calculated resistivities for. 
constant electrode spacing and for expanding electrode spacing are shown in 
Figure 4. 


DISCUSSION OF RESULTS 


The ratios of the AV/AV were greater than one should expect from the 
simple theory for a spherical ore body and indicated very definitely that the 
ore body was situated to the south and east of the hole. Refer to Table 1. 

The apparent average resistivities remained low as long as the theoretically 
spherical potential shells intersected part of the conductor as shown by the 
dashed circles and the graph of Figure 3. It is of interest to note that the calcu- 
lated resistivities increased when the theoretical shells intersected the body 
where it was interrupted by a known non-conducting sill. 

For a small constant electrode spacing the calculated resistivity values were 
greater than for the expanding electrode spacings, indicating that the low 
resistivity values were not caused by conductors close to the hole. Thus the two 
resistivity curves provided a method of separating the effects of conducting 
bodies such as disseminated sulphide zones at the hole from massive bodies 
situated some distance away. 


CONCLUSIONS 


It would appear that the method provides a powerful means of detecting 
conducting sulphide bodies in the vicinity of drill holes and of determining their 
extent as well as their distance and direction from the hole. 
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A CONDUCTING SPHERE IN A TIME VARYING 
MAGNETIC FIELD* 


JAMES R. WAITt 


ABSTRACT 


The secondary magnetic fields are evaluated for the case of a conducting sphere in a relatively 
poorly conducting medium under the influence of a time varying magnetic field. The sinusoidal and 
step function responses are both considered. The responses so calculated are thought to be useful 
in a geophysical prospecting method which utilizes the transient behavior of induced eddy currents 
in a highly conducting ore zone. 





INTRODUCTION 


The reaction of a spherical conducting ore zone to an impressed magnetic 
field is of great practical importance. The general problem of a sphere under the 
influence of an impressed electromagnetic field has been investigated by Debye 
(1) in a rather general manner. March (2) has treated the more specific problem 
of a conducting sphere in the presence of a magnetic dipole. 

This paper deals with the special case of a conducting sphere embedded in a 
relatively poorly conducting medium. The applied field will be assumed to be 
uniform in the region of the sphere. The secondary fields for sinusoidally varying 
primary fields will be developed from first principles. This of course is then a 
special case of the Debye and March solutions when the external medium has a 
small conductivity. The transient solutions can then be obtained directly by 
suitable integrations of the steady state solutions. The secondary field responses 
to a step function type of magnetic applied field are in turn explicitly calculated. 
M.K.S. units are employed throughout. 


DERIVATION OF THE FIELDS 


A sphere of radius R will be considered of conductivity o1, permeability 4, 
and dielectric constant «, in an infinite medium of constants oe, we, and €:. See 
Figure 1. In the neighborhood of the sphere the applied magnetic field is Ho 
exp (wl). The statement that the external applied field is uniform requires that 
the dimensions of the sphere are much less than a wavelength in the external 
medium. This implies that | yeR| <1 where 72 is the propagation constant of the 
external medium given by: 


¥2 = [iosu2w — eouow?]?/?, 
The sphere conductivity 01, however, can be quite large and there is no restriction 
on the parameter | 7R| where 

11 = [iow — eyw?)?/?, 


* Manuscript received by the Editor April 6, 1951. 
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The center of the sphere is taken as the origin of the spherical polar co- 
ordinate system. The polar axis (i.e. z axis) is taken parallel to the applied field. 

The magnetic field everywhere must be expressible in terms of the magnetic 
vector F with a z component only since the resulting circulating or eddy cur- 
rents are only in the ¢ direction. 


H = — (o + iux)F + 1/ive grad div F. (1) 
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Fic. 1. The coordinate system for the conducting sphere. 


The magnetic vector potential Fy for the primary field is given by: 
Fo = Fy, = — (02+ wee) Ho. (2) 


The expressions for the magnetic vector F; inside the sphere and F, outside the 
sphere are then given by: 


> © baln(yir 
Fi =Fi.= >), betsy) P,,(cos 0) ipyw (3) 
0 is 


a © anK (yor) ; 
F, = Fr, = > ————— Pa(cos 6) iusw — (02 + itces)Ho (4) 
rT 


0 


where 7,(2) =(a2/2)"?In4y12(2), Kn(s) =(22/a)"*Knyie(z) and Inyie and Krys 
are modified Bessel functions as defined by Watson (3). The boundary con- 
ditions are that tangential magnetic field and the normal flux density are 
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continuous at r=R. This requires that a,=b,=0 for n>o. The expressions are 
then of the form: 


Fy, = ipywbo sinh y7/r, (ry = R) (5) 
and 
Fy = — (02+ iwes)Ho + ipewa exp (—yor)/r, (r 2 R). (6) 
The boundary conditions are then: 
1 OF; 1 OF» 
— — - nF: = — — 12°F 
, , OF 
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Mi ar? ae yFy = M4 ar? = "Ye Fs 


These equations can then be solved for ao: 








2ui(sinh a — a cosh a) + pue(sinh a — acosha + a? sinh a) 
a = HR’ ; ; | (7) 


2ui(sinh a — acosh a) — 2ue(sinha — acosha + a? sinh a) 





where 





a= [io mie — eyo? }1/2, 





This result is a special case of the classical problem (4) of a plane wave incident 
on a sphere when the radius of the sphere R and the parameters are such that 
| yeR| Kx. 

The case of interest is when the sphere conductivity is relatively high com- 
pared to the external medium and the permeability contrast is zero (i.e. wi=pe 
=yu=4m X10~"). The magnetic and electrical coefficients then reduce to: 


ao = — 3/2H oR*(M + iN) (8) 










where 





; I I cosh a 
BE > a a See G (9) 
Qe 4 a sinh a 
The total external magnetic fields are then given by performing the operations 


Hw) = — 3/2R°H((M + iN)(— 1/r? + 32?/r5) + He (10) 













and 


H,(w) = 3/2R°Ho(M + iN)32p/r°. (11) 


The in phase factor M and the out of phase component N are plotted in Figure 
2 versus the 0 given by 





6 = (o,yw)1/2R 





where it is assumed that displacement currents are negligible (i.e., 71.>«w). 
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Fic. 2. The behavior of the in phase and out of phase components of the secondary magnetic 
fields for a sinusoidal type applied field. 


The corresponding transient problem can be investigated where again the 
approximations are made, for important frequencies in the spectrum, that 
| yoR| “1 and o;>ew. This will be valid for a step function source magnetic field ° 
when the measuring time of the field response is always such that (>oR? and 
>e,/o4. 

The steady state response of the secondary field (either H, or H, component) 
can be written as follows: 

I I cosh a 

H(iw) = x|= +—-— | Ho(iw) (12) 
a? 3 asinha 

where K is a constant. The transient field is again given by a Fourier Integral 
as indicated by the author in a previous paper (5). 
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K 
h(t) = — 


c+ to | I I cosh a 
271 eS c—ico 


: | nian = ae 
a? 3 asinhe 


When the applied field /o(/) is a step function given by: 
ho(t) = Hiu(t) where u(t) = 1 fort >o 


= o for i<o 


then 
Ho(tw) = Ho/ iw. 


Here Ho(iw) can be interpreted as the frequency spectrum of the step function 
applied magnetic field. On the other hand H(iw) is the spectrum of the secondary 
field components. In the above integral C is taken.as a small positive real con- 
stant. ' 

It is convenient to replace iw by s. 


ak fers I coth Bs!/? o re 
h(t) = : i last - = ral ds where B = (om) ‘ (14) 
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The first two terms in the integrand offer no difficulty and the third term is 
handled by expanding into a series 
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This then gives rise to an infinite series of integrals of the following type (6): 
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The transient response is then given by 






F(t) = a = 42+ 1/3 — an 2x — > 2n|x/n erf’ (n/x) — 2erfc(n/x)| (16) 






where 





a = B/2/B = fl2/(gy)/2R 






In the above expressions erf’ and erfc are the first derivative and the complement 
of the Error functions as defined by Jahnke and Emde (7). The secondary field 
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responses at a point P outside the sphere due to an exciting step function source 
field at the sphere is then given by 


h(t) = 3/2R°H,F(t)[1/r3 — 32?/r®|u(?) (17) 


and 
h(t) = 3/2R*HiF(t)32p/r°u(t). (18) 


The function F(é) is plotted in Figure 3. 
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Fic. 3. The transient behavior of the secondary magnetic field for a step 
function type applied field. 


SIGNIFICANCE OF RESULTS 


The steady state solutions will be discussed first. In Figure 2 the function 
M-+iN is plotted versus a frequency parameter 0. The secondary external fields 
are seen by equations (10) and (11) to be both proportional to the complex factor 
M-+iN. The behavior of the secondary field as the circular frequency is varied 
is the same regardless of the position of the observer. It would then seem feasible 
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that this could be the basis of a geophysical method to ascertain the conductivity 
or radius of the embedded sphere. 

It would be necessary in this case to simply obtain the behavior of either the 
in phase or out of phase component of the secondary field whence the parameter 
o\u.R? could be obtained. A local field could be set up by any suitable ungrounded 
wire loop arrangement. The secondary fields are then obtained by receiving 
loops. The primary field, of course, must also be accounted for in measuring the 
in phase component of the secondary field. 

The parameter o,u.R? might also be found directly by measuring the magnetic 
field response following the application of a step function current (i.e. a suddenly 
applied direct current) to the transmitting or sending loops. The applied magnetic 
field will then be also a step function type. The secondary fields are given by 
(17) and (18). The nature of the decay function F(t) is characteristic only on 
the parameter o,u.R? when the above mentioned approximations are valid. 

The conditions of this problem have seemed rather stringent. The solution 
however should be satisfactory for the case of a large massive sulphide body with 
a conductivity say of ¢,=10? mhos per meter in a surrounding country rock with 
a conductivity of the order of ¢,.=10~* mhos per meter. The body might have a 
radius of say 20 meters. The time constant of decay of the secondary magnetic 
field (where F(t) =0.12) is then equal to 2.4 milliseconds. For this case the primary 
field will have a time constant of the order of a few microseconds, that is, it 
approaches a constant value almost immediately following the application of the 
source current. 

The depth and location of the center of gravity of the body on the other 
hand could best be found by varying the geometry factors of the receiving and 
transmitting loops and observing the variations of the magnitudes of the sec- 
ondary fields as is done in more conventional methods of electrical prospecting. 
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ELASTIC WAVES THROUGH A PACKING OF SPHERES* 


, 


FRITZ GASSMANNT 


ABSTRACT 


Based on a theory of porous solids previously developed by the author, the elasticity of a hexag- 
onal close packing of equal spheres is treated. The packing is anisotropic and because of the weight 
of the spheres, also inhomogeneous. The velocities of propagation of elastic waves have been calcu- 
lated for evacuated interspaces and for interspaces filled with a liquid or gas. In the case of evacuated 
or air-filled interspaces, the wave rays and travel times have been computed. The packing which has 
been treated may be of use as a model for a dry or wet loose material such as gravel or sand. Though 
the model is very simplified, the results obtained show some typical effects such as anisotropy, in- 
homogeneity, and a go° angle of emergence. 


THE ELASTICITY OF POROUS SOLIDS 


In a previous paper by the author (Gassmann, 1951) the elasticity of homo- 
geneous isotropic or anisotropic porous solids has been treated. The following 
investigations are based upon the results of the above mentioned paper. There- 
fore these results may first be summarized, all demonstrations being omitted. 

A sample of a porous solid possessing the volume V and the mass m, may be 
composed of a homogeneous isotropic solid material of volume V and mass m 
and of a homogeneous liquid or gas of volume V and mass m filling the pores. 
Therefore, 


SS 


V=V4+V, m=n4+%, 


(1) 


Z 


n = — = porosity. 
V p y 


Macroscopically the porous solid may be considered as a homogeneous elastic 
solid if the pores are assumed sufficiently small. Though the solid material is 
assumed isotropic, the porous solid can behave either as an isotropic or aniso- 
tropic elastic body according to the shape of the pores. All pores are assumed to 
be in connection. Let us consider an initial state of stress and strain (e.g. such as 
the state of rocks under the weight of overlying material) and variations from 
this state small enough to allow the application of Hooke’s law. Two extreme 
cases will be taken into account. In the first case the variations of stress act 
only upon the solid matter, the hydrostatic pressure of the pore content remaining 
constant, due to circulation allowed through the pores and through the surface 
of the body. The porous solid then behaves like an “‘open system” (e.g., slow varia- 
tions of stress applied to a rock sample in free air). In the second case circulation 
of the pore content is not possible, the porous solid behaving like a “closed 


* Manuscript received by the Editor May 21, 1951. 
t Professor of Geophysics and Director of Institute of Geophysics, Swiss Federal Institute of 
Technology, Zurich. 
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system” (e.g., stress variations in the earth’s interior due to earthquake waves). 
Notations: 


x, y, 3 = rectangular coordinates 


S2, Sy, Sz = rectangular components of displacement (Love, 1944), due to stress 
variations. 


Ae; = components of strain variations. 




















OS: OSy OSz 
Ae = ? Aeo Sy Aes = 
Ox dy 0z 
(2) 
OSy OSz OS; OS; OS; OSy 
Aeg = —+ ’ Aes = + ’ Aes = pe 
0z dy Ox 02 oy Ox 


Ap;:, Ap; = components of stress variations of the open and the closed system 
(pressure with positive sign). 


X, a, k=A\+—R (3) 


2 
3 
Lame’s constants and modulus of compression of the solid material. 


k=modulus of compression of the pore content. Stress-strain relations, 
general anisotropic case: 


open system: 


$ 4= 1, 2, , 6 
Api = — Dade}. 
j=1 Eyres ’ 6 
6 (4) 
Aes = — >. ViiADy;, Cig = Cy, Vii = Vi 
j=l 
closed system: 
6 6 
Api = — Do cishe;, Ae: = — Do viiAp; (5) 
j=1 j=1 
Cig = Ci, Vii = Vii 
Relations between the elastic constants: 
Cig = Cig a ty bib; 
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€1 = 62 = €3 = I, €4 = 6 = & =O 
€; 


a; = Yi + Vi2 + ‘s- = 
3k 


b=«&- Bs (Cir + Cin + Eis) 
3k 
D = 1+ a*(a, + a2 + as) 


* 


1+ — (hae ee. 
3k 


D* 


In particular, for the isotropic case: 








- si ea 2 
A,f,k =A+—aG elastic constants of the open system. 
a 
2 
Aw R=A+— pw elastic constants of the closed system. 
2 
From (7) is obtained 
~ k+ 
h=ii kR=k-x e 
k+Q (8) 
0 R(k — k) 
= n(k — R) 


HEXAGONAL CLOSE PACKING OF SPHERES, HOMOGENEOUSLY STRESSED 


Let two homogeneous and isotropic spheres of equal radius R (with Lamé’s 
constants \, #) be in contact under the pressure P (Figure 1). 


B(3X + 22) 








E=-—~ —— = Young’s modulus, 
A+ uh 
x 
¥ = —~—— = Poisson’s ratio. 
2(\ + 2) 
According to Hertz’s theory (Love, 1944), the radius of the circle of contact is 
r= 4/ ss (9) 





4E 


R’ is assumed to be small compared with R, and 





IIe 


et AGO 
S = ——__>=——— = relative displacement of the spheres. (10) 


oe 16F2R 
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Fic. 1. Two equal spheres in contact. 





Fic. 2. The hexagonal close packing of spheres. 
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The hexagonal close packing of equal spheres (Figure 2) now to be considered 
(following an idea of Hara, 1935) consists of different layers (Figure 3). The 
coordinates of the centers of the spheres belonging to the Nth layer (no pressure 


assumed) are 


x 


[4N’ + 2N” + (—1)¥ +1] 
(rz) 


hae [4N” + (—1)¥ +1] 


RN, N’ and N” being arbitrary real integers. 





As initial stress, we assume zero pressure between the spheres of the same layer, 





Fic. 3. A layer of spheres of the hexagonal close packing. 


i.e., the spheres with centres of the same value of z, and a pressure P>O between 
each pair of spheres in contact, belonging to neighboring layers. The correspond- 
ing slight deformation of the lattice (11) of the centers can be neglected in the 
following considerations. If the number of layers and spheres is sufficiently large, 
the packing at its initial state of stress can be considered as a homogeneous 


anisotropic porous solid 





Tv 
porosity » = 1 — == 0.260 | (12) 
| 3V2 
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with the symmetry of Voigt’s class Nr. 26 of crystals. With respect to elasticity 
all directions perpendicular to z are equivalent (“transversely isotropic” sym- 
metry, Love 1944, page 160) and the elasticity is determined by five constants 





Ci, Co, ..., Cs, the matrix of the open system constants ¢;; in (4) taking the 
form 
C; Ci — 2C, C2 ° ° ° 
Ci = 2C5 Ci Ce Oo Oo ie) 
Ce C2 C3 ° ° ° 
: . rs (13) 
° fe) ° C4 ° ) 
° re) fe) re) C4 ) 
° ° ° ° ° C; } 





The é;; can be taken from (4), if for suitable homogeneous deformations of the 
lattice (11) the variations of stress are calculated by means of (10), variations 
As and AP being replaced by ds and dP. With 


F I t 2p I : 6E2P 
ios + WY oe (14) 
12/2 R ds 12/2 (1 — ¥*)?R? 
the matrix é;; is obtained by putting 
Ci = ic Ce = C4 = 4¢, C3 = 166, Cs =o (15) 


into (13). 
THE PACKING OF SPHERES STRESSED BY SELF-WEIGHT 


The packing of spheres introduced by Figures 2 and 3 and the formulas (11) 
may now be considered under the influence of self-weight, the +2 axis pointing 
vertically downward. The interspaces between the spheres may be filled with a 
liquid of constant density ). Each sphere belonging to the Nth layer is in contact 
with three spheres of the (V+1)th layer under the pressure Py and with three 
spheres of the (V—1)th layér under the pressure Py_;. The condition for equi- 
librium of such a sphere contains the vertical components of the six pressures, 
and the apparent weight 47/3R*(p—p)g of the sphere when g=980. 665 cm 
sec’. It is 


bie = 47 me . 
etiam (16) 


If no load is assumed above the first layer, then Po=o and the recurrence formula 
(16) leads to 


OO a os ae 
Ba = 8 Oe Oe 


3V (17) 
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or taking the value of NV from (11): 
w 4 _~ ~ 
Py = — RG — 3)ge. (18) 
a 


Substituting this expression for P in (14), we obtain 


tg / 2 EG — Dg 
¢=—— 
12/2 (1 — ¥?)? 
If the pores are filled with a liquid, the packing of spheres may be called a “‘wet 
system.” If the pores are evacuated [7 =o], the packing may be called a “dry 
system.” If the pores are filled with air, its density and its modulus of compression 
are negligible. Therefore the system possesses the elasticity of a dry system. 





Ve. (19) 





p = (1 — n)p = density of the dry system. 
p = (1 — n)p + np = density of the wet system. 


(20) 


(n = 0.26, see (12)). 


The stress below the plane z=z9>0 (zo large compared with R) is not changed 
if the packing of spheres between the planes z=o and z=2o is replaced by an 
arbitrary layer between the planes z=2, and z=2o, with density p*, depending 
only on 2, if 


is p*(z)dz = pZo. (21) 


od | 


The packing of spheres may be considered as a simple model for a gravel or 
sand, dry (=dry system) or filled with ground water (= wet system). 


PROPAGATION OF ELASTIC WAVES THROUGH THE PACKING OF SPHERES 


The elasticity of the dry system is given by the matrix of the é,;; which is 
calculated by means of (13), (15) and (19) with =o. In order to calculate the 
elastic constants of the closed wet system, first the matrix of the é;; is calculated 
in the same way, but with p 0 in é of (19). Then the matrix of the c;;, found 
by (6), is the matrix (13) with 

a* a* 
Ci re wes Or = ee 
a* 
C3; = sail hell Cy = 46, Cs = 0» 
(22) 


C 
(=m= 1-25, b= 1-8=> b4=b = bg = ). 


iy) 


=) 


‘he partial differential equations of wave motion, the wave length assumed 








CO 
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large compared with R, are obtained by substituting the stress-strain relations 
into the general equations of motion of an elastic solid (Love, 1944): 















































d’s, 0°55 ’s, 0°S 2 Osy 
p = C1 + Cs + C4 + (Ci — Cs) 
ot? Ox? dy? 02? dxdy 
0°s, 
+ (C2 + Cs) ’ 
0x02 
0°Sy _ Oy , Oy d’sy “y Om 
poe Ge re ee ey ir t Ga be a> 
ar? ax? ay? dz? se (23) 
_, OS 
+ C, — C,)-——_ 
Oxdy 
0s, 0°s, 0’s, 0°s, 0°S 2 
pa kg os ae + (C2 + C4) 
ar? ax? y? 02? O0xd0z 
0S, 
+Cs+Cy——>»> 
0 ydz 


where /=time. For the wet system, the C;’s are to be taken from (22), for the 
dry system from (15), putting p=o in (19) and (20). Now let W(z, y, 2, ?) 
=const.=equation of a wave front. 

ow ow ow ow 


— =q@, —=%, —=%, —= qQ. 2 
ot % Ox “ Oy 4 4 % (24) 


The characteristic relation (Levi-Civita, 1932) of the system (23) is 


| Ais — pqo? Aj2 Ais | 
| Ae Ao — pqo” A233 = 0, (25) 
| A43 A 33 A33 — pqo* 

Ay = Cigi’ + Csq2” + C'4q3", Ay» = (Cy = Cs) 4192; 

Ao = Csqi? + Cige? + Cags’, Ao3 = (C2 + Ca) G29, 

A33 = Cuqi’ f: C4q2” + C3q3", Ay3 = (C2 + C4) 9193. 


Assuming a point source (=focus) at «=y=z=o and considering the wave 
propagation in the (x,z) plane, we have by reason of symmetry, g.=0 


= velocity of propagation perpendicular to the (26) 
2 





wave front. 


qi P 
tan a = —» a = angle between the normal of the wave front and the z axis. 
q3 
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The three roots of (25) corresponding to the three different kinds of waves with 


. velocities v=, v2, v3 lead to 








et I 

= (An + Ass + VAii? + Ags? + 4A13? — 2A11A38 

V9” 2p(qi? + gs”) ; (27) 
ee Ao» 


. ee 
p(q1” + 9s”) 

As an example, the velocity v,(v:>v2>3) has been calculated for a packing of 

granitic spheres, dry [m= (z,a)] and filled with water [n=%*(z,a)] for 

a=o° (vertical direction of propagation) and a=go° (horizontal direction), 


using the following values: 

































































“ " dyne 
p = 2.65 —» E= 0.5:10!? —_; vy = 0.25, 
cm? cm? 
o dry, (28) 
gr am dyne 
p=|)1 —wet, k = 2.06:10! 
cm’ cm? 
The results are shown in Figure 4. 
«— velocity w in mec 
3 
2 0 ° 
fi Va 10 E 
ae 
| /s i a vertical 
elocit 
| | £ ——-— horizontal eee 
| ror 
3 
| r -—— velocity wv in mec 
° 
| 8 g 8 g ° 
F — wn | ai 
ios | I al wa +s 
a a 
& + | ‘1400 / £ 
| : 
| | | £ 
:/ a/ : 
| | i ; | 
| | | 
| | [ 
| | | 
| | 
oe y m $ rf ia 
Cia od Gi —— — 
wet. dry dry 








Fic. 4. Velocities of wave propagation in a packing of spheres. 
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The wave rays (=seismic rays) are the bi-characteristics of (23). For the 
rays in the (x,z) plane (qg2=0), the differential equations are: 









































dx oH dz oH 
— = = — = ——"9 
di am «= dt—s qs 
(29) 
— = —_— 5] — = — — se 
dt Ox dt Oz dt ot 
0 10 20 ]40 =e 
standard horizontal distance F 
4 
10 
8 
12 
20 
wo 
ie i 
me | . 
o 30 
® 
Qa 
8 
x | 
3 
2 lao 
o 
a a | 
@ 
40 
45 ; l 
= ~~ ' 2 & 





Fic. 5. Standard ray of elastic waves in a packing of spheres. 


H, corresponding e.g., to 1, is obtained by taking the square root of the cor- 
responding equation (27): 


go = — H(z, 1, 9s). (30) 


The rays corresponding to 2, will now be calculated for the dry system: 
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a = ar WeV 5:2 + 4932) + Vout + 3289: gs? + 144934 , 


‘ j/ “BEV 2g 
OV wp(x — 9?) 


We introduce @ as parameter 





(31) 








Poe 
F(a) = —V a5 + 15 cos 2a + (15 + 9 cos 2a)" + 256 sin* 2a , 


d 
~ = F(a), (32) 


1=L- </z-F(a), 
F(a) 


sin @ 





H=L-Wz- 4 


‘and integrate (29). 

If =o is assumed as focal time, the results are: 
sint y 

0 Fo FOr” 

(auxiliary function, graphically computed) 


T, = T(r) = 141,950 (34) 


T(a) = 5 (33) 


(integration by aid of elliptic integrals) 
Standard ray (Figure 5): 


100 T(a) — [=] é \ 
é= =| (a Fla) cot a 


m al =) 
c= TL F@) (35 








Ioo 


rf’ T(a) (=standard travel time) 
1 


i= 
angle of emergence = go° 
Arbitrary ray: 
a = horizontal focal distance of the ray. 
a a 


= a” 2 = es 
Io0o I0o0 


__ gbié qs/6 
t= see as =e = 0-caen 4 . (36) 
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0.521 
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O 0.372 
® 
| ¢| Sg 
° ee as aa” 
w) S7B ea 
515,,— fells) | 
a t 
oe + Q209 
+ - ' y; 
ee 0.2 A393 ' on 5 
H ' Oo 
Te iS 8 
= aN 1o117 2 > 
0.1 By 200 
| 40.031 
.e] —— 
10 50 100 200 300 
—-+ distance a in m 
Fic. 6. Travel time and apparent velocity of elastic waves in a packing of spheres. 
The ray is symmetric with respect to x=a/2. 
Travel time function: 
2.7408 ; 
a=", *£=4 s=0, %t=r= 7 a5/é = (Figure 6). = (37) 
Deepest point of the ray: 
T a 64 
a=—) x= —) Z= % = —a = 0.4514. (38) 
2 2 T, 
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LEAST SQUARES RESIDUAL ANOMALY DETERMINATION* 


W. B. AGOCS} 


ABSTRACT 


The residual anomaly is defined as the deviation from the mean anomaly surface, or the regional 
surface. The regional surface which best fits the observed anomaly data may be determined by least 
squares. For the case of the simple plane, the equation of the regional would be Z=Ax+By+C, and 
the residual anomaly would be R=G—Z=G—(Ax+By+C), where G is the observed value at the 
station whose coordinates are x, y..The constants of the equation of the regional may be determined 
by using the normalizing equations ) .RAR/dA =o; >. RAR/dB=o; > ,RAR/IC=o. It is shown that for 
a symmetrica! distribution of observational points about a center, the regional value at the center 
would be equal to>.?G,/n and this value subtracted from the central point is the residual. Given the 
coordinates and the observed values, the regional may be determined, and the residual for the area 
calculated rapidly. An example is given of the procedure. 


INTRODUCTION 


The data mapping potential fields is an integrated effect of all sources causing 
an observed anomaly. In applied geophysics the problem is to eliminate or reduce 
to a minimum the effects of deep-seated, non-commercial sources with as little 
distortion of the resultant anomaly as possible. The major field effect which is 
eliminated is the massive regional, and the minor deviations of the observed 
data from this regional is the residual anomaly. Numerous mention has been 
made of the value of eliminating the regional and using the residual in inter- 
pretation procedure, but there is a paucity of general information on the subject. 

The residual anomaly determination is to be studied and discussed in the 
following. The residual anomaly may be defined as the deviation from the mean 
or regional field surface. This may be defined mathematically as 


R=G-Z (1) 


where R is the residual anomaly, G is the observed anomaly, and Z is the regional. 
The methods of determining the regional for an area from which the residual 
-is determined are numerous, and the value of the residual map is dependent on 
the knowledge of the method used to determine the regional and the assumptions 
underlying the choice of method. In the following the various methods in u8e are 
discussed, although combinations of the methods are frequently used. 


REGIONAL AND RESIDUAL DETERMINATION 
Arithmetic Mean 
One of the simplest and most rapid means of obtaining the regional is by the 
* Read by title at the St. Louis Meeting of the Society, April 25, 1951. Manuscript received by 


the Editor June 6, 1951. 
T University of Tulsa, Tulsa, Oklahoma. 
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Fic. 1. The combined field of a buried sphere and a south gradient of 
40 gravity units per kilometer. 


method of arithmetic means. This method is illustrative of the interdependence 
between the chosen aerial distribution of control and the residual deviation. If 
the area of control points chosen is small, the regional value approaches the 
observed value, and the residual anomaly approaches zero, and if the area of 
control is chosen extremely large, the regional approaches zero, and the residual 
value approaches the observed anomaly value. Hence, in this method it is neces- 
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Fic. 2. The regional determined by the arithmetic mean, using the data shown in Figure 1. 


sary to choose an area intermediate between the two extremes which will limit 
the effect of the regional elements but will show the residual with the least 
distortion. 

The value of the arithmetic mean for the regional is placed in the center of 
the area of control, and these mean values or regional values are contoured. The 
contoured results of these mean values is the chosen regional. 
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Fic. 3. The residual anomaly using the arithmetic mean regional. 


The data shown in Figure 1 were used to illustrate the method. Figure 1 
shows the field due to a spherical mass of density contrast unity on which is super- 
imposed a south regional of 40 units per kilometer. The depth to the center of 
mass of the disturbing mass is 1,000 meters, and its radius is 500 meters. The 
maximum deviation due to the anomalous mass is 34.8 gravity units. 

An area of 1.5 kilometers square was chosen, with the exception of the edges 
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of control where the areas are 1.5 by 1 kilometer to illustrate the use of this 
method. Figure 2 shows the area of control with the arithmetic mean for each 
area and the resulting contoured regional. 

Figure 3 illustrates the residual anomaly obtained by subtracting the regional 
from the observed values. 

The arithmetic mean is affected by the anomaly which is being sought. This 
may be observed by referring to Figure 2, which shows a broad, southward 
plunging regional nose. In the interpretation of the regional using this method, 
consideration must be given to the possible distortion introduced into the regional 
by the shallow or local anomalous effects. 

The residual anomaly obtained by the method shows a maximum anomalous 
deviation of 30 units. Furthermore, the positive anomaly is “‘haloed”’ by negative 
anomalies which are caused by the observed low reentrants which were not 
affected by the regional’s smoothing effect. The negative anomalies are pseudo 
effects which must not be considered as structural in the interpretation. 


Contour Smoothing 


One of the oldest methods of determining the residual anomaly is contour 
smoothing. The name of the method describes the procedure for determining 
the regional. The observed anomaly contours are arbitrarily ‘smoothed”’ by the 
interpreter. These smoothed regional contours are then subtracted from the 
observed data. 

There are two possible methods of determining the residual anomaly from 
the observed data. One method is to subtract the regional contours from the 
contours of the observed values. The second method of procedure is to subtract 
the regional contour values from the observed station values. The former method 
is that commonly used in geological work to determine structure prior to regional 
tilting. 

The method is similar to that described previously except that the smoothing 
is done arbitrarily rather than through the use of the observed data, and it 
depends mainly upon the judgment of the interpreter of the data. 


Method of Profiles 


The method of profiles is similar to contour smoothing, excepting the pro- 
cedure for doing the smoothing. Profiles are drawn on the anomaly map in the 
direction of maximum gradients. The distance of separation of the control 
profiles is usually of the order of the width of the local anomalies. Along the 
length of the profiles a constant gradient is assumed. This constant gradient is 
taken as the regional gradient. Along zones of gradient changes in the observed 
data, the regional gradient is adjusted or changed along the profile. 

The regional gradient as indicated along the profiles is contoured smoothly, 
and the resultant is the regional field. 

The residual] anomaly is obtained by subtracting the regional contour value 
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from the observed station values. With closely spaced profiles and averaged 
gradients approximately equal to the observed gradients the residual anomalies 
approach zero due to the elimination of all but sharp, local variations. 


Average Profile Method 


The area covered by the survey is broken into a grid by intersecting profiles, 
regularly spaced. The spacing of the profiles is chosen usually equal to the 
station spacing. The values at the intersections of the profiles are used for the 
calculation of the regional by interpolation of the observed values. 

In this method, a running average of the values is carried, and the average 
value, which is the regional, is placed at the midpoint of the segment over which 
the averages were obtained. At each point of the grid, two average values are 
determined, one for the x direction and the other for the y direction. The average 
of these two average values gives the regional for the point of the grid. 

The average profile method or grid method is a least squares solution for 
the regional. Assuming that the regional Z is a function of x or y, then the 
probability of obtaining the observed value at the point x, y is a maximum when 
D1R,2 is a minimum. It follows that the necessary conditions for a minimum 
for an area of plane surface, where 


R=G-— (Ax+ By+C) 
are that 
>> RAR/IA = 0; >> RAR/dB = 0; >> RAR/IC = o (2) 


since R is a function of the unknowns A, B, C. Equations (2) are the normalizing 
equations. For a station located at the intersection of a pair of profiles, and for 
symmetrically distributed stations on the profiles, with the coordinates of the 
point of intersection «=o and y=o, the constants A and B are equal to zero, 


and the regional 
Z =C=( YG.00)/ (3) 
1 


The residual anomaly at the point of intersection would be 


R=G-2=6—(YGu00)/m (4) 


The average values (ZGn, 2, »)/m are contoured, and the resultant is the 
regional surface. The residual may be determined by subtracting these regional 
values at the grid points as indicated by Equation (4), or by subtracting the 
regional at the points of observation. 

The determination of the residual at the grid points gives a false indication 
of the density of control. Furthermore, the control of the regional will be de- 
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pendent on the control exercised by the original contours on interpolated values 
at grid points. The determination of the residual at station locations from ob- 
served values may be contoured on the basis of the observed data, and false 
anomalies will not be introduced. 


Least Squares Determination 


The best fitting regional to the observed data may be determined by least 
squares methods. The method is simplest to use when the observed gradient 
indicates a uniform plane, although it is possible to fit more complex functions 
to the observed field. 

The procedure followed is similar to that outlined above in the section on 
average profile method. Assuming that a plane would best fit the data, the 
equation for the plane is 


Z=Ax+ By+C (5) 
where Z would be the computed value of the regional for the coordinates x, y. 
A, B, and C are constants which are to be determined. 
The residual or the deviation would be 
R=G—-—Z=G-— (Ax+ By+ C). (6) 
The necessary conditions for the sum of the squares of the deviations or residuals 
to be a minimum are set forth in equation (2). These conditions may be shown 
to be the following in terms of the coordinates of the points of observation and 
the observational values 


> ROR/A = > Ax? + DS Bry + D> Cx — DG = 0, 
do RAR/dB = Yi Any + Di By’ + LiCy — DGy =o, (7) 
>, RAR/AC = Do Ax+ DS By+ DC —- DG=o. 
Equations (7) gives us three independent equations in terms of the coordinates 
and the observed values to solve for the three unknowns A, B, and C. 

The method may be illustrated by the use of the data shown in Figure 1. 
Only the observed data are used, and it is found that the following are the terms 
- for the normalizing equations (7) 

x y G Gx Gy x? y? “y 
1,768 1,326 35,962 287,670 277,632 19,448 11,050 10,608 
For the above determinations 221 stations were used, and the o, o coordinate is 
in the southwest corner of the area of control. A distance of 500 meters between 


stations was taken as a unit of distance. Substituting the above tabular values 
into the normalizing equations, the following are obtained 


19,448A + 10,608B + 1,768C = 287,760, 
10,608A + 11,050B + 1,326C = 277,632, , (8) 


1,768A + 1,326B+ 221C = 35,962. 
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Fic. 4. Residual anomaly obtained by least squares fitting of data. 
The solution of equations (8) yields the following for the regional 


20.03 + 42.49 = Z. (9) 


The resultant residual anomaly is shown in Figure 4. It may be seen that the 
result of the inclusion of the anomaly is to raise the plane of the regional by two 
units above that of the true regional whose equation is 


Z = 20y + 40 
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Fic. 5. Observed field data, 


where Z is the regional and a unit distance is equal to 500 meters, The residual 
anomaly more nearly approximates the true anomaly, but this would be expected 
since the original anomaly was added to a plane regional. 

Another set of data shown in Figure 5 has been used to determine the residual 
present. The equation of the regional here is found to be 


Z = 0.36% — 5.04y + 93-5. 
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Fic. 6. Residual anomaly determined by least squares fitting of the observed data of Figure 5. 


The resulting residual using the above equation for the regional is shown on 
Figure 6. 


Circle or Polygon Summation 


In circle or polygon summation, a template is made of arbitrary radius, and 
four, six, eight or ten points are placed at regular intervals around the circle. 
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The average of the sum of the observed values at the points on the circle is the 
regional value at the center of the circle. 

The method is a least squares fit to the control of circular symmetry, since 
for circular symmetry, with the center of the circle as the origin, the constants 
for the equations of a plane are 


and 


cn (3 o.)/naz. 


1 


Hence the regional value at the center of the circle of summation is the average 
value of the observed values on the circle of control. 

The radius of control of the circle is chosen so that it is equal to or greater 
than the width of the average anomalies. Using this radius minimizes the forma- 
tion of pseudo negative anomalies when averaging over the flanks of positive 
deviations, and also reduces the formation of rim or “haloing” highs about a 
true maximum point. 


CONCLUSIONS 


This paper has set forth the various methods commonly used for determining 
the regional anomaly, and therefrom the residual. The most rapid method for 
determination of the regional and the residual is by arithmetic means. However, 
the average profile, least squares and circle summation or polygon methods 
are variations of least squares fitting of data. The fitting of a regional to the 
observed data by least squares reduces the determination of pseudo anomalies 
which may result from the interpolation of observed data when using profile 
averages or circle summation. The study of the observed data in conjunction 
with the regional that has been determined and the residual determination will 
aid in deciding the validity of structural interpretations. 
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THE WORLD-WIDE GRAVITY PROGRAM OF THE MAPPING 
AND CHARTING RESEARCH LABORATORY OF OHIO 
STATE UNIVERSITY* 


W. HEISKANENT 


ABSTRACT 


This paper deals with the geodetic applications of the gravity method. Gravity anomalies permit 
the absolute undulations of the geoid de 3 the absolute deflections of the vertical to be determined. 
In turn, these, together with astronomical observations, form the basis for a World Geodetic System 
for the control of small scale maps beginning with the scale 1: 100,000. They serve also, together with 
existing triangulations, for the correction of the dimensions of the reference ellipsoid. Additional 
gravity surveys are needed in order to obtain these and other minor objectives. Cooperation among 
scientists and agencies of different countries concerned with these problems is a prerequisite for the 
success of this world-wide program. 


Gravity measurements can be used for geophysical, geological, prospecting, 
and geodetic purposes. Thousands of publications have been written about the 
geophysical and geological applications of gravity methods. The whole science 
of isostasy is essentially based on gravimetric studies. The famous strips of nega- 
tive gravity anomalies in the East Indies, ir the West Indies, to the east of Japan, 
and west of South America have been discovered gravimetrically. The thickness 
of the first layers of the Earth’s interior can be computed gravimetrically as 
well. The gravity method has been an important one for the oil industry and 
also for the mining industry in general. 

When I was invited to Columbus last August, Professor Harding, Director 
of the Mapping and Charting Research Laboratory, asked me to study to what 
practical geodetic purposes gravity anomalies can be used. These investigations 
join closely to the studies which I have done with my co-workers in Finland for 
twenty years. The substance of all these studies is this: 

The gravimetrist has only one tool, gravity anomalies, but this tool is very 
sharp. It can solve many problems, as follows. The gravity anomaly Ag is ob- 
tained in reducing the observed gravity value g to sea level, either by the aid 
of free air, Bouguer or isostatic reductions and subtracting from this reduced 
value go the normal gravity value, y, obtained from the gravity formula. 

The geophysical basis of geodetic applications of the gravimetric method is 
the fact that all the quantities: the undulations of the geoid N, deflections of the 
vertical components £ and 7, and gravity anomalies Ag have the same cause, 
brought about by the same thing, the disturbing masses of the earth’s interior. 
Ag can be observed, WN, &, and 7 can be computed. 


* Presented at the St. Louis meeting of the Society, April 25, 1951. Manuscript received by the 
Editor July 5, 1951. 

t Director of the Finnish Geodetic Institute, Consultant of the Mapping and Charting Research 
Laboratory of Ohio State University, Scientific Director of the Institute of Geodesy, Photogrammetry, 


and Cartography of Ohio State University. 
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The geodesist has the following difficulty. The geodetic measurements are 
reduced to an equipotential level, the geoid, and his astronomical observations 
are referred to the geoid as well; but he is compelled to carry out his geodetic 
computations along another more simple surface, the ellipsoid of revolution in 
order to get the coordinates of the control points. He must therefore have some 
method for converting the measurements from the geoid to the ellipsoid, i.e., for 
computing the quantities N, & and yn. He has not had such a method so far. 

As the classic arc measurements cover only a tiny part of the continents and 
they fail, at least up to now, totally in the oceans, the arc measuring method can 
give only the general size of the earth but not its exact shape. Only in some few 
parts of the world have the quantities NV, &, and » been computed by means of 
this classic method, and even there the obtained values have not been correct. 

The only method, which can give N, &, and 7 is the gravimetric method. In com- 
puting the undulations N of the geoid we need only to use the famous Stokes’ 
formula, derived more than a century ago in 1849. The Stokes’-Vening Meinesz’ 
formula (1928) gives the deflection of the vertical components & and 7, (we sign 
them now with the subindex g). These gravimetric quantities NV, &, and 1, 
are absolute ones. They indicate the real undulations of the geoid and real angle 
(deflection of the vertical) between the normal of the ellipsoid used and the 
plumb-line. 

What is needed is a fairly detailed gravity anomaly map of the vicinity of 
the points where it is desired to compute the quantities NV, £,, and n,. In addition 
we need the general gravity anomaly maps over the world, which give the gravity 
anomalies in broad lines. As exploration geophysicists know, there are large areas 
in the United States, Canada, Europe, India, and small regions in many countries 
of the Middle East and Central and South America where detailed gravity 
anomaly maps exist. 

In spite of the existing millions of gravity measurements in the world, how- 
ever, we need for our purposes more /Jocal gravity surveys in many places as well 
as regional gravity measurements in some continents and oceans. If such new 
gravity measurements were to be carried out at nearly equal distance from one 
another, about 3,000 to 4,000 stations could give the general gravity field fairly 
well. 

As to the local gravity field, about 100 gravity stations are needed inside the 
circle of 200 km radius in the vicinity of each deflection of vertical station. If more 
gravity data than this are available we, of course, use them all. 

Under these conditions we can—omitting some minor results—obtain the 
following four objectives: . 

1. That of computing the undulations N of the geoid with an accuracy of 
2 to 5 and the deflection of the vertical components & and 7, with an ac- 
curacy of 0.7” to 1.0”. This can be done, e.g., at the initial points of the 
existing many different geodetic systems (North American, South Ameri- 
can, European, Indian systems, etc.). In this way we can convert all these 
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systems to the same World Geodetic System. We have only to tilt and 
either uplift or lower each of these systems so much as the quantities &, ng 
and WV indicate. 

2. That of converting the astronomical coordinates and azimuth from the 
geoid to the reference ellipsoid when we know the quantities & and ny at 
the astronomical points. That will be accomplished by means of the equa- 
tions: 

ou - §, 
A=) — a, sec ® 
A = A’ — 7, tan ® 


where ®’, \’, A’ are the astronomical latitude, longitude, and azimuth 
referred to the geoid, and where ®, A, A are the corresponding geodetic 
quantities referred to the reference ellipsoid. Thereby we would establish 
the geodetic coordinates of all astronomic points in the World Geodetic 
System and, incidentally, the distance between any such points along the 
ellipsoid used, if necessary. 

These distances can be easily converted from the ellipsoid to the geoidal 
distances. The point error will be about 1.0” to 1.5” or 30 to 45 m, and 
the error of the distance +/2 times larger, or about 45 to 65 m. 

3. That of controlling maps of scale 1:100,000 and maps of even smaller 
scale without any triangulation. This accuracy will not decrease with the 
distance, as is the case in the triangulation. The point accuracy noted in 
the second objective of 30 to 45 m corresponds to an accuracy of 0.3 to 
0.45 mm on 1:100,000 scale maps and to an even smaller amount on maps 
of smaller scale. 

4. That of correcting the dimensions of the used reference ellipsoid (as Lam- 
bert in 1947 pointed out) when comparing the gravimetrically computed 
quantities £ and n, at the astronomical points of the existing arc measur- 
ings with the corresponding quantities and 7 obtained from the arc 
measurings. 

If the differences (£,—£) and (y,—7) are small along the whole arc, the 
ellipsoid used is good, but if (£,—£) and (n,—7) increase or decrease system- 
atically along the arc, we must correct the ellipsoid so that the systemati- 
cal part of (§,—£) and (n,—7n) disappears. 

In collecting the existing gravity material and in analyzing them, such famous 
scientists as Vening Meinesz, Lejay, de Graaff Hunter, Hirvonen, Hess, Ewing, and 
Woollard have promised their wholehearted cooperation. Ewing’s student, Dr. 
Worzel, will continue the world-wide gravity measurements at sea, and Woollard 
and his students will continue the global gravity measuring trips to the different 
continents and to the ocean islands. The ‘‘grand old man” of American geodesy, 
W.D. Lambert, is now a consultant of Ohio State University Laboratory. I my- 
self will conduct the analyzing work. 
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Vening Meinesz, president of the International Union of Geodesy and Geo- 
physics, has promised all possible help also from this Union. 

We wish to get help for the “Columbus team.” When we shortly will send 
members of the Society of Exploration Geophysicists a circular asking for gravity 
data, or still better, for gravity anomaly maps, I hope that members will be 
willing to let us have such gravity data. The laboratory doesn’t need the excep- 
tional values but only the mean values of squares, say of 50 to 100 square miles. 
All data will, of course, be used confidentially. 

If gravity data contributed are not referred to sea level and are not tied to 
the general gravity systems, the “Madison team” of Professor Woollard will be 
ready to carry out necessary connection measurements. The world-wide gravity 
measuring trips of Woollard and of his students have basic significance for all 
geodetic application of the gravity measurements. Only by means of such world 
measuring trips can we reduce gravity values of different countries to the same 
system. 

If it be presumed that necessary financial support will be available in addition 
to cooperation among the agencies and scientists concerned, the “Columbus 
team” will be able to accomplish the aims outlined here in the next three years. 
The geodetic problems involved are perhaps the most important in this century. 

The International Union of Geodesy and Geophysics at its Ninth Congress 
(August 21 to September 1, 1951 at Brussels) adopted the following resolution 
concerning this world-wide program: 


The International Union of Geodesy and Geophysics, considering 


(a) that a more and more precise determination of the size and shape of the Earth and the 
fixing of positions on the same reference system constitutes one of the essential aims of the 
Association, 

(b) and that, in spite of the gaps which exist in the world gravity system, a better knowledge of 
the size and shape of the Earth could be attained by the use of existing gravity data in associ- 
ation with existing triangulation, 

(c) and that some indication might also be obtained concerning the differences between various 
systems of reference, 


records 


with great satisfaction the establishment of the large program of research by the Mapping and 
Charting Research Laboratory of the Ohio State University of Columbus (Ohio) and strongly ap- 
proves its objectives, 


and considers 


it necessary that, in order to make this program conform with the aims of the Association close 
collaboration should be established with appropriate persons and institutions in order to study 
(a) gravity values, (b) such methods of reduction, and (c) such methods of calculation as will 
insure the greatest possible accuracy of the results and allow new gravity data to be used for 
future revision and extension. 
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O. F. RITZMANN* 
GEOCHEMICAL PROSPECTING 
U. S. No. 2,551,449. J. R. Menke. Iss. 5/1/51. App. 6/15/49. Assign. Nuclear Development As- 
sociates, Inc. 
Method for Locating Deposits. A geochemical-prospecting system in which samples are analyzed 
for an indicative element and the relative abundance of its isotopes measured and plotted. 


MAGNETIC PROSPECTING 
U.S. No. 2,549,845. J. C. Mouzon, E. S. Lamar and J. Bardeen. Iss. 4/24/51. App. 6/26/45. 


Ship’s Magnetic Submarine Detector. A submarine-detecting system for a ship having oppositely- 
connected detecting coils located fore and aft normal to the longitudinal axis of the ship, and com- 
pensating coils perpendicular to the detecting coils and to each other in series with the detecting 


coils. 


U.S. No. 2,549,857. E. O. Schonstedt. Iss. 4/24/51. App. 6/27/46. 


Cable-Sus pended Aerodynamic Body. A bird for housing an airborne magnetometer in which the 
magnetometer head may be rotated go° about a longitudinal axis and the center of suspension main- 
tained above the center of gravity to maintain stability in flight. 


U. S. No. 2,550,719. R. Raspet. Iss. 5/1/51. App. 10/14/47. 


Temperature Compensation Means for Magnetic Field Balances. A magnetic field balance whose 
moving system has adjustable masses on threaded arms extending perpendicular to the magnet axis, 
the rods being made of dissimilar metals to shift the mass in a direction parallel to the magnet axis 


with temperature changes. 


U.S. No. 2,551,596. G. O. Haglund. Iss. 5/8/51. App. 5/21/46. 


Aerodynamic Body for Carrying Detection A pparatus. A bird for housing an airborne magnetome- 
ter to which the towing cable is attached through gimbals located at the cross-sectional center of 


form and above the center of gravity. 
U. S. No. 2,555,209. V. V. Vacquier and G. Muffly. Iss. 5/29/51. App. 11/1/43. Assign. Gulf Re- 
search & Development Co. 


Method and Apparatus for Measuring the Values of Magnetic Fields. A gimbal-mounted aircraft 
magnetometer stabilized in two planes with the stabilizations controlled by one or more magnetically- 


sensitive devices. 


U.S. No. 2,556,199. F. W. Lee. Iss. 6/12/51. App. 12/31/47. 


Dynamically Balanced Weight System. A magnetic field balance which is mounted on gimbals 
maintained horizontal by a pendulum and having compensating magnets in the pendulum and also 
in horizontal arms. 


SEISMIC PROSPECTING 


U.S. No. 2,547,703. A. J. Hermont and J. C. Toups. Iss. 4/3/51. App. 7/10/48. Assign. Shell De- 
velopment Co. 


Seismic Exploration System. A seismograph avec system having a multi-stage signal transmission 
channel with nonlinear resistances in a bridge circuit to control the input of each stage, and a control- 
signal channel with rectifiers to supply bias current to the bridges. 


* Gulf Oil Corporation, Patent Department. 
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U.S. No. 2,548,990. K. W. McLoad. Iss. 4/17/51. App. 10/5/48. Assign. Socony-Vacuum Oil Co., 

Inc. 

Invertible Geophone. A geophone in which the parts whose relative motion is detected are inde- 
pendently suspended from the case and arranged so that gravity causes one of them to rest against 
a stop when the geophone is in one position and the other one to rest against a stop when the geo- 
phone is in an inverted position. 


U.S. No. 2,551,417. C. H. Carlisle. Iss. 5/1/51. App. 6/27/46. Assign. Standard Oil Development Co. 


Apparatus for Seismic Exploration. A cable and detector assembly for under-water seismic opera- 
tions having a relatively rigid cable onto which two similar halves of a geophone case may be clamped 
and with the geophone system rotatively mounted in the case with its center of gravity off the axis 
so that the geophone system always remains right side up. 


U.S. No. 2,551,857. A. N. Stanton. Iss. 5/8/51. App. 9/6/49. 

Vertical-Vibration Pickup. An electromagnetic geophone having a horizontal cylindrical case 
with the cable passing through a hole on its axis and an annular magnet with suspended moving coil 
which is made to move by an annular mass on a tube with an annular plate whose lower edge engages 
the coil so that the device will operate in any angular position about a horizontal axis. 


U. S. No. 2,554,905. J. E. Hawkins and J. R. Cornett. Iss. 5/29/51. App. 6/1/46. Assign. Seismo- 
graph Service Corp. 


Seismic Signal Amplifier. A seismograph avc system having two cascaded attenuator networks 
each with variable impedance rectifiers whose bias is obtained from the amplifier output. 


U.S. No. 2,555,806. G. D. Mitchell, Jr. Iss. 6/5/51. App. 3/18/47. Assign. Texaco Development 

Corp. 

Seismic Prospecting Method, Including Generation of a Cylindrical Wave Front. A method of 
seismic prospecting in which a number of shots along a line parallel to the line of geophones are fired 
simultaneously so as to generate a cylindrical wave front which arrives at all of the geophones at the 
‘same time. : 


U.S. No. 2,556,299. C. B. Scott. Iss. 6/12/51. App. 9/12/47. Assign. Standard Oil Development Co. 


Method of Generating Directional Seismic Waves. A method of seismic shooting in which a vertical 
distribution of explosive is placed in a shot hole and detonated at the top, the explosive having a 
detonation rate equal to the seismic velocity in the earth adjacent the shot hole as determined:by a 
trial shot. 


U.S. No. 2,557,714. P. S. Williams. Iss. 6/19/51. App. 9/13/47. Assign. Standard Oil Development 
Co. 


Shot Pulse Recording in Reflection Seismography. A method for improving the identification of 
reflections by recording the incident impulse with a geophone placed in a borehole in the path of the 
incident energy and connected to an amplifier whose attenuation function is the same as that of the 
earth. 


U. S. No. 2,557,888. R. W. Olson. Iss. 6/19/51. App. 3/7/49. Assign. Geophysical Service, Inc. 
Attenuating Circuit. An avec circuit for a push-pull amplifier in which two resistors are connected 
in series with the grid of each amplifying tube and the points between the resistors returned to 
ground through opposed diodes which are biased by the control voltage. 
WELL LOGGING 
U. S. No. 2,547,875. S. Krasnow. Iss. 4/3/51. App. 10/29/36. Assign. Schlumberger Well Surveying 
Corp. 


Apparatus for Taking Physical Measurements in Boreholes. A system for remotely indicating a 
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well-logging parameter in which the parameter is made to change the current in a saturable reactor 
in an oscillator circuit in the well so as to change its frequency, and the oscillator output is converted 
to acoustic vibration of a wire line for transmitting up the hole and indicated by a vibration detector 
in contact with the line at the surface. 


U. S. No. 2,547,876. S. Krasnow. Iss. 4/3/51. App. 1/5/44. Assign. Schlumberger Well Surveying 

Corp. 

A pparatus for Investigating a Plurality of Physical Values in Boreholes. A system for remotely 
indicating two well-logging parameters in which one parameter is made to change the frequency of 
an oscillator in the well by changing the current through a saturable reactor in its circuit and the 
other parameter controls the frequency of a second oscillator in the well which modulates the signal 
_of the first, the combined signal being impressed on the wire line as an acoustic vibration which is 
picked up at the surface and separated into the two components. 


U. S. No. 2,547,950. B. D. Lee and G. Herzog. Iss. 4/10/51. App. 5/14/49. Assign. The Texas Co. 

Electrical Analogue. An electrical model for electric well logging having a wedge-shaped tank of 
electrolyte with a partition parallel to the sharp edge of the wedge, the partition being made of 
insulating material with conducting inserts so that current can flow through the partition but not 
along it. 


U. S. No. 2,549,109. J. W. McPhee. Iss. 4/17/51. App. 3/19/49. Assign. Lane-Wells Co. 
Radioactive Locating Means. A method of radioactive depth marking by using a casing collar 
having a radioactive material homogeneously incorporated in it. 


U. S. No. 2,550,004. H-G. Doll. Iss. 4/24/51. App. 12/22/43. Assign. Schlumberger Well Surveying 

Corp. 

Method of Establishing Markers in Boreholes. A method of producing a depth marker inside casing 
which is related to a radioactive marker outside the casing by running a radioactivity detector and 
an adjacent electrode and passing a large current into the casing from the electrode so as to polarize 
the casing, the polarization being subsequently located with a natural potential log. 


U. S. No. 2,550,005. H-G. Doll. Iss. 4/24/51. App. 4/22/48. Assign. Schlumberger Well Surveying 
Corp. 
Well Logging Method and Apparatus Utilizing Periodically Variable Spontaneous Potentials. A 
system for producing fluctuating pressure in a fluid-filled well by means of a reciprocating pump 
piston and observing the alternating electro-filtration potential produced at various horizons. 


U.S. No. 2,552,428. A. B. Hildebrandt. Iss. 5/8/51. App. 9/21/49. Assign. Standard Oil Develop- 
ment Co. 
Sleeve Sup ported Electrode for Well Logging. An electric-logging electrode system having a spring- 
expanded sleeve of electrically-insulating material carrying on its outer surface electrodes which are 
pressed against the formation but which are insulated from the body of drilling fluid in the borehole. 


U.S. No. 2,553,233. A. Blanchard. Iss. 5/15/51. App. 2/17/45. Assign. Schlumberger Well Surveying 
Corp. 
Method for Determining the Presence of Gas in Dispersions. A method of determining gas in drilling 
mud by subjecting the mud to supersonic vibration and measuring the potential developed between 
a pair of electrodes immersed in the mud. 


U.S. No. 2,555,350. J. F. Bayhi. Iss. 5/15/51. App. 4/17/47. Assign. Standard Oil Development Co. 

Method and A pparatus for Detecting Corrosion. A system for detecting corrosion of pipe caused by 
electrolysis which has a magnet rotating in a horizontal plane close to the inside surface of the pipe 
and also has longitudinally-spaced electrodes which contact the pipe and separately recording varia- 
tions in magnetic flux and potential with depth. 
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U.S. No. 2,554,174. H-G. Doll. Iss. 5/22/51. App. 6/13/42 and 10/8/47. Assign. Schlumberger Well 
Surveving Corp. 
Well Logging Method and Apparatus. A system for determining permeability of formations by 
forcing fluid alternately against one side of the hole and the other and measuring the alternating dif- 
ferential electrofiltration potential developed between electrodes on each side of the apparatus. 


U.S. No. 2,554,844. L. M. Swift. Iss. 5/29/51. App. 3/22/46. Assign. Well Surveys, Inc. 

Casing Collar Locator. A well-logging sonde having a non-magnetic section in its case with magnet 
pole pieces above and below it and a coil for detecting variations in flux caused by casing collars in an 
electronic circuit for superimposing the collar signals on the logging signal. 


U.S. No. 2,557,168. J. J. Arps and O. W. Ward. Iss. 6/19/51. App. 11/28/49. 

Continuous Electric Logging While Drilling. A method of electric logging while drilling in which a 
relaxation oscillator in the drill is connected to insulated electrodes and each oscillation trips a relay 
which disconnects an insulated section of the drill string, the change of resistance of the drill string 
to ground being detected at the surface. 


U.S. No. 2,557,488. C. L. White. Iss. 6/19/51. App. 3/28/47. Assign. Phillips Petroleum Co. 


Apparatus for and Method of Determining Permeability of Earth Formations Penetrated by Well 
Bores. A method of determining permeability by placing a liquid opposite the formation and a lighter 
liquid above it and following the interface down the well by means of a line and float of intermediate 
density which floats on the interface. 


U.S. No. 2,558,427. H. C. Fagan. Iss. 6/26/51. App. 5/8/46. Assign. Schlumberger Well Surveying 
Corp. 
Casing Collar Locator. A magnetic differential type casing joint locator having two permanent 
magnet circuits each including part of the casing so as to form a magnetic bridge with a pick-up coil 
between the magnets to detect unbalance caused by the joint. 


MISCELLANEOUS 


U.S. No. 2,552,890. J. D. Eisler. Iss. 5/15/51. App. 3/9/44. Assign. Stanolind Oil and Gas Co. 


Elevation Indicator. An elevation meter for use on a vehicle having a disc pendulum carried on a 
transverse horizontal axis by torsion elements and reflecting a beam of light to a photocell after re- 
flection at an angle proportional to the pendulum deflection. 


U.S. No. 2,557,021. P. S. Williams. Iss. 6/12/51. App. 9/12/46. Assign. Standard Oil Development 
Co. 


Hydraulic System for Measuring Differences in Elevation. An elevation measuring system in which 
the points are connected by a liquid-filled tube and difference in pressure observed. 


U.S. No. 2,557,158. C. F. Teichmann. Iss. 6/19/51. App. 11/30/48. Assign. Texaco Development 
Corp. 
Radioactive Exploration. A system for radioactivity surveying from an aircraft in which a direc- 
tional detector is adjusted by means of a cam to point toward a fixed spot on the ground as the craft 
flies over it and after completing the measurement is directed to a new spot. 


U.S. No. 2,557,691. F. Rieber. Iss. 6/19/51. App. 3/19/49. Assign. Geovision Inc. 


Electronic Function Generator. A system of portraying reflecting horizons on a c-r tube screen in 
which the velocity-depth function is introduced by having a photocell pick up light from an auxiliary 
c-r tube through an empirically-made absorption wedge and modulate the scanning of the portraying 
tube. 
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Science and Common Sense, by James B. Conant, Yale University Press, New Haven, 1951, 371 pp., 
$4.00. 


In this book, Dr. Conant, who is President of Harvard University as well as a distinguished or- 
ganic chemist, has endeavored to explain the methods of experimental science for the layman. Draw- 
ing largely on material which he presented over a three-year period to Harvard undergraduates 
taking a survey course in Natural Science, he addresses this book to “the intelligent citizen who as a 
voter may, to an increasing extent, be interested in congressional action on scientific matters” as 
well as to “the lawyer, banker, industrialist, government official, politician or newspaperman.” 
Although the author has intended it for those with little or no training in the sciences, the book should 
make worth while, stimulating reading for many actually engaged in scientific activity. 

If everyone who refers at one time or another to the “scientific method,” were pressed to define 
the term, there would probably be a wide range of variation in the replies. To some it is the organiza- 
tion and classification of information; to others it is the impartial and rational derivation of conclu- 
sions from facts. Still others would consider it the synthesis of many separate observations into a 
statement of ultimate reality. To Dr. Conant, the scientific method is something more specific; it is 
the formulation, as a result of experimentation and observation, of a conceptual scheme, such as a 
hypothesis or theory, which yields deductions that can be tested by further experiment and observa- 
tion. All the other definitions are inadequate because they describe equally well many other spheres 
of activity, ranging from history to the detection of criminals which would not be classified as science. 
On the basis of Conant’s definition, science is dynamic and requires continual progress from experi- 
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ment to deduction to experiment. If all laboratories were to close, science would cease to exist, no 
matter how great the accumulation of scientific information in the libraries. Science, in such an 
event, would turn into dogma. 

The author’s approach is essentially historical. He illustrates his explanations of the scientific 
method by fairly detailed accounts of classic experiments conducted by such pioneers in physics, 
chemistry and biology as Torricelli, Boyle, Lavoisier, Volta, and Pasteur. These examples are chosen 
for their simplicity since it is much easier to explain, in a book for the layman, the development 
of the concept that the atmosphere exerts a pressure than to present the origins of the quantum 
theory in a way that appeals to him. 

A chapter on earth sciences entitled “The Study of the Past” points out the methods and ways of 
thinking which are peculiar to geology, also why geology is a science and why history is not although 
both aim to give us a record of the past. A brief section on geophysics emphasizes the point that there 
is little difference between the methods of geophysics and the methods of physics. 

In discussing gravity measurements in this section, the author falls into the error of referring to 
“variations in the gravitational constant” from place to place. In a brief reference to seismic prospect- 
ing he indicates that “formations” associated with petroleum can be located by observing the speed 
of seismic waves through them. “Structures” would be a more accurate term. Both slips are minor 
and do not detract in the least from the validity of the discussion on the geologic sciences. 

The most interesting chapter to most readers engaged in commercial geophysics should be the 
one entitled “The Impact of Science on Industry and Medicine.” Here he makes the highly cogent 
point that the chief object of industrial research should be to reduce the degree of empiricism in the 
processes involved. The less empiricism required, the better is the conceptual structure by which 
predictions can be made and new applications can be deduced. The analysis here reminds the re- 
viewer of Paul Weaver’s stimulating essay, “The Relative Place of Empirical and Analytical Meth- 
ods of Geophysical Prospecting,”’ published in Geophysics in July, 1942. This broader view of the 
purpose of applied research will probably not be accepted by many executives who look for more 
tangible products or improvements to justify research expenditures. But in the long run, reduction 
of empiricism usually pays off better than the type of “research” program which is limited to more 
immediate practical objectives. 

Dr. Conant’s book is recommended to all readers of this journal as a lucid and thought-provoking 
discussion of a subject that concerns aJl geophysicists, both professionally and as citizens who must 


help establish public policy regarding the place of science in our national life. 
; M1Ton B. Dorin 





“Die geophysikalische Erschliebung des Emslandes,” by Dr. Hans Cloos. Erdél und Kohle, Hamburg, 
Germany, No. 5, Vol. 4, May, 1951. 


This article is an excellent case history describing the discovery of oil in this area and includes a 
discussion of the economic, geologic, geophysical, and, to a limited extent, philosophical factors. Its 
organization is patterned very closely to many of the case histories appearing in the Society’s Geo- 
physical Case Histories volume, which the author states in the introduction is his intent. The seventeen 
illustrations are an important factor in aiding the reader in following the detailed geologic interpreta- 
tion of the reflection and refraction seismic as well as the magnetic and gravity data because these 
are discussed more thoroughly, particularly from the stratigraphic and paleogeographic aspects, than 
in the average case history. Dr. Cloos also presents from an historical standpoint the discovery and 
development of this area in which there was no commercial interest prior to 1938; perhaps because 
there was no subsurface information. However, the original refraction work, which did not succeed in 
outlining any commercial structures, served to arouse interest in the area, and several reflection 
seismic surveys and a more detailed gravity survey were subsequently conducted. 

The early regional magnetic and gravity surveys were not adequate to be of commercial signifi- 
cance. However, academically they furnished a picture which could be used as a basis of discussion 
regarding the basement configuration and variation in thickness of the sediments. A more detailed 
commercial gravity survey that was conducted subsequently showed the presence of small anomalies 
superimposed on the regional effects. The small amount of torsion balance data available is also dis- 
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cussed and the conclusion drawn that it was of no commercial value except for outlining a few special 
trends. 

A tribute to the quality of American geophysical instruments is made in this paper because of its 
statement that when no significant seismic reflection data were obtained in the early 1940’s it was 
hoped when American equipment became available that the area could be worked successfully. 
However, when this equipment became available, although improved results were obtained, superior 
interpretation techniques were necessary to utilize the new records to the best advantage. 

The general impressions obtained from reading this article are that geophysics and geology have 
no political boundaries, that comparable problems occur on a global scale, and that the techniques 
which can be used to solve a certain geological-geophysical situation in one area will be equally 
successful in another provided similar factors are involved; also that the more subsurface geologic 


information available, the better will be the resulting interpretation of the geophysical information. 
RICHARD A. GEYER 





Practical Oil Geology (6th ed.), by Dorsey Hager, McGraw-Hill Book Co., Inc., New York, 1951, 

589 pp., $7.50. 

Dorsey Hager has revised his comprehensive volume on the relationship of geological work 
and the oil industry. In it there is a great deal of material which is reminiscent of Lahee’s classic 
volume on field geology. However, Hager deals more extensively with petroleum exploration, the 
drilling of wells, and the work of the petroleum geologist. 

As a handbook for the practicing petroleum geologist, its value is limited. It would serve very 
well for students as an introduction to petroleum geology, and for those whose work necessitates an 


understanding of petroleum geology in all its many phases. 
RoBERT W. GARWICK 


General Crude Oil Company 





Physical Geography, by Arthur N. Strahler, John Wiley and Sons, New York, 1951, 442 pp. (double 
columns), $6.00. 


Principles of Human Geography, 6th edition, by Ellsworth Huntington, revised by Earl B. Shaw, 
John Wiley and Sons, New York, 1951, 773 pp., $6.25. 


Both of these books present geography in a somewhat different light from that in which many of 
us, whose last formal contact with it was in grade school, tend to regard the subject. Neither makes 
any effort, except for purposes of illustration, to catalog information on the natural or cultural features 
of specific localities. Both are more concerned with the general principles governing the various kinds 
of environment in which man lives. Professor Strahler’s book treats the physical factors, such as 
climate and physiography, which control environment. The late Professor Huntington’s book, on 
the other hand, deals with the dependence of human activity and development upon environment. 
Taken together, the two books cover both the cause and the effect of the environmental conditions 
found in areas of human habitation. 

Physical Geography is an exceptionally attractive, well-written treatise on the physics of the 
earth’s surface. It is intended for college students who have had little or no background in physical 
science or mathematics. As Professor Strahler points out in the introduction, physical geography is a 
body of basic principles selected from geodesy, astronomy, cartography, geomorphology, geology, 
hydrology, meteorology, botany, and oceanography. There are few persons, even among the techni- 
cally trained, who are so familiar with all these branches of earth science that they cannot obtain a 
better understanding of many natural phenomena from an elementary, non-mathematical presentation 
such as this. Although most of the sciences from which physical geography is derived are divisions of 
geophysics, the average geophysicist, like most others in scientific work, will probably find much 
outside his particular specialty that is new to him in a book of this type. 

Strahler’s book is divided into three parts. The first deals with topics relating to the earth’s 
shape and rotation, such as illumination, time, and tide. The second deals with maps, landforms, and 
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geologic structures affecting topography. The emphasis given this phase of the subject reflects the 
author’s primary field of interest, which is geomorphology, and this section alone could constitute a 
good beginning text in this field. The third part covers weather, climate, vegetation, and soils. The 
photographs and line drawings are unusually numerous, well planned, and effectively reproduced. 
Full use is made of maps and block diagrams, particularly in the section dealing with geomorphology. 

The Huntington book considers geography not so much a physical science as a social science 
although such topics as weather, climate, maps, and physiography are treated from the physical 
standpoint. The scope of the book can best be conveyed by listing the headings for its ten major 
divisions: Man’s Geographical Relationship, Man’s Relation to the Earth as a Globe, The Funda- 
mentals of Climate, Man’s Relation to Land Forms, Physiography and Human Progress, Man’s 
Relation to Soil and Minerals, Man’s Relation to Climate, Man’s Regional Relationships, Man’s 
Relation to Man, and Countries of the World. Emphasis is given throughout to the explanation of 
differences between various racial, national, and regional groups on the basis of geographic distinction. 
The last part of the book, Countries of the World, summarizes the histories of the various nations in 
terms of their individual geographic influences. 

Although on the surface there might not appear to be any close connection between human 
geography and geophysical prospecting, there is actually an intimate relationship, especially as re- 
gards the planning of field operations. To plan a seismic or gravity survey in a new area requires 
knowledge of climatic factors, topography, transportational facilities, and local availability of food 
and fuel. While this book could not in itself serve as an adequate source of such information for all 
areas in which geophysical operations might be intended, it should focus attention on the many 
problems that inevitably arise in any activity where environment is as controlling a factor as it is in 
prospecting. The book contains a fairly complete bibliography, with titles arranged according to 
location, which should be useful when problems arise in specific areas. 

Mitton B. Dosrin 





Geography of the Pacific, edited by Otis W. Freeman, John Wiley & Sons, New York, 1951, 573 pp. 


If anyone wishes to know the elevation and the square area of certain Pacific islands, or other 
types of factual data on population, commerce, and climate, this book is recommended highly for 
these purposes. In short, it is a compendium of facts, figures, photographs and maps written by 
thirteen specialists and is suitable as a reference book for all age groups from junior high school 
through college and adult education. However, it is not a mature “geography book” from the stand- 
point of discussing and integrating the complex inter-relationships among the physical, human, com- 
mercial, and political factors which comprise the mature approach to a study of geography and 
geopolitics. This fact was particularly apparent to the reviewer because of having recently reviewed 
such a “geography book,” Geography in the Twentieth Century, in the July 1951 issue of Geophysics. 

Admittedly, it is difficult to have continuity or a general theme in a book written by a dozen or 
more people; nevertheless, absence of such does not facilitate the reading of this book from cover to 
cover. In addition, all the contributors are from the United States, and the areas involved include 
countries which have their own competent geographers whose contributions might have resulted in 
achieving a more international and authoritative quality. The criterion for selecting some of the areas 
included also is not quite clear. For example, why should two chapters be devoted to Australia and 
none to Japan and Formosa? Furthermore, none of the islands off the western shores of South and 
Central America is represented. 

A rather comprehensive bibliography is found at the end of each section. However, the almost 
complete absence of any documentation using these references is striking. A brief description in 
general terms of the geology of most of the areas discussed in this book is presented and should be of 
interest to geologists and geophysicists. The sections on meteorology, climatology, and topography 
for each area should also be useful to members of these professions when they have occasion either to 
visit or prepare exploration programs for any of the countries discussed in this book. 

RIcHARD A. GEYER 

































REVIEWS 711 


Einfiihrung in die Physik der festen Erde (Introduction to the Physics of the Earth) by Dr. Hellmut 
Berg, published by S. Hirzel Verlag, Stuttgart, 1949, 295 pages. 


This book represents another attempt to condense the information contained in some of the 
available compendia and periodical literature of geophysics. Admittedly this is one of the objectives 
as stated by the author in the preface but it is difficult to see why cognizance has not been taken of 
the progress made in both pure and applied geophysics as described in the literature of the past ten 
years or so. There is very little to be found in this volume which does not appear in books published 
previously of this type such as, for example, Tams “Grundziige der physikalischen Verhaltnisse der 
festen Erde” (“Fundamentals of the Physical Conditions of the Earth”) which appeared in two vol- 
umes, 1932 and 1037. 

The section on instruments in the chapter on geomagnetics stops with the Schmidt vertical bal- 
ance. No mention is made whatever of the airborne magnetometer which should prove to be a major 
factor in increasing our knowledge of the geomagnetic characteristics of the earth not only from the 
prospecting standpoint but from pure geomagnetism as well. In the section on the temperature of the 
earth the publications of Francis Birch and his contemporaries are conspicuous by their absence. In 
seismology the work of Ewing and his colleagues in the field of marine seismic studies is also ignored, 
as are the theoretical advances in wave propagation made by Pekeris. 

One of the advantageous features of books of this type which have appeared in the past such as 
Tams; Meisser “Praktische Geophysik” (“Applied Geophysics”) 1943; and Jung “Angewandt Geo- 
physik” (“Applied Geophysics”) 1948, has generally been the presence of extensive bibliographies 
which permit the reader to augment the information available in these digests. However, even this is 
missing and the bibliography; which includes “important geophysical periodicals,” occupies less than 
a page. Upon examining the list of periodicals cited, a reason for some of the deficiencies mentioned 
above may be postulated. For example, such periodicals as Transactions of the American Geophysical 
Union, Geophysics, and the Geophysical Supplement of the Monthly Notices of the Royal Astronomical 
Society do not appear. 

RicHARD A. GEYER 





The Tectonics of Middle North America (east of the Cordilleran System), by Phillip B. King, Prince- 
ton University Press, Princeton, New Jersey, 1951, 203 pp., $3.75. 


This latest publication of Phillip King arose from a “need for a modern description of the tectonic 
features of the United States for classroom use and as an explanation of the tectonic map of the 
United States” as stated in the author’s foreword. The book is therefore essentially a textbook, at 
about senior collegiate level, and primarily useful in elucidating the features of the United States 
tectonic map. The reader should keep the author’s purpose constantly in mind in order to evaluate 
the contents of the volume properly. 

Mr. King is thoroughly qualified to write this book. He has been a geologist with the Geological 
Survey since 1930 and has made well-known contributions to geologic progress wherever he has 
concentrated his efforts. His work in West Texas is especially familiar to petroleum geologists. More- 
over, he has already published ‘An outline of the structural geology of the United States” and he 
edited the “‘Tectonic Map of the United States.” 

The volume is divided into sections according to geologic provinces. The Central Stable Region, 
the Interior Lowlands, Paleozoic structures south and southeast of the Central Stable Region, and 
the Coastal Plains are discussed. In each stratigraphy and structure are considered. However, there 
is a distinct difference in approach to the writing of the section on Paleozoic structures as compared 
to the other sections. The Paleozoic structures are not only described but an attempt is made to ar- 
range the geologic processes scientifically. That is, a sequence of events in time and space is de- 
veloped so that prediction, part of the scientific method, can be made of structural and sedimenta- 
tional conditions in unknown areas. This geologic philosophizing is essentially absent in the other 
three sections which are restricted to descriptions of tectonic features. As a consequence, approxi- 
mately one hundred pages are devoted to Paleozoic structures and thirty to each of the three remain- 
ing sections, 
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The contributions of the three smaller sections on the Craton, Interior Lowlands, and Coastal 
Plains are therefore confined to the amplification of the tectonic map. These sections achieve that 
purpose. They may result in appearing fragmentary to a reader expecting a comprehensive discussion 
of the areas. In some instances, little geologic investigation has been carried out in the area, and the 
picture presented must rely on meager evidence. In other cases, a great amount of data is available 
but only the apparently significant results can be reported in a book of this type. Such spot sampling 
can lead to misconception. Furthermore, each bit of evidence must be thoroughly reviewed because 
so much importance is attached to it in a book of this type. For example, in the discussion of the Gulf 
Coastal Plain the statement is made that thinning gulfward of the Tertiary deposits is suggested by 
seismic and gravity observations. The late review of the geophysical evidence by the 1950 G.S.A. 
symposium on Gulf Coast sedimentation did not offer any such suggestions. The axis of the Gulf 
Coast geosyncline is still undefined. 

The section on Paleozoic structures is the most interesting to read because it does go beyond the 
mere descriptive phase. In general, King tries to synthesize the sedimentation in a geosyncline, the 
structure in a tectogene, and the accompanying igneous activity into a unified geologic process. He 
attacks this p-oblem by dividing the Paleozoic structures into four salients, Newfoundland, Northern 
Appalachians, Central and Southern Appalachians and the Wichita-Ouachita systems. He studies 
each one relatively closely and finds in each supporting evidence for the time and space relations 
necessitated by his thesis. Instead of attempting to present this complete thesis here, some incidental 
points of interest to the reviewer are noted below. 

(1) Strong evidence is presented against the concept of simultaneous world orogenies, relatively 
sudden and of short duration. Instead, the case for more-or-less continuous deformation at different 
times and in different provinces is strongly defended. 

(2) Apparently isolated basins, such as the Narragansett Basin, found within orogenic belts are 
developed as perhaps the latest phase of the orogeny. 

(3) The role of the Gulf of St. Lawrence is very interesting especially to those connected with the 
petroleum industry. The counterparts of the structures on Newfoundland are found displaced much 
to the westward on the mainland so that the Gulf of St. Lawrence ‘‘must be underlain by major 
structural features of obscure character.” Evidence for the displacement is also found in the sub- 
marine topography in the vicinity of the Grand Banks. King points out the striking similarity of this 
structural linkage with that between the Southern Appalachians and the Ouachita structures. 

(4) The inclusion of the Ouachita as an Appalachian salient leads to the classification of the 
Ouachita sediments as eugeosynclinal and the Wichita sediments as miogeosynclinal. However, the 
Wichita structure system is at right angles to the Ouachita structural axis. This is explained, of 
course, by the thrusting of the Ouachita system over the Wichita system. But in every other case the 
eugeosynclinal and miogeosynclinal structures are parallel. All in all, the evidence for the thesis 
proposed by King is most abundant in Newfoundland and becomes less and less as one goes south. 
In the Ouachita area, the evidence is indeed so meager that the ideas proposed are controversial. 

As stated in the béginning of this review, this book will undoubtedly become a standard textbook 
for senior collegiate geological courses. This reviewer sincerely welcomed the opportunity to take a 
refresher course in the geologic structural history of middle North America, and the refresher 
course was especially welcome because the material is presented very concisely and logically. A 
maximum amount of information is gained with minimum effort, and that alone is sufficient recom- 
mendation for reading this volume. 

NELSON C, STEENLAND 





“Radium and Deep-Sea Chronology,” by Hans Pettersson, Nature, v. 167 June 9, 1951, p. 942. 


This letter to Nature, coming from the Oceanografiska Institutet in Géteborg, Sweden, gives 
interesting new results on the sedimentation problem in the sea. The Swedish Deep-Sea Expedition 
has recently obtained new data on the concentration of uranium and radium in sea water. These 
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data confirm earlier results of a large deficiency of radium in equilibrium with uranium. The author 
had previously suggested this deficiency might be explained by the precipitation of ionium, the parent 
of radium. This precipitation would also account for the high concentrations of ionium and radium 
in red clay and some other surface deposits in the sea. By computing the amount of ionium precipi- 
tated and the concentration of ionium in the bottom sediment, the author has estimated the sedimen- 
tation rate to give the necessary dilution for this latter value. The sedimentation rate thus computed 
is o.5 mm per thousand years for a core taken near Tahiti. It is suggested that this figure can be 
changed to 0.4 due to the diffuseness of the radium values from which the ionium concentration in 
the core was computed. The sedimentation rate is also computed by obtaining the decrease in ionium 
concentration with depth from the radium measurements. If this decrease is assumed to be due only 
to radioactive decay, the age versus depth curve is obtained for the sediments since the half life, or 
rate of disintegration, for ionium is known. The rate of sedimentation by this older method agrees 
with the above value of o.4 mm per thousand years. Reports of the Swedish Deep-Sea Expedition 
will carry a more detailed account of this work. 
Lynn G. HOWELL 





“Measurements of Gravity in Ireland,” by A. H. Cook, Department of Geodesy and Geophysics, 
Cambridge University, Geophysical Memoirs No. 2, Part 1, 1950. 


Prior to the measurements described by the author, no observations of gravity had been made 
in Ireland. Measurements in Wales had shown that gravity increased toward the Irish Sea, whereas 
submarine observations in the Irish Sea gave somewhat lower values, and it was of interest to know 
how the gravity field of Ireland was related to that of Wales and the Irish Sea. In addition the values 
of gravity would be useful for geologic and geodetic purposes. 

The apparatus used was the Cambridge pendulum apparatus designed by Sir Gerald Lenox- 
Conyngham which has been described in other publications (Lenox-Conyngham, 1928; Bullard, 
1933; Bullard and Jolly, 1936; and Bullard, 1936). Two pendulums of period about 1.01 second swing 
together with equal amplitude and opposite phase so that sway disturbance is eliminated and ground 
movements do not perturb the mean of the two periods. To reduce damping and so enable small 
amplitudes to be used, the pendulums swing in an air pressure of about 30 mm Hg. The method of 
calculating the results is described in detail in this paper. In all, four stations were occupied by the 
pendulum apparatus. 

While the pendulum observations were being made Mr. H. I. S. Thirlaway ran a gravimeter 
traverse around the four stations, a distance of approximately 400 miles. Difference in gravity be- 
tween stations ranged from about 22 mgal to 143 mgal with the probable error being about +0.5 mgal. 
A report on the gravimeter survey is to be published as Part 2 of these memoirs. It is interesting to 
note that the standard deviation for both the pendulum and gravimeter is 0.5 mgal. It is regrettable 
that Mr. Thirlaway did not have a modern gravimeter for this survey. With such equipment this 
probable error would have been reduced by a factor of five or even ten. 

The final results which are listed below are derived from data taken with the pendulum and the 
gravimeter. These results are expressed as differences of gravity from Cambridge, and using least 
squares they are: 


Dublin, Dunsink Observatory +120.84+0.56 
Sligo, Courthouse +197.28+0.82 
Galway, University College + 96.21+0.64 
Cork, University College — 22.51+0.87 


The value of gravity at Cambridge (Pendulum House) has been taken as 981.2650 cm/sec?. No 
probable error for this value is given. 
An excellent bibliography is included with this paper. 
EUGENE FROWE 
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* Proceedings of the Cambridge Philosophical Society, Vol. 47, No. 3 (July, 1951) 

*Quarterly Journal of the Geological Society of London, Vol. 106, No. 2 (May 18, 1951) 

*La Ricerca Scientifica, Vol. 21, No. 3 (Mar., 1951) 

Review of Scientific Instruments, Vol. 22, Nos. 4, 5, 6 & 7 (Apr.—July, 1951) 

*Revue de I’ Institut Frangais do Petrole et Annales des Combustibles Liquides, Vol. 6, No. 5 (1951) 

*Revista di Geofisica Applicata, Vol. 12, No. 1 (1951) 
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Books 
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Basin and Southeast New Mexico, Permian Basin Geophysical Society, Midland, Texas, 1951 
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HISTORICAL NOTE 


Dr. Joseph A. Sharpe has called attention to a paper of considerable historical interest which 
appeared in Volume IV (1920) pp. 83-85 of the Bulletin of the American Association of Petroleum 
Geologists. Dr. Sharpe states, ‘“This appears to be the first very clearly stated suggestion concerning 
the use of the seismic reflection method for the mapping of underground structure.” The paper, en- 
titled “Suggestions of a New Method of Making Underground Observations,” is by Dr. J. A. Udden 
and is here reprinted in its entirety. 


Bulletin Amer. Assoc, Pet. Geologists IV, pp. 83-85. 1920. Paper read at Dallas meeting March, 
18-20, 1920. 


BULLETIN OF THE AMERICAN ASSOCIATION OF PETROLEUM GEOLOGISTS 
SUGGESTIONS OF A NEW METHOD OF MAKING 
UNDERGROUND OBSERVATIONS 


J. A. UpDEN* 


In the north-central part of the State of Texas, where wells are being drilled to the depth of 4,000 
feet and more in order to penetrate the Bend formation (Pennsylvanian), which is the chief oil-bearing 
horizon in that part of the State, it has been found that the structure of the Bend and that of the 
underlying Ellenburger limestone (Ordovician) is quite generally the same. The plotting of under 
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FicureE 1. Reflection of earth waves, started at surface, from the contact between the Bend 
formation and the underlying Ellenburger limestone. E, point of explosion; S, position of seismo- 
graph; R, point of reflection of earth wave; RD, depth from surface to Ellenburger limestone. 


ground structure by drillers’ logs or by observations on cuttings from the Bend is very difficult, for 
the reason that there are small variations in the rocks of the Bend, not only vertically but also hori- 
zontally, and that these small variations lie near the limit between shale and limestone. Whether a 
certain part of the formation is reported as shale or as limestone will in large part depend upon the 
judgment and experience of the driller, or other observer. It has been exceedingly difficult to trace a 
definite horizon in the Bend with any considerable certainty. Exposures of the Bend indicate also 
that there are horizontal changes in this formation. For those who use drillers’ logs in the making 


* Dr. J. A. Udden, Director, Bureau of Economic Geology and Technology, Austin, Texas. 
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of structural maps, it has been found expedient to use the contact of the Bend and the underlying 
Ellenburger as the key horizon. Contours are drawn on the upper surface of the Ellenburger. 

The top of the Ellenburger is a plane marking the greatest change of all the sediments in this 
part of the State. The Ellenburger itself is about 1,000 feet thick, where it is present in its full de- 
velopment, and almost all of this is one solid limestone. The Bend, which overlies it, is, as already 
stated, quite variable in its composition. In the San Saba country it has considerable limestone. West 
from here it becomes more and more shaly, and the same change is noticeable in a northerly direction, 
at least at some points. The Bend is overlain by the Strawn, Canyon and Cisco divisions of the 
Pennsylvanian. These consist largely of shale with some sandstones and limestones, and may collec- 
tively be characterized as much less uniform in composition than the Ellenburger. In the course of 
development in the Ranger-Desdemona field and other fields in this part of the State, it has become 
evident that the structure of the Bend does not always parallel the structures in the outcropping 
rocks and that therefore the prospecting of the Bend by the common geologic methods is not entirely 
satisfactory. 

It has suggested itself to the writer that it ought to be possible, with present refinements in 
physical apparatus and their use, to construct an instrument that would record the reflections of earth 
waves started at the surface, as they encounter such a well-marked plane of difference in hardness 
and elasticity as that separating the Bend and the Ellenburger formations. 

Some such instrument as a seismograph, especially constructed to be sensitive to vertical waves, 
possibly a slight adaptation of some of the present seismographs would be amply suitable for the 
purpose. A seismic wave might be started by an explosion at the surface of the earth, and a record 
of the emerged reflection of this wave from the upper surface of the Ellenburger limestone might be 
registered on an instrument placed at some distance from the point of explosion. The record would, 
of course, be a component of the direct and the reflected waves. It ought to be possible to notice the 
point at which the first reflection from the Ellenburger appears on a continuous record. The time of 
the record being known, it ought to be possible to figure the depth to the plane of reflection. It 
seems to me that by the making of a number of observations at points where the distance to the upper 
surface of the Ellenburger is known, it would be possible to learn to interpret the record correctly. 
It might be found necessary to make new adaptations on the working parts of the instrument. The 
difficulties ought not to be insurmountable. After these are once mastered, observations could be made 
anywhere, probably with quite rapid progress, and the entire upper surface of the Ellenburger might 
be mapped over wide areas with relatively small effort. The topography of the Ellenburger could, so 
to speak, be contoured directly by any number of observations desired. With such a topographic 
map of the surface of the Ellenburger, it seems to me that millions of dollars’ worth of drilling could 
be eliminated, which will otherwise have to be done before all of the pools in the Bend shall have 
been located. 

Why not pay twenty thousand dollars for “theoretical” experiments, if by so doing a million in 
expenses on “practical’’ work can be saved? 


GEOPHYSICS IN HIGHWAY CONSTRUCTION 


Dr. Sigmund Hammer has called attention to a very important paper relating to geophysics as 
used in highway construction. The paper is ‘Development of Geophysical Methods of Subsurface 
Exploration in the Field of Highway Construction” by R. Woodward Moore, Highway Engineer, 
Physical Research Branch, Bureau of Public Roads. It is published in Soil Exploration and Mapping 
sssued as Bulletin No. 28 (1950) of the Highway Research Board, and it is obtainable through the 
Bureau of Public Roads, Department of Commerce, Washington 25, D. C. 


PROFESSOR SYDNEY CHAPMAN’S LECTURE 


The presence in the United States during the academic year 1950-51 of Professor Sydney Chap- 
man of Oxford University permitted the Distinguished Lectures Committee of the Society to secure 
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the services of this outstanding authority on the magnetism of the earth as lecturer before most of 
the local sections of the Society. 

During the eight day period, March 26 through April 2, Professor Chapman delivered his lec- 
ture, entitled “Survey of the Present Problems of the Magnetism of the Earth,” tomembets of the geo- 
physical societies and their guests at Denver, Tulsa, Houston, Dallas, Fort Worth, and Shreveport; 
conferred with individual geophysicists in these cities; and visited geophysical laboratories and manu- 
facturing plants in Denver, Tulsa, Houston, and Dallas. - 

In the following condensation of the lecture, prepared from a written text provided by Professor 
Chapman, particular emphasis has been given to Professor Chapman’s views concerning the present 
théories of the earth’s magnetism. 


SURVEY OF THE PRESENT PROBLEMS OF THE MAGNETISM OF THE EARTH 
Historical 


The Greeks first observed magnetic and electric attractions, but it was not until the twelfth 
century that manuscripts noted the directive property of the magnet. The French soldier-scientist, 
Petrus Peregrinus, first described the dipolarity of the magnet in a famous letter. Petrus believed, 
however, that the compass pointed towards the true pole of the heavens. A Portuguese, Juan de 
Castro, made detailed observations of the earth’s magnetic declination in 1543-45. The magnetic 
dip was discovered by Robert Norman in 1576, and he concluded that the controlling influence was 
within the earth. 


The Earth a Magnet 


In 1600 William Gilbert published his treatise de Magnete in which he expressed his discovery 
that the terrestrial globe is itself a great magnet. Gilbert thus fathered the science of terrestrial 
magnetism. He clearly recognized the big problem: why is the earth a greai magnet, and the secondary 
problem why is tts magnetism not quite regular. 


The Secular Variation 


Within thirty years, Gellibrand had discovered the secular magnetic variation, and a few decades 
later Edmond Halley, the astronomer, proposed a theory for secular variation. He pictured a fantastic 
earth formed of an inner core surrounded by one or more shells, each magnetized differently, rotating 
at different rates about different axes. 


Magnetic Maps 


To check his theory, Halley made the first scientific voyage—a magnetic survey of the North 
and South Atlantic Oceans. On his return he published the first isogonic chart, and a few years later 
Halley produced a world magnetic chart. About 1890 van Bemmelen, using early naval records, con- 
structed declination charts for epochs from 1500 A.D. onwards at 50 year intervals. 


The Dipole Field 


By 1836 the French Admiral Duperrey had drawn maps showing the lines of horizontal magnetic 
force on the earth’s surface, and by that time Poisson had developed the theory of magnetic potential. 
The great German mathematician Gauss confirmed mathematically that the earth’s magnetic field 
proceeds from within, and that, to a first approximation, the field is that of a uniformly magnetized 
sphere, or a magnetic dipole, at outside points. Gauss showed for the first time how to measure the 
magnetic intensity in terms of units of length, mass, and time, and he set up the first magnetic ob- 
servatory. 


The Non-Dipole Field 


Although the earth’s field is approximately a simple dipole field, divergences from this are im- 
portant enough to form, together with the main field and the secular variation, the third main problem 
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of the subject. The non-dipole field, namely the vector difference between the earth’s field and the best 
fitting dipole field, can be represented by maps which show its regional character and complexity. 


Tsoporic Maps 

The secular change of the earth’s field can likewise be represented on maps by lines of equal rate 
of change, or isopors. At present, the regions of decreasing intensity exceed those of increasing in- 
tensity, and iv: the last century the dipole field has decreased by about five percent. 


Local Anomalies 

Some of the anomalies of the earth’s field are certainly due to magnetic deposits not far from the 
surface of the crust. But other major regional anomalies are present which gradually change their 
force and location, and it is difficult to believe these can be explained by crustal magnetism. 


Paleomagnetism 

The researches of the Department of Terrestrial Magnetism of the Carnegie Institution of 
Washington have extended the time-span of our knowledge of the earth’s magnetism. The results 
obtained are consistent with the idea that for 50 million years or more the earth’s magnetic polarity 


has had an average direction nearly coinciding with that of its axis of rotation. These conclusions as . 


to the stability of the general direction of the earth’s magnetization are supported also by similar 
researches on sediments taken in cores from the bed of the ocean. A contrary indication is given by 
some volcanic dikes in South Africa, Northern England, and elsewhere. These dikes show strong 
magnetization of a sign opposite to that of the present day field, and the apparent contradiction is as 


yet quite unexplained and mysterious. 


Rotational Theories 

Kelvin in 1883 speculated on the cause of the earth’s magnetism, and concluded that it must be 
connected with its rotation and its bigness. But even rapidlv rotating bodies, in the laboratory, show a 
degree of magnetization so small that the effect defied detection until discovered by the delicate 
experiments of S. J. Barnett some twenty-five years ago. This “Barnett effect” is too small to explain 
the earth’s field. exo 

Other rotational hypotheses of a fundamental character have been proposed to explain the earth’s 
dipole component along the geographic axis. For example, the centrifugal force caused by the earth’s 
rotation will cause the planetary electrons in the atoms of the earth to lie farther from its center, 
on the average, than the positively charged nuclei, thus producing a net electric field. But this effect 
fails to account for the earth’s field by more powers of ten than the Barnett effect. 


Solar and Stellar Magnetism 

The discovery of magnetism in the sun by Hale at Mt. Wilson, first in sunspots, and then the 
general magnetic field, gave added stimulus to attempts at explanation of cosmic magnetism, par- 
ticularly as the sun’s field was related to its rotation in the same way as for the earth. For example, 
Blackett’s hypothesis implied that the rotation of electrically neutral matter would produce a mag- 
netic field as if the matter were negatively electrically charged and by its rotation would constitute 
an electrical current system. Bullard pointed out that if this were so, the outer layers of the earth 
would contribute their share to the field of these non-existent but magnetically effective quasi- 
electric currents. Therefore, downward through these layers the earth’s field would, on this hypothesis, 
change in a manner different from that usually supposed (the inverse cube law). Specifically, calcula- 
tion showed that the vertical force, not far below the surface, should still increase in the same way 
on the new hypothesis as on the old, but that the horizontal force should at first decrease downwards, 
instead of increasing. Blackett actively promoted measurement of the earth’s field in mines, and it 
gradually appeared that the hypothesis would not work. The experiments led to the definite proof 
that the outer non-magnetic layers of the earth do not contribute to the field. Meanwhile the funda- 
mental type of theory of magnetization by rotation received a hard blow from Babcock’s discovery 
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of a star with a polar field of several thousand gauss which reverses its magnetization almost sym- 
metrically every nine days. 


Electric Currents as the Cause of the Earth’s Field 


Another possible explanation of the earth’s magnetism ascribes it to electric currents. Currents 
in bodies of large size decay slowly, but even so the half life of a current system such as could produce 
the earth’s field is short compared with the age of the earth. Hence, if electric currents produce the 
earth’s field, they must be continuously maintained. Larmor in 1919 suggested that electric currents 
such as seem likely to be the cause of the earth’s and the sun’s magnetic fields may be due to dynamo 
action by convective motion in their fluid interiors in the known presence of their magnetic fields. 
In the earth’s liquid core (the diameter of which is 0.55 times that of the earth), convection is likely 
to be produced by radioactive transformations even if the concentration of radioactive matter there 
is much less than the observed amount in the crust of the earth. Motion of matter in a magnetic field 
necessarily sets up electromotive forces which impel electric currents if the matter is conducting. 
Larmor suggested that convection in meridian planes might thus set up electric currents around the 
earth’s axis which would maintain the component of the dipole field along the axis of rotation, or even 
cause it to grow from small beginnings such as the Barnett effect would produce, precisely as in a 
self-exciting dynamo. Cowling, however, has shown that convective motion in meridian planes, 
symmetrical about the magnetic axis, cannot maintain the axial dipole field. 

During the last decade, Elsasser has devoted much effective effort to the theoretical study of the 
earth’s magnetic field, and more recently Bullard has made further advances, partly along lines con- 
sidered by Elsasser. The present position seems to be that the component of the main dipole field 
which lies along the earth’s geographical axis is still not definitely explained, but that there is hope of 
explanation in a dynamo theory on the basis of an unsymmetrical convective motion. 


The Non-Dipole Field and Secular Variation 


The secular variation, though slow by human reckoning, is extremely rapid compared with the 
geological changes in the earth’s crust. It points to changes in the earth’s interior which are so rapid 
as to be most naturally ascribed to the liquid core, where there is sufficient fluidity to permit rela- 
tively rapid change. In nature most motions are to some extent turbulent, and it may be supposed 
that any large scale convection in the core is accompanied by smaller scale eddies superposed thereon. 
These will have their own additional dynamo effect, superposing smaller scale systems of electric 
current on the main system. The non-dipole field may be taken as indicating the presence of such 
subsidiary electric current systems in the core; and as the eddy motion and the resulting electric 
currents will change, this will produce the secular variation. It seems likely that the effective eddies 
and their electric current systems are those that lie near the outer surface of the core, because chang- 
ing magnetic fields produced deep in the core would induce secondary electric currents in the outer 
layers, preventing these fields from extending outside to any appreciable degree. 

By spherical harmonic analysis Vestine has approximately indicated the distribution of electric 
currents near the surface of the core that could be responsible for the earth’s surface field. McNish has 
shown that the non-dipole field could be produced by about 14 dipoles, the moment of each being 
approximately one-hundredth that of the main field, situated near the surface of the core. Their dipole 
moments indicate something as to the intensity of the equivalent electric current systems, each ex- 
tending over about one-fourteenth of the surface of the core. Bullard has estimated the dipole moment 
and electric current eddy that could produce the rapid secular magnetic variation over and near the 
southernmost part of Africa. 

The electric current eddies and conductivity of the core being thus estimated, the electromotive 
forces may be inferred. According to the dynamo theory, these depend on the product of the speed of 
motion in the convective eddies and of the magnetic intensity in the core. According to the inverse 
cube law, the magnetic intensity at the surface of the core should be about 4 gauss. If this is also 
valid inside the core, the eddy speeds required are of the order of a meter an hour. This, however, is 
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only the component speed perpendicular to the magnetic field, and there may be a considerably 
greater component along the field. 


The Toroidal Field Inside the Core 


There is some reason to suppose that the magnetic field inside the core is perhaps ten times as 
intense as that just outside, and that its character also is different from that inside a uniformly 
magnetized sphere, for which the lines of force (inside and outside) lie in meridian planes. It is thought 
that inside the core the lines of force are much twisted around the axis, so that the field has a con- 
siderable “toroidal” component around the axis, which does not extend at all outside the core. This 
twisting of the lines of force is attributed to a slight differential rotation of the core relative to the 
solid outer mantle of the earth, whose rotation is being very gradually slowed down by tidal friction. 
In a medium as electrically conducting as the core is supposed to be, the lines of magnetic force re- 
main always attached to the same elements of the fluid, and therefore share in the distortion that a 
differential rotation would produce. A convective eddy motion would also twist the lines of force, 
increasing the intensity of the magnetic field at the expense of the kinetic energy. 


Conclusion 


At the present time the theories of the earth’s main magnetic field, its major anomalies, and secu- 
lar variation are in a decidedly tentative state. The two latter are attributed to changing electric 
current systems maintained by the dynamo action of changing hydrodynamic eddies near the surface 
of the earth’s core. The main field is thought, but not yet proved, to be produced by a main system 
of convective motion, not confined to meridian planes, in the earth’s core. The westward movement 
of the secular variation and the decrease of the earth’s dipole moment during the past century are 
thought to be transient features. 


Plea to Exploration Geophysicists 


Exploration geophysicists can assist the progress of the “‘pure’”’ science of geomagnetism ina 
number of ways. For instance, they can deposit magnetic data with geomagnetic institutions or 
libraries, or they can publish such data in summary or outline. In addition, they can use their tech- 
nical experience in the design of instruments to aid, by advice and discussion, in the much needed 
modernization and re-equipment of magnetic observatories. 





























CONTRIBUTORS 


W. B. Acocs received the B.S. in E.M. majoring in 
Mining Geology in 1934, the M.S. and Ph.D. in Physics in 
1944 and 1946, respectively, all from Lehigh University. 
During 1943 and 1944, he was Instructor in Physics and 
Geophysics at Lehigh University. From 1934 to 1940, he 
was a geophysicist for the Gulf Research and Development 
Company. The major part of this period with G.R.&D.C. 
was spent in Venezuela as Party Chief—Geophysicist on 
research refraction and reflection seismograph parties. 

In 1945, 1946 and 1948, 1949, Dr. Agocs was Chief 
Geophysicist for the Venezuelan Atlantic Refining Company. 
During 1947 and 1948, while at Lehigh University, and 
during 1950 he was a geophysical consultant. Since Septem- 
ber 195c, Dr. Agocs has been Associate Professor of Physics 
and Geophysics at the University of Tulsa, and a Geo- 
physical-Geological Consultant. 

Dr. Agocs is a member of the Institute of Radio Engi- 
neers, American Institute of Physics, Sigma Xi, and So- 
ciety of Exploration Geophysicists. 





W. B. Acocs 


KENNETH E. Bure received his B.S. degree in Electri- 
cal Engineering from the University of Texas in 1926. Upon 
graduation he was employed by The Geophysical Research 
Corporation to do experimental geophysical work. Mr. Burg 
was computer and observer on one of the original experi- 
mental reflection crews and was in charge of several electri- 
cal prospecting crews prior to the time he was made Party 
Chief on a large refraction crew in 1928. 

In 1930 he resigned to join the newly-formed Geo- 
physical Service Inc. As Party Chief with that organiza- 
tion, he worked in the Gulf Coast, Mid-Continent, and 
Pacific Coast areas, and in Canada. He was placed in 
charge of the Southern Louisiana and Mississippi areas as 
District Supervisor in 1936, where he was instrumental in 
developing the techniques and equipment used in that area. 
In 1941 he was placed in charge of the Research Depart- 
ment in the Dallas Laboratories of Geophysical Service 

KENNETH E. Burc Inc. 
Mr. Burg resigned in 1942 to join the Stanolind Oil and 
Gas Company as Seismograph Technician. In 1943 he was made head of the Seismograph Department, 
and in 1944 promoted to Geophysical Supervisor in charge of all the geophysical operations of the 
Stanolind Oil and Gas Company. He rejoined Geophysical Service Inc. in January 1947, as consulting 
Geophysicist. In t950 he was promoted to Vice President, Technical, and in charge of equipment, 
field and office techniques and research. 

Mr. Burg is a member of the Society of Exploration Geophysicists, American Association of 
Petroleum Geologists, American Geophysical Union, Seismological Society of America, The Dallas 
Geophysical Society, Institute of Radio Engineers, and Tau Beta Pi Honorary fraternity. He has 
served on the Research Committee of the American Association of Petroleum Geologists, and has 
served the S.E.G. as a member of the regional program committee, as chairman of the Standing 
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Commitee on Membership, in the preparation of the Cumulative Index and as Secretary-Treasurer 
of the Society. He has also served as Vice-President of the Dallas Section and Chairman of the Pro- 
gram Committee for the first Dallas Regional Meeting. 












































SyDNEY CHAPMAN studied engineering at Manchester 
University, mathematics at Cambridge, was chief assistant 
at the Royal Observatory at Greenwich, Fellow and Lec- j 
turer at Trinity College, Cambridge, and professor of 
mathematics and natural philosophy at the Universities of 
Manchester, London and, at present, Oxford. Professor 
Chapman is a Fellow of the Royal Society of London, For- 
eign Member of the National Academy of Sciences, Wash- 
ington, President of the International Association of Ter- 
restrial Magnetism and Electricity, and in previous years ( 
served the Royal Meteorological Society, the London 
Mathematical Society, the Royal Astronomical Society, 
and the Physical Society of London as President. After a 
year as Research Associate at California Institute of Tech- 4 
nology, he appeared before the Denver, Tulsa, Dallas, Fort 
Worth, Shreveport, and Houston Local Sections in a dis- 
tinguished lecture tour sponsored early this year by the 
SyDNEY CHAPMAN Society of Exploration Geophysicists. 














J. Harvey Cross was born May 7, 1924, in Chicago, 
Illinois. In September, 1940, he entered Oberlin College, 
Ohio and attended school until he was drafted into the 
Armed Forces in March, 1943. He served approximately 
three years as a rifleman with the 11th Air Borne Infantry 
Division, 511th Regiment, including 18} months in the 
Southwest Pacific Area. After an honorable discharge in 
November, 1945, he entered the University of Texas where 
he held a teaching fellowship. After he received a B.S. de- 
gree in physics in June, 1948, he accepted the position of a 
Junior Physical Engineer with the Magnolia Petroleum 
Company. In April, 1949, he resigned this position to re- 
enter Texas University for graduate work in geophysics. 
During the past year he has been working on his thesis as a 
Research Scientist Class I under a Navy research contract 
(N6onr-266, T. O. VIII). This project is also being sup- 
ported by yearly research grants from the Humble Oil and 
Refining Company and the Shell Oil Company. He ac- J. Harvey Cross 
cepted a position with the Humble Oil and Refining Com- 
pany as a Junior Exploration Geophysicist after having graduated with a M.A. degree in June, 1951. 

Mr. Cross is a member of the Society of Exploration Geophysics, American Geophysical Union, 
and Sigma Pi Sigma. 
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Fritz GAssMANN received his diploma in Mathematics 
and Physics in 1923, and his Ph.D. in Mathematics in 1925, 
both from the Swiss Federal Institute of Technology, Zurich 
(Switzerland). In the period 1925-42, he was at first as- 
sistant at the Swiss Federal Earthquake Service, then 
teacher and rector of the Kantonsschule Aarau and instruc- 
tor for Geophysics at the Swiss Federal Institute of Tech- 
nology. Further training in Geophysics led him to Jena, 
Naples, and Freiburg. Since 1942 he has been professor 
of Geophysics and director of the Institute of Geophysics 
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at the Swiss Federal Institute of Technology, Zurich. 


CLARK GOODMAN 








Fritz GASSMANN 


CLARK GOODMAN was born September 9, 1909, in Mem- 
phis, Tennessee. He received the B.S. in Chemical Engi- 
neering at California Institute of Technology in 1932, and 
the Ph.D. in Physics at Massachusetts Institute of Tech- 
nology in 1940. From 1932 to 1933 he served as research 
chemist for Union Oil Company; from 1933 to 1936 in the 
same capacity for the Kelco Company; from 1940 to 1942 
as research associate, 1942 to 1947 as assistant professor, 
and from 1947 to the present time as associate professor at 
Massachusetts Institute of Technology. He was consultant 
for the Manhattan District from 1943 to 1946; Technical 
Aide, Division 17, NDRC in 1945 and 1946; and Senior 
Physicist, Clinton Laboratories, Oak Ridge, Tennessee, in 
1946. 

Dr. Goodman is a Fellow of the American Physical So- 
ciety and the Geological Society of America. 








W. HEISKANEN 


CONTRIBUTORS 


W. HEISKANEN was born in Kangaslampi, Finland in 
1895. He received the degree of Ph.D. from the University 
of Helsinki in 1924. From 1921 to 1928 he served as geod- 
esist with the Finnish Geodetic Institute. He has been 
Assistant Professor of Geodesy at the University of Hel- 
sinki since 1926, and was Professor of Geodesy at Finland’s 
Institute of Technology from 1928 to 1949. He has served 
as Director of the Isostatic Institute of the International 
Association of Geodesy since 1949, and as Reporter on 
Isostasy of the Association since 1930. 

Dr. Heiskanen was a Member of the Finnish Parlia- 
ment from 1933 to 1936. He has been a member of the Fin- 
nish Academy of Sciences since 1928; member of the Board 
of the Finnish Cultural Foundation since 1937; honorary 
member of the Dutch Geographical Society since 1938; 
member of the American Geophysical Union and American 
Geographical Society since 1946; member of the Norwegian 
Academy of Sciences of Bologna since July 1950; and mem- 
ber of the Academia Nazionale dei Lincei in 1951. He has 


served as President of the Finnish Geographical Society, Finnish Astronomical Society, and Finnish 
Society of Christian Culture. He was Research Associate with the California Institute of Technology 
in 1948. Among the honors he has received are: Commandeur, Finnish White Rose, Estonian White 


Star, and Hungarian St. Stephan. 


Dr. Heiskanen has published five textbooks and about seventy-five scientific investigations on 


geodesy, gravimetry, and astronomy. 


DarRRELL S. Hucues received his A.B. and M.S. de- 
grees from the University of Kentucky in 1926 and 1928 
and his Ph.D. degree from the University of Chicago in 
1931. During the years 1931-1933 he held a National Re- 
search Fellowship at the California Institute of Tech- 
nology. After one year of teaching and one year on produc- 
tion research with the Kettleman North Dome association, 
he joined the Shell Oil Company as chief of a gravity party. 
During the next ten years he worked on the interpretation 
of gravity, seismic, and electric data, and the design and de- 
velopment of instruments. In September 1945, he joined the 
staff of the Physics Department of the University of Texas 
as Professor of Physics. Since that time he has been en- 
gaged principally in research on physical properties of rocks 
and minerals at high pressure and temperature. 
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CONTRIBUTORS 


HaArROLp R. Prescorrt received his B.S. degree in Elec- 
trical Engineering from the University of Kansas in 1927. 

He has been engaged in geophysical work since 1927, 
when he was employed by the Marland Oil Company in 
connection with gravity pendulum development. Employ- 
ment continued with Continental Oil Company, following 
its merger with the Marland Oil Company, in the capacity 
of Development Geophysicist. In 1949 he became Assistant 
Manager of the Geophysical Department of Continental 
Oil Company and in June 1951 was promoted to Manager 
of the Geophysical Section. 

He is a member of the Society of Exploration Geo- 
physicists, the S.E.G. Houston Section, the American 
Association of Petroleum Geologists, and the Houston 
Geological Society. 


managing director. 





D. J. SALT 


E. J. SrULKEN was born and reared in Texas. He stud- 
ied mathematics and physics at The University of Texas, 
receiving there his B.A. degree in 1934 and his M.A. degree 
in 1935. From 1935 to 1940 he was employed by Geophysi- 
cal Service, Inc. in the capacity of computer. During the 
following two years he served at Allen Military Academy, 
Bryan, Texas, as Instructor of Mathematics and Physics. 
Since 1942 he has again been associated with Geophysical 
Service Inc. serving successively as Chief Computer, Seis- 
mologist, and Chief Seismologist. Here he has worked 
chiefly on seismic computing methods and on special prob- 
lems of interpretation. 

Mr. Stulken is a member of the Society of Exploration 
Geophysics, The Mathematical Association of America, 
The Dallas Geological Society, and the Dallas Geophysical 
Society. 








HAROLD R. PRESCOTT 


D. J. Satt was graduated from the University of 
Toronto in 1948 in Engineering Physics. The summer of 
1948 he spent on the Seward Glacier in Alaska and the 
Yukon, conducting seismic operations. In November 1948 
he joined Mining Geophysics Corporation interpreting 
aeromagnetic surveys and conducting field work. In March 
1951, he organized Geo-Explorers Ltd., of which he is now 


Mr. Salt is a member of the Association of Professional 
Engineers of Ontario, the Arctic Institute of North America, 
and the Society of Exploration Geophysicists. 
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ACTIVE 


Charles Harry Carlisle (D. H. Gardner, H. G. Patrick, M. M. Slotnick) | 
Charles Burton Cox (Vincent Miller, R. M. Bradley, C. L. Brownlow) 

Tobias Flatow (D. H. Gardner, H. G. Patrick, M. M. Slotnick) 

Theodore Burroughs Fryer, Jr. (H. C. Bemis, R. C. Dunlap, G. J. Hulsewe) 

Robert Alvey Hall (J. C. Waterman, H. R. Thornburgh, Robert Dyk) 

Robert Linders Heineck (R. W. Mossman, H. M. Thralls, H. H. Andrews) 

Weikko Aleksanteri Heiskanen (Sigmund Hammer) 

Ralph C. Howard (John E. McGee, J. T. Murrell, W. B. Lee, Jr.) 

Fred Thomas Kafka (Leo J. Peters, L. W. Gardner, C. H. Dresbach) 

John Frederick McFarlane (Earle W. Johnson, T. O. Hall, Jack Martin) 

Frederick Allan McKinnon (A. D. Dunlap, H. W. Stoneman, Tom D. Mayes) 
Melvin John Rieger (Roy L. Lay, Vernon E. Briard, R. N. Buryere) 

Benjamin Boags Robinson, Jr. (C. R. Lowe, David Bell Campbell, Ralph S. Jackson) 
Ralph Thomas Smith (J. E. Lee, Jr., V. G. Feather, J. C. Palmer) 

Paul Harris Steed (E. V. McCollum, T. J. Thomas, Craig Ferris) 

Donald Desmond Utterback (Paul S. Lewis, John D. Marr, Karl W. Abel) 

Cornelis Jan Velzeboer (Ch. Dresbach, A. van Weelden, P. Mrosovsky) 











ASSOCIATE 


Robert Berry Arnold (R. G. DeGood, E. J. Handley, L. D. Ervin) 

Adrian Anthony Becker (John M. Crawford, Roland F. Hughes, Bird G. Swan) 
Owen Burtis Bennett (A. I. Innes, Andrew Gilmour, Albert Taylor) 

Herry Stinson Birdseye (Leland C. Adams, Joseph Fishback, K. C. Thompson) 
Louis David Blom (A. Klaveness, C. C. Zimmerman, H. W. McDonnold) 
Leonard Keith Buller (Sidon Harris, R. H. Dana, Elmer F. Blake) 

fulian Martin Busby (Erdelyi Fazekas, David Reed, William B. Pine) 

Troy Norman Crook (W. B. Duty, R. A. Lowery, Robert S. Duty, Jr.) 

Julius Perry Edge (N. K. Moody, R. G. Harbach, R. F. Aldridge) 

Betty Gaines (W. E. Jasper, O. B. Jackson, L. F. Fischer) 

Marvin Grace (A. B. Hamil, Gordon Mace, Sam Pena) 

Robert James Graebner (R. C. Dunlap, Jr., Paul E. Swenson, T. P. Ellsworth) 
John E. Hall, Jr. (R. F. Bennett, R. D. Roberts, C. B. Smith) 

James Kenneth Harrison (A. Klaveness, C. C. Zimmerman, H. W. McDonnold) 
Donald Edwin Heilner (Kenneth E. Burg, Robert C. Dunlap, Jr., Homer L. Heggy) 
Ray Joseph Hope (T. O. Hall, Chester Sappington, Ronald Wright) 

James Richard Hunt (H. E. Itten, C. J. Long, P. H. Garrison) 

George Roy Jones (E. J. Handley, R. G. DeGood, Orville Strozier) 

Patrick Francis Kelly (S. E. Giulio, R. F. Weichert, R. J. Watson) 

Eugene Louis Kessler (Karl Dyk, M. B. Widess, H. E. Itten) 
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j Thomas Carroll King (R. M. Dreyer, T. G. Arnold, Craig Ferris) 

Joseph Cowden Knight (Thomas O. Hall, John P. Paschall, M. R. Henderson) 
Morris Andrews Long (R. T. Gallagher, S. A. Spencer, W. B. Agocs) 

Vincent D. Mackle (L. F. Athy, H. R. Prescott, W. C. Kimball) 

Richard George Martin (Carl W. Blakey, Thomas A. Manhart, B. W. Beebe) 
Robert Barnard O’Connor, Jr. (F. Goldstone, C. H. Gregory, W. Hafner) 
Budd Orlich (Thomas A. Manhart, B. Warren Beebe, Carl W. Blakey) 
Everett Horace Parker, Jr. (J. G. Harrell, A. E. McKay, Sidon Harris) 

Sam Walker Patterson (R. G. Reiser, C. M. Wert, G. R. Belt) 

Thomas Bruce Portwood, Jr. (R. F. Weichert, C. N. Hurry, W. W. Clark) 
Louis Victor Provensal (Roy L. Lay, A. A. Hunzicker, B. H. Treybig, Jr.) 
Aubrey Henry Rabensburg (L. L. Nettleton, C. B. Claypool, J. A. Culbertson) 
Robert William Riedel (W. A. Longacre, L. O. Bacon, B. H.,.Boyum) 

Harold Herbert Sawyer (Judson Mead, Andrew H. McNair, Charles Deiss) 
Eugene Theodore Schnieder (R. W. Mossman, H. M. Thralls, H. H. Andrews) 
John David Sistrunk, Jr. (U. A. Rowe, Chester Sappington, E. W. Johnson) 
Melvin Mathein Smith (Harold K. Herr, Deane J. Wolff, H. L. Mendenhall) 
Mohammed Wafick (J. R. Cornett, Sam D. Rogers) 

Warren Alton Wilson (J. G. Harrell, A. E. McKay, Sidon Harris) 
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Jacobus Willem De Bruyn (A. van Weelden, B. Baars, V. P. Ulrich) 

Richard Longley Fentem (A. I. Innes, Andrew Gilmour, Albert Taylor) 

Kasmier Dominic Kubik (A. B. Hamil, G. M. Mace, Jr., B. K. Johnson) 

Wallace Edward Pennington (A. I. Innes, Andrew Gilmour, A. W. Taylor) 

Robert Frank Pomeroy (Leo J. Peters, Sigmund Hammer, Thomas A. Elkins) 

John Joseph Schneider, Jr. (Chester P. Sappington, Thomas O. Hall, Roy F. Bennett) 
Merrill Smith (Roy L. Lay, C. R. Wallace, A. A. Hunzicker) 





TRANSFERS TO ASSOCIATE 


John Mike Donna (James B. Macelwane, Kenneth E. Burg, Fred J. Agnich) 
Harold Lloyd Krivoy (H. V. W. Donohoo, Henry R. Joesting, Philip Donnerstag) 













































SOCIETY ROUND TABLE 


NOMINEES FOR 1952-53 EXECUTIVE COMMITTEE 


The following brief biographical sketches are presented to acquaint members of the society with 
nominees appearing on the official ballot which will be mailed to all voting members at a later date. 


FOR PRESIDENT 











Curtis H. JOHNSON GERALD H. WESTBY 


Curtis H. JOHNSON was graduated from the University of California at Los Angeles in 1932 
with an A.B. and Honors in Physics, and took postgraduate work in acoustics also at U.C.L.A. From 
1933 to 1938 he was geophysicist and later chief geophysicist with Rieber Laboratory. From 1939 
to the present date he has been employed by General Petroleum Corporation where he currently 
holds the position of assistant chief geophysicist. 

Mr. Johnson is a member of the Society of Exploration Geophysicists, The American Association 
of Petroleum Geologists, and the American Geophysical Union. He has served the S.E.G. as vice chair- 
man of the arrangements committee for the national convention in Los Angeles in 1947, chairman 
of the program and arrangements committee for the 1947 Pacific Coast Regional Meeting and chair- 
man of the Pacific Coast program committee for the 1948 annual meeting. In 1948 he participated in 
the organization of the Pacific Coast Section of the Society of Exploration Geophysicists and was 
elected its first president. He is at present Vice President of the Society of Exploration Geophysicists, 
and is General Chairman of the program and arrangements committee for the Twenty-Second Annual 
Meeting which will be held March 23-26, 1952, in Los Angeles. 

His other recent activities on behalf of the society include membership on the Standing Com- 
mittee on Honors and Awards, 1949-50 and 1950-51, and as member-at-large on the Standing Com- 
mittee on Program and Arrangements for the 1951 annual meeting. 

GERALD H. WEstBy was born in Stoughton, Wisconsin, of Norwegian-English parents, on 
August 4, 1898. He received his elementary education in the LaCrosse, Wisconsin, public schools and 
was graduated from high school there in 1916. He attended the University of Chicago and, because of 
the influence of Professor R. D. Salisbury, majored in geology. He was graduated from the University 
of Chicago with his class in 1920, since credit was given for army service in the Coast Artillery Officers 
Training School, Fort Monroe, Virginia, where he received a commission as Second Lieutenant. 

In the summer of 1920, Mr. Westby taught physiography and geography at the University of 
Montana, and early in the fall left to accept a position with S. Pearson & Son, London, England, to 
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do geological work on various assignments throughout the world. From late 1920 to early 1924, he 
was resident geologist for S. Pearson & Son in Algeria. After completing three months of research 
work at Stanford University, Mr. Westby accepted a position as field geologist with the Empire 
Oil and Refining Company, with field headquarters at Denver, Colorado. In the fall of 1926 and in 
early 1927, while still employed by Empire, he attended the first class in geophysics to be offered at 
the Colorado School of Mines. In late 1927 he was made chief geophysicist of the Empire Oil and 
Refining Company and maintained that position until September 1933. 

Since 1933 Mr. Westby has been associated with Seismograph Service Corporation, Tulsa, 
Oklahoma, and has been president of the company since 1935. He is also president of Seismograph 
Service Corporation of Canada, operating in that country, Seismograph Service Corporation of Dela- 
ware, operating in Mexico and Venezuela, and chairman of the board of directors of Seismograph 
Service Limited, an English subsidiary operating in the Sterling Bloc areas. 

Numerous articles by Mr. Westby have been published in the trade magazines and several in 
the early issues of the Society publications. 

Mr. Westby is a member of the Society of Exploration Geophysicists, the Geophysical Society 
of Tulsa, The American Association of Petroleum Geologists, the American Geophysical Union, the 
American Association for the Advancement of Science, and the American Institute of Mining and 
Metallurgical Engineers. He was mid-continent representative on the Founding Committee of the 
Society of Exploration Geophysicists and is, therefore, a charter member. He served the Society as 
secretary-treasurer in 1935-1936 and as vice president three times (1931-1932, 1932-1933, and 


1936-1937). 
FOR VICE-PRESIDENT 





O. C. CLIFFORD, JR. Roy L. Lay 


Roy L. Lay received his B.S. degree in Electrical Engineering from Rice Institute in 1928. Upon 
graduation from Rice he was employed by the Geophysical Research Corporation as an observer and 
was later promoted to computer and party chief. 

Mr. Lay was employed by The Texas Company as party chief in April, 1932 and during the next 
several years initiated this company’s reflection work in Texas, Louisiana, California, and the Rocky 
Mountain area. He was promoted to field supervisor in 1937, and during the next year made a super- 
visory trip to Colombia and Venezuela, and was in charge of the seismic, gravity and magnetometer 
work in Egypt. In July, 1943 he was promoted to assistant chief of the geophysical division and in 
September, 1945 was advanced to the position of assistant manager. Gn April 1, 1950 he became 
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manager of the geophysical division, the position which he now holds, having jurisdiction over The 
Texas Company’s geophysical operations in the United States and Canada. 

He is a member of the American Association of Petroleum Geologists, Houston Geological So- 
ciety, and has been a member of the Society of Exploration Geophysicists since 1935. He was Second 
Vice-President of the S. E. G., Houston Section, in 1947-1948, President for the 1948-1949 term, 
and served on the executive committee of that section from 1949 to 1951. 

O. C. CLIFFORD, JR., received his first two years collegiate work at the University of Wyoming. 
He was graduated from the University of Chicago in 1924 with a B.S. in Geology, remaining at that 
institution the following year engaged in graduate study. He was employed in 1925 by Humphreys 
Corporation on one of the earliest torsion balance parties. Between 1926 and 1933, he held positions 
in gravity, magnetics, seismics, and subsurface geology for a variety of companies with experience in 
all of the oil producing states of the Mid-Continent. In 1934 he was employed by Continental Oil 
Company as control geophysicist with headquarters in Ponca City, Oklahoma. In 1936, he joined 
The Atlantic Refining Company as district geophysicist for the Gulf Coast area of Texas. In 1943, 
he was transferred to the main office in Dallas as chief geophysicist, the position he currently holds. 

Mr. Clifford is a member of the Society of Exploration Geophysicists, The American Association 
of Petroleum Geologists, and the American Geophysical Union. He has served the S. E. G. as a 
District Representative, and as a member of the editorial board of the Case Histories volume. He is a 
past president of the Dallas Geophysical Society. He has written for the A.A.P.G. bulletin and is 
represented as a co-author in Geophysical Case Histories, Volume One. 


FOR SECRETARY-TREASURER 





Cart L. BRYAN James L. Morris 


Cart L. BRYAN was born in Guadalupe County, Territory of New Mexico, on August 10, 1908. 
Three years of undergraduate work in electrical engineering at the University of Texas were followed 
by employment in June 1928 in the geophysical department of Gulf Production Company at Houston. 
He was transferred to Shreveport, Louisiana in 1930 as assistant torsion balance supervisor. From 
1932 to 1942 he served various divisions and affiliates of the Gulf Oil Corporation in numerous ca- 
pacities on torsion balance, pendulum, gravimeter, and seismograph field crews on domestic and 
foreign assignments. 

He entered the U. S. Navy in 1942, served three years aboard aircraft carriers in the Pacific 
Ocean Area, during which time he attained the rank of Lieutenant Commander, and was transferred 

















— IT 





SOCIETY ROUND TABLE 731 


to inactive duty in January 1946, at which time he rejoined the Gulf organization. Shortly thereafter 
he was assigned his present position as chief geophysicist for Gulf Refining Company at Shreveport. 
In 1950 he received the B.S. degree in geology from Centenary College of Louisiana. 

Mr. Bryan has been a member of the Society of Exploration Geophysicists since 1939. He is also 
a member of The American Association of Petroleum Geologists, American Geophysical Union, 
Ark-La-Tex Geophysical Society, Shreveport Geological Society, and Sigma Gamma Epsilon, honor- 
ary geological fraternity. He has served successively as Secretary-Treasurer and President of the 
Ark-La-Tex Geophysical Society. 

James L. Morris was graduated from West Virginia University in 1936 with an A.B. degree in 
Geology and received the degree of Geological Engineer (Geophysical Option) in 1938 from the Colo- 
rado School of Mines. He was employed upon graduation by The Pure Oil Co., beginning as a mag- 
netometer operator. Subsequent to further geological and geophysical experience in various divisions 
of the company he served as district geologist, then as a geophysical supervisor of seismic and gravity 
operations. In 1948 he was promoted to division geologist of the Texas producing division of The 
Pure Oil Company and since 1949 has served as exploration manager in that division. 

He is a member of the American Association of Petroleum Geologists, Ft. Worth Geological 
Society (Vice President in 1944), the Society of Exploration Geophysicists and the Ft. Worth Geo- 
physical Society (District Representative 1949-50 and Past President). 


SOCIETY OF EXPLORATION GEOPHYSICISTS, HOUSTON SECTION 
SECOND AWARD OF HONORARY MEMBERSHIP 


CITATION 
T. I. Harxins 


The Award Committee of the Houston Section, S.E.G. has asked me to give a resume of E. E. 
Rosaire’s activities in the geophysical field and it was decided that it would be worth while mention- 
ing some of his background before entering the field of geophysics as well as those subsequent to it. 

It gives me a great deal of pleasure to be able to do this because I, along with numerous others, 
have the fee ling that Gene did more to help geophysics to attain its present position than any other 
one individual. . 

He received his M.S. degree from Chicago in 1921 and his Ph.D. from the same school in 1926. 
He was in the Development Branch of the Western Electric from 1922 to 1925 and was employed 
by the Geophysical Research Corporation and sent to Texas in the spring of 1925 to organize two 
field parties contracted to the Gulf Production Company (now Gulf Oil Corporation). While he was in 
charge of the second party the Moss Bluff salt dome was found. This was the first American seis- 
mograph crew to find and locate a salt dome. Three weeks later this same crew discovered the Port 
Barre dome in south Louisiana. 

Rosaire was transferred to Oklahoma in 1926 to direct Amerada Petroleum Corporation’s re- 
flection work. He came back to the Gulf Coast in 1927 where he had charge of all Geophysical Re- 
search Corporation’s contract work. He remained with G.R.C. until 1932 and during this time re- 
flection dip shooting was developed in the Gulf Coast. In 1932 he helped to organize the Independent 
Exploration Company which consisted of a number of ex-G.R.C. men. As head of the Independent 
Exploration Company he originated and developed the idea of using traverse misclosures to detail 
faults associated with Gulf Coast structures. He resigned as president of the Independent Exploration 
Company in 1937 to enter the soil analysis and electrical prospecting field. He originated the idea of 
analyzing well cuttings to predict the nearness of the bit to the oil sands. He presently lives in Dallas 
and is a partner of the Geochemical Surveys, Inc. 

In 1930 the Society of Petroleum Geophysicists was organized in Houston, Texas, and Gene was 
one of the charter members and helped to organize this association. The late Donald Barton was 
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elected president, E. E. Rosaire, vice-president, and John Weinzerl, secretary-treasurer. Gene was 
elected president in 1934. He helped to organize and write the charters and by-laws of the Houston 
section of the Society. He has always been very active and interested in all phases of the Society’s 
work, 

Gene was an Air Force Lieutenant during the First World War. He was elected to the Order of 
Quiet Birdmen in 1927. This is an association of ex-military fliers who have achieved some special 
distinction in the war other than in the line of duty. His particular achievement was that of being the 
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leading German ace in the American army. He w-: ckc:i five U. S. planes in one day, which was his 
biggest day. 

He originated four handed chess and has writien a book on rules and how to play the game. He 
did considerable work on tax research in the City of Houston and helped to reorganize the Tax Re- 
search Association. 

The above facts are only the highlights of Gene’s career and there are numerous others that 
should be mentioned but it was the purpose of this citation to be brief and touch only on his out- 
standing achievements. Again I want to say that this has been a pleasure to me, and I want to thank 
the Committee for favoring me with the opportunity of preparing this résumé. 


EASTERN REGIONAL MEETING 


The Executive Committee of the Society of Exploration Geophysicists has authorized an Eastern 
Regional Meeting which will be held in Pittsburgh, Pennsylvania on October 25, 1951. 

The one-day program is planned to be of general interest to geophysicists and geologists and will 
not be limited to local subjects. The geophysicists will be guests of the Pittsburgh Geological Society 
for an evening meeting at which Dr. W. C. Krumbein will speak. 

The meetings will be held in the Mellon Institute Auditorium and registration will be in the ad- 
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jacent foyer. Blocks of rooms for this meeting have been reserved by two hotels; the Webster Hall, 
across the street from Mellon Institute, will be Headquarters Hotel and has a few single rooms from 
$5.50 to $7.00 and double rooms from $8:00 to $12.00. The Schenley Hotel, one block from Mellon 
Institute, has double rooms from $8.00 to $12.00. When requesting reservations be sure to notify the 
hotel that you will attend the SEG Eastern Regional Meeting. : 

Everyone registered at the meeting will be invited to tour the Mellon Institute of Industrial 
Research and the Gulf Research & Development Company laboratories on the following day, Friday, 
October 26th. 

R. B. Ross has been named General Chairman of the meeting and with Dr. Frank Press is ar- 
ranging the program. Local arrangements are being made by L. F. Melchior assisted by W. B. Robin- 
son and David A. Davis. 


ANNUAL PACIFIC COAST MEETING 
November 1 and 2, 1951 


The Pacific Coast Section, Society of Exploration Geophysicists, will hold its annual fall meeting 
on November 1 and 2 at the Ambassador Hotel in Los Angeles. The meeting will again be held in 
conjunction with the annual sessions of the Local AAPG and SEPM groups and coordinated programs 
will be presented. Flint Agee, of United Geophysical Co., Pasadena, vice president of the southern 
division of the Pacific Coast Section, will be in charge of program and arrangements for the S.E.G. 

More than 200 exploration scientists on the West Coast are members of the Pacific Coast Section. 
Officers of the society are: John Sloat, Union Oil Co. of California, Bakersfield, vice president, North- 
ern District; Flint Agee, vice president, Southern District; O. K. Fuller, Jr., Humble Oil & Refining 
Co., Los Angeles, Secretary-Treasurer; John D. Hale, Seaboard Oil Co., Bakersfield, district repre- 
sentative to the S.E.G.; and V. E. Prestine, Western Geophysical Co., Los Angeles, district representa- 
tive to the S.E.G. 


ANNUAL MIDWESTERN REGIONAL MEETING 
November 19 and 20, 1951 


The Baker Hotel, Dallas, Texas, will be the scene of the annual joint meeting sponsored by the 
Dallas Geophysical Society, Fort Worth Geophysical Society, Permian Basin Geophysical Society, 
Ark-La-Tex Geophysical Society, and Geophysical Society of Tulsa. The Tulsa Section will be in 
charge of general arrangements through a committee composed of E. M. McNatt, J. P. Garner, 
K. M. Lawrence, W. H. Courtier, Frank Searcy, and Karl Dyk, Chairman. Most of the papers pre- 
sented during the two-day meeting will concern geophysical problems encountered in the midwestern 
oil producing region of the United States, including a featured address by Dr. B. B. Weatherby. 
More than one thousand are expected to attend. Detailed announcements will be mailed to all 
members of the society well in advance of the meeting, along with hotel reservation forms. Those who 
await receipt of the official announcement will be assured of obtaining better accommodations, 
since the hotels have reserved space for delegates who apply on the official form. 


TWENTY-SECOND ANNUAL MEETING 
March 23-27, 1952 


The next Annual Meeting of the Society of Exploration Geophysicists will be held March 23-27, 
1952 in the Biltmore Hotel, Los Angeles, California. The society will again meet jointly with The 
American Association of Petroleum Geologists and Society of Economic Paleontologists and Mineralo- 
gists. Curtis H. Johnson, Vice President, is general chairman of program and arrangements for the 
S.E.G., and Furman A. Grimm is vice chairman in charge of local arrangements. 

The official announcement of this meeting will be mailed in the near future to all members of 
the society, along with the official hotel reservation form. Members are urged to await receipt of the 
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officia] announcement before requesting accommodations, inasmuch as the hotels will reserve their 
more suitable rooms for those whose requests are made on the official form. 


ANNUAL GULF COAST MEETING 
May 29 and 30, 1952 


Tentative dates have been set by the Houston Section for a two-day meeting on May 29 and 30, 
1952 at the Rice Hotel in Houston, Texas. Although the actual plans of the meeting have not been 
formulated, this announcement is made for the benefit of those companies and individuals who wish 
to schedule the attendance of their personnel in advance. Development of program plans and ar- 
rangements will be announced in future issues of Geophysics, and an official notice will be mailed to 
all members of the society well in advance of the meeting. 

Eugene W. Frowe, of Robert H. Ray Company, vice president of the Houston Section, is in 
charge of program and arrangements. 


DEMAND INCREASING FOR OUT OF PRINT ISSUES 


Recent expansion of curricula to include courses in exploration geophysics has created a further 
demand for complete files of Geophysics in universities throughout the world. Almost daily the business 
office must report to new subscribers that certain back numbers of the Journal are not available. 
Members of the society and subscribers to Geophysics are urged to report to the business office any 
of the following numbers for which they have no further use, in the interest of making available this 
essential literature on exploration geophysics. 


1938 Volume III Numbers 1, 3 and 4 
1940 Volume V Number 1 

1941 Volume VI Numbers 1 and 4 
1942 Volume VII _ Number 1 

1943 Volume VIII Numbers 1, 2, 3 and 4 
1944 Volume IX Numbers 1, 2, 3 and 4 
1945 Volume X Numbers 1 and 2 
1947 Volume XII Number 1 


PERSONAL ITEMS 


K. A. RoBertson, formerly manager of the Calgary office of United Geophysical Company, is 
now president of Precision Exploration Company with offices at 409 Palace Building, Tulsa, Okla- 
homa. 


Joun W. Daty, West Texas supervisor of General Geophysical Company, spoke at a meeting 
of the Canadian Society of Exploration Geophysicists August 24 on the subject, “How Seismograph 
Interpretation is Related to Changes in the Sedimentary Section.” The meeting was held in the 
Palister Hotel, Calgary. 


LoreNz SHOCK announces the opening of the Midland Geophysical Company with offices at 233 
Capitol Building, Midland, Texas. 


Jacques ScCHOEFFLER, field geologist with Societe Cherifienne des Petroles, Rabat, Morocco, 
is in the United States under the sponsorship of the Economic Cooperation Administration. While 
here he will obtain training in geophysics at the University of Texas and Colorado School of Mines, 
and at the same time he is visiting various contracting and producing companies for what the ECA 
calls “‘in-plant”’ study of oil exploration in the United States. Upon returning to Rabat Mr. Schoeffler 
will become chief geophysicist of his company. He has visited the business office of the society and 
has applied for membership. 
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Rosert Dyk, formerly senior geophysicist for General Petroleum Corporation, has joined Geo- 
physical Service Inc., as manager of the Rocky Mountain division with headquarters in Denver, 
Colorado. Prior to this change he was located in Bakersfield, California, and was president of the 
Pacific Coast Section of S.E.G. 


R. G. SOHLBERG, formerly in charge of the geophysical division of Richmond Exploration Co. 
(Standard Oil Co. of California) at Maracaibo, Venezuela, has resigned his connection with that 
organization and has accepted the position of manager for United Geophysical Company of Canada, 
with offices at Calgary, Alberta. 


H. K. Herr, formerly with Kerr-McGee Oil Industries, is now seismic party chief for the Ray- 
flex Exploration Company at Miles City, Montana. 


W. J. BusTEED, former Secretary-Treasurer of the Colorado School of Mines Society of Student 
Geophysicists, is now with Geophysical Service Inc., at 6000 Lemmon Avenue, Dallas 9, Texas. 


F. REx WALLING, Jr. has resigned from Stanolind Oil & Gas Co. and is now connected with the 
Grant Oil Tool Co., and MacClatchie Manufacturing Co., at Odessa, Texas. 


STANFORD L. Myers has been appointed equipment and safety supervisor of Southern Geo- 
physical Company with offices in Fort Worth and Midland, Texas. 


Rosert L. GEYER has resigned from the faculty of the School of Geology, Louisiana State Uni- 
versity, to accept the position of geophysicist with the Stanolind Oil & Gas Company. His address 
is 1815 North Cheyenne, Tulsa 6, Oklahoma. 


J. S. WELBOAN, JR. has been promoted to Assistant to Division Manager of the geophysical 
division, producing department of The Texas Company. Employed shortly after graduation from 
the University of Oklahoma in 1935, Mr. Welboan has served the Texas Company successively as 


computer, party chief, and supervisor. 


A. E. McKay, formerly assistant chief geophysicist for The Atlantic Refining Company, and 
from June 1950 to July 1951, executive vice president of Southern Geophysical Company, announces 
the organization of Continental Geophysical Company with offices at 6100 Camp Bowie Blvd., Fort 
Worth, Texas. 

E. F. K. Zarupski has resigned from Seismograph Service Ltd., London, to accept the position 
uf party chief with Century Geophysical Corporation in Canada, at 615 8th Avenue West, Calgary, 
Alberta. 

J. L. Keze ver has been named executive vice president in charge of field operations of Standard 
Seismograph Company, Oklahoma City, it was announced by the company president, C. C. Copy. 
Kezeler recently resigned from the position of chief geophysical engineer with Century Geophysical 
Corporation to accept the new position. 


S1pon Harris, president of Southern Geophysical Company, announces publication of the first 
issue of ‘Southern Record,” company house organ to promote good will among Southern employees 
and to keep widely separated crews informed on activities of fellow workers. 


C. T. RANkIN, party chief with Amerada Petroleum Corp., has been transferred to Glasgow, 
Montana. 


Jor L. Moore has resigned from Honolulu Oil Corporation to establish offices as consultant in 
Midland, Texas. 


CHARLES C. BATEs, consultant with A. H. Glenn & Associates, has moved from Washington, 







































736 SOCIETY ROUND TABLE 





D. C., to College Station, Texas, to complete work toward a Ph.D. degree in geological oceanography 
at A. & M. College of Texas. 


RosBert K. Youne, party chief for Continental Oil Co., is now located in Duncan, Oklahoma. 


Donan M. Davis has opened consulting geological offices at 8223 Moraga Drive, West Bel-Air, 
Los Angeles 49, California. Mr. Davis was formerly with the Union Oil Company of California and 
Richfield Oil Corporation. Prior to this he was with United Geophysical Company in the capacity 
of party chief. He will specialize in problems relating to the Sacramento Valley of California. 


Puri P. Gasy has resigned from the California Standard Company, and is now engaged as a 
partner in Delta Exploration Company and as manager of the Canadian branch office located at 
536 8th Avenue West, Calgary, Alberta. 


W. Day WA pRE?P has taken the position of senior map draftsman in the exploration department 
of Sunray Oil Corporation at 1509 First National Building, Tulsa, Okla. 


R. H. Dana has been promoted from vice president to executive vice president of Southern 
Geophysical Company. W. D. Barrp has moved from treasurer and supervisor to vice president, 
and R. E. Davis from supervisor to assistant vice president of the company. 


THE NATIONAL RESEARCH COUNCIL announces that its Committee on Geophysics, Advisory to 
the Office of Naval Research will meet 7 November, 1951. The Committee will advise the Office of 
Naval Research on proposals received by it in the field of geophysics. Those wishing support in basic 
research in this field should submit five copies of their proposal and a letter of endorsement signed 
by the President or other duly authorized financial officer of their institution. These should be sent 
in each case to the Commanding Officer of the ONR area branch office for that district. In addition, 
at least one advance copy should be sent to the Geophysics Branch, Office of Naval Research, Depart- 
ment of the Navy, Washington 25, D. C. To insure consideration at the forthcoming meeting of the 
Committee, such proposals should be received by ONR in Washington not later than 1 October, 1951. 


THE AIR FORCE ANNOUNCES that men in civil life, either veterans or non-veterans with a degree 
in civil engineering, geology, or photogrammetric engineering can now enter the Air Force directly 
as officers, with a topographic engineer commission. This announcement was made at First Air Force 
headquarters at Mitchel Air Force Base, New York, by Major Charles D. Morat, Jr., Director of 
Military Personnel Procurement. 

Major Morat said that experience to qualify a man for this commission must be in these fields: 
engineering, surveying, cartography, topographic drafting, geodetic surveying, or other related fields. 
Experience in lithographic and reproduction work is also desirable. More specifically, an applicant 
must have technica] knowledge of survey work, such as determining exact locations and measure- 
ments of points, elevation lines, areas and contours on the earth’s surface; or, should be experienced 
in photogrammetric and cartographic draftsmanship; or, photographic laboratory methods, photo- 
transferring and lithographic processes, or the operation and maintenance of lithographic equipment. 

Appointments will be made in grades from first lieutenant through lieutenant colonel, depending 
on age. Inquiries should be addressed to the Department of Military Personnel Procurement, Head- 
quarters First Air Force, Mitchel Air Force Base, New York, or phone Garden City 3-4000, extension 
3202. 
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431, Midland, Tex. 


Second Vice President: C. Newton Page, Honolulu 
Oil Corp., Box 1391, Midland, Tex. 

Secretary: Lou G. Cornish, Seasonee Service 
Corp., Box 1825, Midland, 

Treasurer: G. A, Grimm, Tide Water Associated 
Oil Co., 204 Ridglea, Midland, ‘ 
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TEXAS | 


SEISMOGRAPH COMPANY 





Rheroagh PETROLEUM 
practical GEOLO GY 


up-to-the-minute 
Kenneth K. Landes 


© Complete: covers all three divisions of 
petroleum geology—geological occur- 
rence, geographical distribution, tech- 
niques for locating oil and gas deposits. 


®@ Up-to-the-minute: includes data on 
newly discovered fields, and on oil and 
gas papers published as late as July, 1950. 
@ Practical: emphasizes how oil or gas 
is trapped, and how these traps may be 
discovered. 

@ Well-illustrated: contains 241 illustra- 
tions, of which 34 are full-page size and 
have full-page legends. 


1951 660 pages $10.00 





For sale by: 
Society of Exploration Gicaiantidies 
Box 7248, Tulsa 18, Oklahoma 














ANNOTATED BIBLIOGRAPHIES 
OF ECONOMIC GEOLOGY 


Available—Vols. I-XVII (1938-1944) ; 
Vol. XIX (1946); Volume XX (1947) ; 
Volume XXI (1948); Volume XXII 
(1949) 

XVIII (1945) in preparation. 

The intervention of World War II pre- 
vented receipt of foreign journals and 
disrupted the work of the Committee 
and Abstractors. The numbers have, 
therefore, been greatly delayed. 


General Index to Vols. I-X 


Price $5.00 per volume anywhere 
in the world 


ECONOMIC GEOLOGY JOURNAL 
INDEX TO VOLS. I-XL 
Index to Vols. I-XX (1906-1926) ..$3.00 
Index to Vols. XXI-XXX (1927- 
1935) 2 
Index to Vols. XXXI-XL (1936- 


Available from: 
Economic Geology Publishing Company 
100 Natural Resources Building 
Urbana, Illinois 











AEROMAGNETIC SURVEYS LIMITED 
WORLD-WIDE AIRBORNE MAGNETOMETER SURVEYS 


1450 O'Connor Drive 


Toronto, Canada 


Associated with 


Photographic Survey Corp. Ltd., Toronto, Canada © Geotechnical Corp., Dallas 
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With Kodagraph Autopositive Paper, you can make 
positive copies of your seismic records without the 
usual negative step. Your own direct-process or blue- 


In a hur r y print equipment and field chemicals will do the job. 


i 9 And the wide exposure latitude and room-light proc- 
for copies y essing of Kodagraph Autopositive Paper save you 
time and money—every copy is exactly right. A leaflet 
giving full details is yours for the asking. Eastman 
Kodak Company, Industrial Photographic Division, 
Rochester 4, N. Y. 


Kodagraph Autopositive Paper 
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CHECK LIST 


FOR GEOPHYSICAL SUPPLIES 








Electronic Components Loading Poles 
Geophysical Cable Test Equipment 
Field Supplies Hawthorne Bits 


Leece-Neville Alternators 


INDUSTRIAL <> ELECTRONIC SUPPLY CO. 











134 Leslie 1124 East Forth 

Dallas, Texas Tulsa Oklahoma 

Phone RA-4563 Phone 3-8122 
By William H. Twenhofel AKES a comprehensive presentation of the standard 
Professor Emeritus of Geology, methods of sampling for various types of sediments 
University of Wisconsin and rocks, gives outstanding methods for accurate analy- 
and Stanley A. Tyler ses, and shows miscellaneous forms of graphical repre- 
a vineveage | = Geology sentation of the characteristics of sediments and sedi- 
ee Te mentary rocks. Includes new and important treatments 
183 pages, 24 tables, illustrated with of coal and oil shales, and of igneous rocks as sources of 
charts, graphs, and halftones sediments. In all, the book shows in adequate detail how 
to analyze sediments with greater accuracy and to appraise 

$3.00 the results of the examination with more surety. 

For sale by 
SOCIETY OF EXPLORATION GEOPHYSICISTS 

Box 7248 Tulsa 18, Oklahoma 





| 


JOURNAL OF APPLIED PHYSICS 


This is a monthly journal designed particularly for those applying physics 





in industry and in other sciences. It publishes reviews of recent progress in 


applied physics, original research papers, news, and advertisements. 





Subscription Price Ue Ee. Foreign 
To members of American Institute of Physics ...... $10.00 $11.00 
IE A iidce Ni detisavktaradnnbiesdeenieel 12.00 13.00 

Single copies—$1.00 
Address 
AMERICAN INSTITUTE OF PHYSICS 
57 East 55 Street New York 22, N.Y. 
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ROTARY & SPUDDER 
SHOT HOLE & CORE 
DRILLING 


MF P0.BOX 325 PHONE 2-3452 
LUBBOCK, TEXAS 




















, 


AERO SERVICE CORPORATION 
and 


FROST GEOPHYSICAL CORPORATION 


announce the formation of 


FROST AIRBORNE SURVEY CORPORATION 
of Tulsa, Oklahoma 


For the conduct of aeromagnetic surveys under patent 


rights of Gulf Research and Development Corporation. 
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A 
DEPENDABLE 
GRAVITY 
METER 
FOR 
ACCURATE 
SURVEYS 





THE STABILITY and regularity of 
drift, the precision, and the 
ruggedness built into the North 
American Gravity Meter to with- 
stand field handling without 
taking a “set” when jolted or 
jarred, make it one of the most dependable meters available. 


These factors combine to give accurate and dependable readings the 
first time . . . saving the time and cost of remaking the survey. 


Where weight or size is important you will find the North American 
Gravity Meter to be extremely compact and light. One man can conveniently 


carry it on his back .. . it fits easily into a small boat or canoe .. . it can be 
read from its mounting in sedan, jeep, or even helicopter . . . it fits well 
in the space available in a diving bell. 

When you next order... or make ...a gravity meter survey, a North 
American Gravity Meter is your assurance of an accurate and economical 
survey. 
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| TIME PROVEN 
unitizen Fotlable 
TIC 24 SEISMOGRAPH 


ulvere- you ti nd- 


‘bath geographically and. geophysically from 
‘the glacier ‘regions of Alaska to. the 
“tropical jungles of New Guinea 





IN FHIS-RUGGED WATER TIGHT 
STAINLESS ‘STEEL CASE IS... 


| e 25 Trace Camera 

200 Ft. 6\In. Width Paper 

e 6 Vt. Governor Controlled Motor 
e Accurate Timing Fork System 

e 100.Cycle. Timing Check 

¢ Communication System 

e Automatic Trip Circuit 





x 
‘ 


Approxishately “40. sets of TIC. port- 
able ree CIS are in” eperatiqn 








_ res 
Instrument Co. 


3732 WESTHEIMER, HOUSTON, TEXAS 














“Some 15,000 employees 
in all branches of our 
organization are taking 
advantage of the 
Payroll Savings Plan...’ 


> 


SPENCER LOVE 
Chairman of the Board 
Burlington Mills 


“The Payroll Savings Plan is one of the important personnel services 
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offered to the employees of our company. At the present time some 15,000 


e 
sp 


maples in all branches of our organization are taking advantage of this 
endid plan for systematic savings. In times of national emergency this 


plan assists in stabilizing the economic life of the employee, the community 


and the nation.” 


Last call for the 1951 Defense Bond 
Campaign! 

hile the campaign was scheduled 
for six weeks, beginning Labor Day, 
the accounting period will include all 
Payroll Savings Plan bond purchases 
and enrollments in September and 
October. 

If you haven’t conducted a person- 
to-person canvass to put a Payroll 
Savings Blank in the hands of every 
man and woman in your company 
there is still time to join the thousands 
of companies which have added | 
a million employees to the Payroll 


The U. S. Government does not pay for this advertising. The Treasury De- 
partment thanks, for their patriotic donation, the Advertising Council and 


GEOPHYSICS 


Please mention GEopHysiIcs when answering advertisers 


Savings Plan through person-to-person 
canvasses. 

Phone, wire or write to Savings Bond 
Division, U. S. Treasury Department, 
Suite 700, Washington Building, Wash. 
ington, D. C. Your State Director will 
give you all the help you need—appli- 
cation blanks, posters, envelope stuff- 
ers, personal assistance. 

Give your employees an opportunity 
to save for their future and at the same 
time, help to maintain America’s eco- 
nomic security—put an application 
blank in their 4 
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New A.A.P.G. Book! 1951 


POSSIBLE FUTURE 
PETROLEUM PROVINCES 
OF NORTH AMERICA 


UP-TO-DATE PICTURE OF UNDISCOVERED RESOURCES. SYMPOSIUM REPRE- 
SENTING THE WORK OF 17 GEOLOGICAL SOCIETIES, 13 FEDERAL, STATE, AND 
PROVINCIAL GEOLOGICAL SURVEYS, OR ANALOGOUS ORGANIZATIONS, REP- 
RESENTATIVES OF 33 OR MORE OIL COMPANIES AND 8 UNIVERSITIES, AND 
24 CONSULTING GEOLOGISTS. 


e First published in A.A.P.G. Bulletin, February, 1951 
® Complete symposium, 6 x 9 inches, bound in cloth 
© 360 pages, 153 line drawings 
PRICE: $4.00 POSTPAID 
(Member price, $2.50) 


THE AMERICAN ASSOCIATION OF PETROLEUM GEOLOGISTS 
BOX 979, TULSA 1, OKLAHOMA, U.S.A. 














STRUCTURE OF TYPICAL 
AMERICAN OIL FIELDS 


A SYMPOSIUM ON THE RELATION OF 
OIL ACCUMULATION TO STRUCTURE 





VOLUME III BY MANY AUTHORS EDITED BY J. V. HOWELL 


This is the newest (1948) book in the Structure Series, It contains descriptions 
of oil and gas fields characterized by distinctive and peculiar features of 
stratigraphy, structure, and accumulation. 


516 pages 6 x 9 inches 219 illus. Clothbound 


PRICE: $4.50 POSTPAID 
(Member price, $3.50) 


THE AMERICAN ASSOCIATION OF PETROLEUM GEOLOGISTS 
BOX 979, TULSA 1, OKLAHOMA, U.S.A. 
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Braithwaite & Caldwell 


Enid, Oklahoma,U. S. A. 


Al Braithwaite Walter T. Caldwell 
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CHAPTER EIGHT in the Fascinating Story of the Search for Oil 


In 1938 The electrical well logging method, 


pioneered and developed by Conrad and Marcel Schlum- 
berger as early as 1928, was thoroughly proved and ac- 
cepted during its first decade of service, and by 1938 “a 
Schlumberger” had become synonymous with the search 
for new oil reserves. 

“This method, which consists essentially in measuring 
the electrical resistivity and potential along the walls of a 
drill hole by introduction of a series of electrodes, not only 
gives exceedingly important information as to oil, gas and 
water content of formations penetrated, but also yields a 
log which, for correlation purposes, is often superior to 
any other type of log and which in specific cases may be 
of even greater utility than a complete core section,” re- 
ported E. De Golyer in a lecture delivered at Princeton 
University in 1939. 

The early 
electrical 
were plotted by 
hand, 


logs 


as pic- 





tured in this illustration of recording in California in 1932. 
Present day logs are recorded photographically on film 
and four or more curves can be recorded simultaneously 


on the same log. 


195 1 @ @ @ Since new and improved instruments are 


still the key to future progress in the endless search for new 
oil reserves, General maintains extensive instrument labo- 
ratories of its own to develop and improve dependable 
equipment to meet specific needs for various areas. In the 
hands of capable crews who are qualified by training and 
experience to accurately compile and interpret seismograph 
data, General instruments have helped many operators to 
determine and locate conditions favorable to finding 

new oil _ re- 
Let 
General’s crews 


serves. 


help you. 


GEOPHYSICAL COMPANY 


GULF BLDG. 


A PRINT OF THIS 


HOUSTON 


ILLUSTRATION SUITABLE FOR FRAMING WILL BE SENT ON REQUEST 
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Transformers of 
Uniform Performance 
Developed through 
Years of 
Precision Manufacture 








Thermador Geophysical 






Transformers feature: 


(1) hermetic sealing, (2) 






hum-bucking construction, 






(3) close tolerances. 










WRITE FOR CATALOG 
Address: Geophysical Department 


THERMADOR 


ELECTRICAL MANUFACTURING CO. 
5119 District Boulevard 
Los Angeles 22, California 





Distributed by 


A. R. BEYER & CO., 1318 | RADIO, INC., 1000 S. Main | WILKINSON BROS., 2408 | MONTAGUE RADIO. 220 
Polk Ave., Houston 2, Tex. | St., Tulsa 3, Okla. Ross Ave., Dallas 1, Tex. | Willow Street, Beaumont, Texae¢ 
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W BIG CASH PRIZES IN 
the NEW KOENIG 


TRUE EXPERIENCE CONTEST! 
$100 - $75 + $50 - Lots of $10 and $5 Prizes 


You've saved many a day with your handy, helpful KING Winch up front on your 
Jeep or truck. You were mighty glad your faithful KING Front-Mount was there, 
and so are we. . 








We want to give you some hard cash money just for telling us about a real whing- 
ding of an experience with your KING Winch . . . how it got you out of the soup. 
You'll cash in on the gravy. 


It’s as easy as shooting holes in the sand. No formal rules to follow, entry blanks to 
fill out or carton tops or winch wrappers to send in. All we want is your story and 
a photograph of your KING Winch in action. 


$100 First Prize for best story and action shot of incident. 


$75 Second Prize © $50 Third Prize © $10 for all Honorable Mention stories 
and negatives. © $5 for all Honorable Mention stories and prints. 





The number of Honorable Mention awards will be 
determined by the judges, dependent upon the amount 
of acceptable material received. 


SIMPLE RULES 


Everyone is eligible to enter this cash money contest, 
from truck driver to party chief, exploration company a ee ge er err 








or oil company president. 


Entrants will be accepted from anywhere in the world, 
wherever KING Winches are used. You can bet your 
walking shoes too, that wherever there are geophysical 
crews, there are also KING Winches. Stories must be 
submitted in the English language. 


All entries will be judged on the interest of the inci- 
dent, and not on the literary merit of the entry. 


Representatives of our advertising agency and major 
publications will judge the contest, and decisions of 
the judges are final. 


This big cash money contest closes at midnight Decem- 
ber 31, 1951. 


ENTER NOW! 
® 





over 500 words on how your KING 
Front-Mount Winch saved the day 
for you. You don’t need to be 
literary . . . just give the facts. 

2. Get a glossy print or negative 
of your KING Winch installation 
(preferably acticn shot), to illus- 
trate your story. If you do not 
have a photograph of your winch 
installation at the time of the in- 
cident, you can stage a shot under 
similar conditions. 

3, Send your story and photograph 
or negative, w:th your*name and 
address, and your own profile, if 
possible, to address below. 

All entries become our property. 


None will be returned. 





IRON WORKS 





2214 Washington Avenue 


Houston 10, Texas 
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%& PHOTO GRAVITY 


REPUBLIC EXPLORATI 3 


TULSA, OKLAHOMA Mit 
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I. you’RE looking for minerals or oil you 
can go about it in two ways: send out ground parties to 
make surface studies . . . or better still, direct your ground 
exploration with an air survey*. 


Surveys over land or water areas by AERO with the air- 
borne magnetometer are fast. They cost a fraction of 
ground surveys. They spotlight areas of greatest interest. 


Te iia 


For example, take the recent iron ore discovery in Penn- 
sylvania—one of the most important finds in recent years. 
Ground parties had explored the area continuously since 
1917; there were no surface indications. But AERO’s survey 
pointed to a massive ore body 1500 feet below the surface! 
“The discovery was made possible by a new aerial survey 
method using the magnetometer,’ says the mining company. 


Sapticeanptee een 


What was the cost of this survey? 7¢ an acre! 


Yes, 7¢ an acre! And the 900,000 acre survey was flown in 
only 30 days. It’s typical of how AERO cuts costs and speeds 
exploration programs. Moreover, AERO’S magnetic records 
are accurate to better than one part in 10,000—a depend- 
ability achieved through AERO’s more than a quarter 
million miles of aerial exploration. 


With the airborne magnetometer AERO has helped uncover 
important new finds of chrome, asbestos, ilmenite and 
other minerals . . . and delivered important geologic infor- 
mation for oil companies. 


BLNS 
rtereratate®, 
SSSR 


4 


xD 
o 





So 
S555 


Call on AERO’S worldwide mapping experience 
when you plan your next exploration project. 


AERO 


SERVICE CORPORATION 


PHILADELPHIA 20, PENNSYLVANIA 
*There’s a third way. You could Oldest Flying Corporation in the World 


> 


KKK 


© 





KK 
ORR 
KKK 


SRK OL 








TOPOGRAPHIC MAPS ® PLANIMETRIC MAPS 
PRECISE AERIAL MOSAICS 

AIRBORNE MAGNETOMETER SURVEYS 
RELIEF MODELS 

COLOR PHOTOGRAPHY 











hire a mystic and his divining rod! 
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HERE’S WHY DRILLERS CALL 
“BLUE DEMONS” ROCK CUTTERS 


This wide variety of rock cuttings shows why Hawthorne “Blue 
Demon” Bits efficiently drill 90% of the medium and hard forma- 
tions encountered by geophysical crews. Drillers save on roller bits, 
particularly in standstone and cemented formations. 


“Blue Demons” surpass all similar bits in drilling exceptional 
footage and holding gauge size in harder formations because of 
exclusive welding process and cutting tips. Each bit finger cuts only 
a small part of the hole and takes only a small part of the load. 


These same “Blue Demons” drill soft formations far faster than 
conventional bits. The result is faster drilling, fewer round trips 
and drastically reduced bit costs. 


Drillers on 900 drills the 
world over are drilling faster, HAWTH ORNE 
with fewer bit changes, - “BLUE DEMON” 
reduced drilling costs, wit 
long lasting Hawthorne “Blue ROCK CUTTER BITS 
Deane. te ' REPLACEABLE BLADES 





Use “Blue Demons” for ” 497 
rock bit efficiency at “fish tail” 1 $ -10 Sl 
cost. WRITE FOR ILLUSTRATED CATALOG 





HERB SP SMUAUS ALA 


P.@. BOX 7366, HOUSTON 8, TEXAS 





INC. 
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Post-war model for 
use with either 
vertical or horizontal 
systems. 


The HEILAND RESEARCH CORPORATION has been 
appointed exclusive American sales and service rep- 
resentative for the internationally known SCHMIDT- 
ASKANIA line of magnetic prospecting equipment. 


*« Field magnetometers * Calibration coils 


* Photoelectric remote variation recording equipment 


Write for dependable instruments 


complete details 


HEILAND 
RESEARCH 
CORPORATION 


132 East Fifth Ave. 
Denver, Colorado 
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ilures with 


ECTRIC PLANTS 


When power fails, business comes to a stand- 
still. Materials in process are wasted, men 
and machines stand idle, lives may he en- 
dangered. Prevent the disastrous effects of 
power failure with a low-cost,easy-to-install 
Onan Standby System. The generating plant 
starts automatically when commercial pow- 
er fails, takes over the power load wit “a 
interruption, stops automatically when 

er is restored. Many models: 400 to 3 900 
watts, A.C, 


arrison 


Equipment Ca, Ine. 


1422 San Jacinto e« HOUSTON e Phone ATwood 6315 
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Te, = 
OSE AIR B0ATS- kine, - 


WINCHES ON ALLTRUCKS 
PORTABLE -. INSTRUMENTS EA 


300 GALLON WATER TANKS ON DRILLS 
PONTOONED EQUIPMENT lig x7 
ELECTRIC te. SHOT HOLE LOGGERS | 


PROVEN + ECONOMICAL ”~ DEPENDABLE + FAST 


EXPERIENCE 
SUPERVISORS, 16 YEARS + PARTY CHIEFS, 10 YEARS + FIELD CHIEFS, 12 YEARS 


(BASED ON 6 CREW AVERAGE, 1960) 


MARINE 


EXPLORATION CO. 


3732 WESTHEIMER ROAD LLOUSTON PHONE JUSTIN 1521 
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CERTAIN VARIATIONS OF 


POROSITY 


Produce inviting and 





profitable targets for 
EXPLORATION by 


GRAVIMETRIC 
METHODS 


FOR AN EXPLANATION 
WRITE OR CALL 


KLAUS 
EXPLORATION CO. 


TELEPHONE 2-155! BOX 1617 
LUBBOCK, TEXAS 
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Pa ea a ce his motto was, 
“Git thar fustest. with - cote In the field of oil ex- 






k led sae dive. you reliable information as 
sc Ds aorae apaiate eal use of more 
oe » saving you time and money! 






Next ste; & as , many y successful oil men do, let a Mayes- 


Bevan Gravity Meter Survey take the vital first step in your 
oil exploration. 
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INCE 1932 Independ- 
ent has had from two 
to eighteen crews in the field at all 
‘times... operating from the swamps of 


Louisiana to the Rocky Mountains, in D. D. Doty has been work- 

coastal tidewaters, tropical jungles, and pore in: through the Inde- 

ou: deater eae pendent organization for 

. J y 18 years. Now a senior 

The down-to-earth experience gained party chief, his service is 
perience gai “ ‘ 

through successfully overcoming special typical of the extensive 


experience of all key per- 


problems and adverse conditions pro- 
sonnel, 


vides the “know how” that gets every 

job done efficiently and makes Inde- 

pendent’s service profitable to our cli- 
ents. You are invited to 
consult with us about your 
exploration program. 








OE 


Geophysical Surueys 1973 WEST GRAY % HOUSTON, TEXAS 


ESTABLISHED 1932 
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For the Right Explosive 
for Seismic Work 





Party Chiefs jy _« q 
Depend on ATLAS |E e @ 














Every seismic shot represents a big investment 
in time and money. Many “doodlebugs” protect 
this investment—make sure of getting explosives 
and blasting supplies that are just right for the 
job—by consulting Atlas when they plan their 


shooting. 





They know that Atlas technical men are thor- 
oughly familiar with seismic work, and have the 
experience to recommend the explosive and det- 
onator that are just right for a particular job. 
What is more, “doodlebugs” know that Atlas 
service is tops. They know they can depend on 
Atlas to supply the right explosive, at the right 
place, at the right time . . . for seismic work on 


dry land, in swamps or marshes, or under water. 








SU Wilmington 99, Delaware 





Makers of dependable explosives for seismic work 


PETROGEL* TWISTITE* MANASITE* 


high-velocity dynamite fast-coupling device extra-safe detonators 











* Reg. U. S. Pat. Off. 
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(la uh id tnbide Mat count! 


Quality in a watch shows up in the accuracy of its readings. And that 
accuracy depends upon the balanced adjustment of precision parts. 





Quality and reliability in a seismic survey show up in the accuracy 
of its findings, too. The best available equipment; trained, experienced 
personnel; and the skilled interpretation of data are important parts 
of a good survey . . . but we believe that in addition to these parts, 
the close, personal attention of competent, supervisory management 
personnel must provide the balance and coordination that insure the 
highest degree of positive results. 


This is the service we offer management interested in integrity of 
operation and scientifically accurate and dependable seismic surveys. 


DALLAS, TEXAS 


M.C. Kelsey E.F.McMullin J. F. Rollins G.W. Fisher J.C. Spiers 
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For Results Use 


RUSKA MAGNETOMETERS 


Ruska Magnetometers are of latest 
design, based on most recent de- 
velopments. They excel in precision, 
workmanship and materials. Easy 
to operate and maintain, Ruska 
Magnetometers have a fine ap- 
pearance and a durable finish. 














@ Vertical 
@ Horizontal 
@ Recording 


Level vials for Ruska Prospecting, Observatory and Expedition instruments are rigidly tested 
before and after mounting. The Ruska Level Tester (Cat. No. 5995) used for this purpose 
is read to one second per scale division; one-tenth second can be estimated. 


ASK FOR ILLUSTRATED CATALOG 


f | S K A INSTRUMENT 
CORPORATION 
4607 MONTROSE BLVD. HOUSTON 6, TEXAS 
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Geophysical Engineering 


Seismic..Gravity.. Magnetic 
Surveys 


ROBERT H. RAY JACK C. POLLARD 


2500 BOLSOVER ROAD e HOUSTON, TEXAS 
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Special Explosives 
— for Geophysical 
Exploration 


For seismic exploration on land or 
under water, Hercules has devel- 
oped a full series of special explo- 
sives to give accurate recordings 
under the varying conditions en- 
countered by field crews. These 
explosives are quickly available 
throughout the oil-producing 
areas. Write for name of nearest 
distributor. 


VIBROGELS® 


... these special gelatins are de- 
signed to give accurate records and 
uniform results under severe condi- 
tions. 


GELAMITE® S$ 


... economically replaces ‘‘Vibro- 
gels’”’ for some seismic work. 


SPIRALOK® 


... time-saving rigid cartridge as- 
sembly. 


VIBROCAPS® 


. . no-lag seismic detonators. 


TITAN® CAPS 10V and 20V 


... for use where extra-strong det- 
onators are needed. 


HERCULES POWDER COMPANY & 


INCORPORATED 


917 King Street, Wilmington, Delaware 
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4 fact-packed books for the Geophysicist ... 


New facts about 
three - quarters 
of the earth 


How to apply 
our new un- 
derstanding of 
earth waves 


* “the finest sin- 
gle volume 
available” 


Written from the 
modern _physi- 
co - chemical 
viewpoint 


MARINE GEOLOGY 
’ By PH. H. KUENEN, University of Groningen 


With the rapid increase in our understanding of marine geology, facts 
about the ocean floor are becoming steadily more important to geophysi- 
cists. The latest geological information about submarine landforms and 
their interpretation, sedimentary cover of the sea floor, the structure of 
atolls, and other topics of special interest to geophysicists are among 
the subjects covered. 


In addition to reviewing facts gathered in recent studies by oil com- 
panies, oceanographical institutions, and a few individuals, Dr. Kuenen 
gives a clear picture of the far greater number of problems still to be 
solved. 1950. 568 pages. $7.50. 


EARTH WAVES 
By L. DON LEET, Harvard University 


Dr. Leet reviews our present knowledge of elastic waves propagated 
through the earth, and shows how these facts can be used in prospecting, 
mapping earth structures, and other fields. The data he includes on two 
new wave types, on the interpretation of refraction profiles, and on 
microseisms have never been published before, or have been published only 
in periodicals of limited circulation. The four chapters cover: the measure- 
ment of earth waves; observed types of earth waves, transmission of earth 
waves, and microseisms. A Harvard Monograph in Applied Science... . 
1950. 122 pages. $3.00. 


ECONOMIC MINERAL DEPOSITS, Second Edition 
By ALAN M. BATEMAN, Yale University 


Dr. Bateman gives a comprehensive account of the occurrence and form 
of mineral deposits, based on his first-hand knowledge of these deposits in 
all parts of the world. He completely revised the second edition to em- 
phasize the interpretation of mineral deposits and give many practical 
conclusions for ore finding. 

* Vernon E. Scheid of the University of Idaho called the second edition: 
. . . the finest single volume on mineral deposits that is available to the 
general student in America.” ... 1950. 916 pages. $7.50. 


Introduction to 
THEORETICAL IGNEOUS PETROLOGY 


By ERNEST E. WAHLSTROM, University of Colorado 


THE MINES MAGAZINE said: 

“The author has stressed the quantitative approach to the subject. He 
presents a survey of the theory and practice of geophysics and of physical 
chemistry—especially the phase rule... . 


“This excellent book goes a long way toward remedying that situation 
(the lack of a reliable reference book) and will fulfill a great need... . 
for those concerned with the practical application of petrological knowl- 
edge.” 1950. 365 pages. $6.00. 
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TRIAD Geoformers Will Meet or 
Exceed the Most Exacting Require- 
ments of Your Geophysical 
Equipment 


TRIAD GEOFORMERS are trans- 
formers expressly designed for the 
exacting requirements of geophysical 
prospecting equipment. They are not 
merely adaptations of standard com- 
mercial type transformers. From 
“Trialloy” shielding to ‘Climatite” 
vacuum impregnation and Triad her- 


metic sealing, they are designed and 
tooled for geophysical work. 

Write for Catalog GP-51. It will help 
you in the design and constructicn 
of your geophysical equipment. 
Sample quantities of Triad Geoform- 
ers are readily available from distrib- 
utors listed below. 





Sold and Stocked by: 


A. R. Beyer Co., 1318 Polk St., Houston, 
Texas—Electronic Supply Co., 2615 . East 
Foothill Blvd., Pasadena, California—Radio 
1 hy Inc., 1000 So. Main St., Tulsa, Oklahoma— 
Wilkinson Bros., 2406 Ross Ave., Dallas, 


= MFG CO Texas. 





2254 Sepulveda Blvd., Los Angeles 64, Calif. 
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Charge Anchors vs. 
eee Shot Hole Casin g 


The Charge Anchor embodies means and methods where- 
by all of the charges may be located at the desired points 
a in an uncased bore and separately anchored therein and 
Soemagoasite Geoteent fired one at a time at desired intervals from below upwardly; 
thus dispensing with casing in the hole and correspond- 
ingly reducing the time and expense necessary to case the 
bore and to thereafter withdraw the casing. 
Better results can be obtained in loose or soft formations; such as sand, silt, gravel, soft 
clay, etc. More profiles at less cost and less work; in some cases Anchors will cut the cost 
as much as 75% per hole. 
Often it is difficult to use casing because sand fills the pipe after the first shot is fired. We 
can defeat this problem. 
If you want two shots or more in soft formations, call us for a demonstration in the field. 
There will be no charge except regular sales price of our Anchors; and if the method does 
not work, there will be no charge. 
This method is past the experimental stage. It has been in use by numerous crews for a 
period well over five years. Names of satisfied customers can be supplied upon request. 
This method is ideal for formations found in South Texas, Louisiana, & Mississippi with 
charges up to 25 pounds per shot. Anchors will not work in rock, hard shale, and sand 
stone. Call us collect. Tell us your problems. We can soon tell you IF we can be of service 
to you. 
Take advantage of this valuable method. 


Mr. O. V. GEE, Representative JERRY D. WILLIAMS, Owner 
in Lake Charles, La. Box 707, Chickasha, Okla. 
Phone 66233 Phone 3617 
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Eliminate Electric Power 
Line Hazards In 
Exploration Shooting 


? 











With this Firing Line 
Catcher shots can be fired 


caught, and the bar is con- 





trolled in any desired di- 






near electric power lines rection to a limited height 






safely. and radius from the bore. 






Construction of Catcher: We have forty-seven 





Catchers in use at this time 





A bar approximately six 
| feet in length with hooks 


arranged alternately 





—several have rented for 






more than four years of 






continuous service with 





around the bar. It has a 





outstanding results. 





cross arm and _ control 


weight attached to the top, 








, : ’ Patent Number 2,472,997 
an Soe ee neveling weight. Firing Line Catchers are not for sale. 
’ The catcher is designed to be expelled They are placed with crews on a rental 
s from the hole under hard shots and basis of $25.00 per month and a mini- 
hole blow. At upward movement the mum of four months rental. 


circuit is positively broken, wires are Call on us for a free demonstration. 


EXPLORATION ANCHOR COMPANY 
BOX 707, CHICKASHA, OKLAHOMA 
Phone 643 or 3617 


Please mention GropHysics when answering advertisers 


Ye — 
o 













GEOPHYSICS the Journal of the Society of Exploration Geophysicists 


 Pelly 





GROPHYSICAL ENGINEERING COMPANY 





ee INTONIO TEXAS 


SEISMIC ° GRibIEY ° MAGNETIC SURF EXS 


. of world-wide geophysical 


engineering experience go into 


every job we undertake. 
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geophysical 


our proven marine and land cables 


adapted to your requirements 7 


WRITE FOR CATALOGUE 


md . aAUuUU ac lur 4 C c i} 4 ty 


5616 LAWNDALE * HOUSTON 3, TEXAS 
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Dr. William H. White, Associate Professor, University of 
British Columbia, Canada, reproduced by permission. 


WATTS Microptic Alidade 


Glass circle can be read to one minute and includes Beaman Stadia reduction 
factors. Telescope, scale and bubble readings all taken from same head position. 
Construction robust yet exceptionally light. 

Write for List GS/66 to: 


JARRELL-ASH CO. 
165 Newbury Street, Boston, Mass. 


WAT ITSA. tte. e ENGLAND 
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Computing procedures are recommended 
to the client to be followed in areas of 


Vz = Vo + KZ steep dips. The scheme, roughly indicated 
»_ Kz in outline above, may be used in steep dip 
Z- 2Vz areas where a linear increase in instan- 

Vv. Kt taneous velocity with depth holds to an 
Re K SINH 3 acceptable approximation. Appropriate 


schemes, with necessary charts, are devised 
for any known velocity function pertaining 
to any particular area. 


The latest proved developments in instrumentation and seismic 
theory are available to you when you employ Midwestern. 
We welcome your inquiry regarding services of our contract 


seismograph field crews. 





MIDWESTERN 


GEOPHYSICAL LABORATORY 


"TUN UY 
3400 S. Harvard a Tulsa, Oklahoma 


Manufacturers of Precision Seismograph Instruments 
Operators of Contract Seismograph Field Crews 
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é An easily read, clearly illustrated text on 


Modern Geophysical Techniques 








EXPLORATION GEOPHYSICS 


by J. J. Jakosky, Sc. D. 


Thirteen fact-packed chapters 
fully cover all contemporary 
methods; plus permit, trespass 
and insurance problems. A 
basic textbook for every 
geologist, geophysicist, engi- 
neer and physicist concerned 
with exploration, well logging 
and production. Adopted by 
many leading universities. 


In 1200 pages and with 707 
illustrations, the 1950 revised 
Expleration Geophysics 
covers the entire field of 
exploration by modern geo- 
physical methods. It is con- 
cisely and clearly written by 
an internationally known geo- 
physicist, in close collabora- 
tion with 39 other leading 
authorities. 





Send your money order or check for $12.50 for a copy of Exploration Geophysics on 


TRIJA PUBLISHING COMPANY, TRIJA BUILDING, LOS ANGELES 24, CALIFORNIA 


3 H 
v 1 
H 5-day approval. If you are not fully satisfied, merely return the book in its original 
| condition and your money will be promptly refunded. ' 
t j 








CONTENTS 


CONTRIBUTORS AND COLLABORATORS, Preface 
6 pages 
CHAPTER I 18 pages 


INTRODUCTION—History of. geophysical meth- 
ods. Contemporary workers and developments. 
Trends. Development of future methods. 

CHAPTER II 52 pages 
GEOLOGIC AND ECONOMIC BACKGROUND 
OF EXPLORATION GEOPHYSICS—Factors 
governing application and choice of methods. Field 
technique. Use of Geophysics in prospecting for 
petroleum and metallic minerals, water supply, and 
engineering applications. 

CHAPTER III 186 pages 
MAGNETIC METHODS—Theory, instrumenta- 
tion, field operations, calculation and interpretation. 
Airplane, helicopter and ship-borne operations. 

CHAPTER IV 190 pages 
GRAVITATIONAL METHODS—Theory, instru- 
mentation, field operations, calculation and inter- 
pretation. "Pendulum, torsion balance, gravity meter 
measurements. Land and under-water meters. Div- 
ing bell. Leveling and photogrammetric mapping. 

CHAPTERS V and VI 202 pages 
ELECTRICAL METHODS—Theory, instrumenta- 
tion, field operations, calculation and interpretation. 


——— equipotential resistivity and inductive 
met 
CHAPTER VII 300 pages 


SEISMIC METHODS—Refraction and reflection 
techniques. Theory, instrumentation, field opera- 
tions, calculation and interpretation. Dip, correla- 
tion, pulse, spot and geological correlation, con- 
tinuous profiling. Fault mapping. Low velocity layer 


corrections, Well velocity shooting. Graphical com- 
puters, Instrumental analysis, Drills, Explosives. 
Air shooting. Off-shore seismic operations. Radio 
surveying. 

CHAPTER VIII 28 pages 
CHEMICAL METHODS—Physical _ principles. 
Sampling, analysis, field operations and interpreta- 
tion. 

CHAPTER IX 21 pages 


THERMAL METHODS—Theory of heat flow. 
Vertical gradients and lateral variations. Interpreta- 


tion. 
CHAPTER X 29 pages 
RADIOACTIVITY METHODS—Particles and 
uanta. Radioactive elements. Instrumentation. 
ield techniques. Interpretation. 
CHAPTER XI 106 pages 


BORE HOLE INVESTIGATIONS-—-Resistivity, 
self potential, temperature, photoelectric, loggin 
techniques. Dipmeter, photoclinometer, sidewal 
sampling, section gauges, radioactive markers. 
Radioactivity well logging. Examples and interpreta- 
tion. 

CHAPTER XII 21 pages 
PHYSICAL PRINCIPLES APPLIED TO PRO- 
DUCTION PROBLEMS—Theory of oil produc- 
tion. Reservoir pressure. Methods for determining 
fluid level. Solution of pumping problems. 

CHAPTER XIII 20 pages 

LAND TENURE, PERMIT AND TRESPASS 
PRACTICES: INSURANCE for Geophysical op- 
erations, PATENTS governing Geophysical Instru- 
mentation and Exploration. 


INDEX—Name and Place. Subject. 17 pages 
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YOUR FIRST STEP TOWARD 
SUCESSFUL COMPLETIONS 


The clear, definite character data 
obtained with Reliable seismograms are not the result of chance. Reliable’s 
precise instruments, based on years of laboratory and field testing, in 
the hands of our experienced crews, working under the direct supervision 
of one of the partners, mean the best in high quality work for every 


exploration dollar spent. 


Reliable’s 32-trace seismograms, 16-traces simple and 16- traces 
mixed, are now the clearest obtainable. But even the best instruments 
must be used properly, and that knowledge and skill are part of our 
service to you. Your assurance that you can depend on Reliable is the 
high interest and the skill which Reliable’s crews bring to your job, and 


the unusually close supervision. 


Write for the availability of crews that will assure you of a job 


well done. 








Glenn M. McGuckin Perry R. Love 
Phone 108 Yoakum, Texas P. O. Box 450 
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PORTABLE 
SEISMOGRAPH 
INSTRUMENTS 














IN CANADIAN MUSKEG AND BRUSH COUNTRY, geophysical 
prospectors have proved both the performance and portability of 
SIE seismographs. Yes, these SIE seismograph instruments can be 
carried easily anywhere a man can go... and performance features 
and characteristics are equal to full-sized truck-mounted instruments. 
Before you undertake difficult areas . . . jungle, swamp, mountains or 
snow... get all the facts about the SIE P-11 Portable Seismograph. 
It gets the records! Write for complete details. 


SOUTHWESTERN INDUSTRIAL ELECTRONICS 


POST OAK ROAD HOUSTON 
434 SEVENTH AVE. EAST — CALGARY, ALBERTA, CANADA 


CO. 


19, TEXAS 


2831 
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For Safe, Efficient, 


Economical Drilling | A new volume which will be indis- 
oe specify pensable to field geologists 
SPANG STRUCTURAL GEOLOGY 
CABLE TOOLS OF NORTH AMERICA 
anatnecanehens nel By A. J. Eardley, University of Utah 


Spang Cable Tools 
stand out as the top 
performers in their 
field for: 


@ Shot Hole Drilling 


© Geological 
Exploration 


This is the first volume to describe in 
detail the structural evolution of an en- 
tire continent: formation and constitu- 
tion of mountain systems, basins, arches, 
and volcanic archipelagos, beveling of 
highlands, filling of basins. Double- 
column pages, 11144” x 814”, provide 
the equivalent of two ordinary volumes, 
and permit large maps and other draw- 
ings. Some 750 illustrations and 16 
paleotectonic maps in full color, 


@ Petroleum 
Production 


@ All types of Cable 
and Churn Tool 
Drilling 





Try Spang today! 620 double-column pages $12.50 
S Fe A ii i CABLE For sale by 
TOOLS SOCIETY OF EXPLORATION 
N NY. BUTLER, PA. GEOPHYSICISTS 
SPANG & COMPANY, LER, PA Box 7248 Tulsa 18, Oklahoma 


For Sale by Dealers Everywhere 




















A practical guide to... 


the principles and techniques used in oil and 
gas exploration 














pt for practicing petroleum geologists, eng:neers, and 
others concerned with oil exploration and oil production, 
this comprehensive book gives a clear, thorough coverage of 
the principles, processes, and methods of practical im- 
portance and use in petroleum geology. Avoiding long, 
detailed descriptions of oil fields and their stratigraphy, 
the book emphasizes those characteristics of oil and gas 


Just Published! 


fields which make locating them easier and quicker. Principles of 
Packed into this 509-page reference book is a wealth of leg orm 
much needed information on a wide range of pertinent 

subjects—from the origin and accumulation of oil, im- GEOLOGY 


portant geologic structures, through classification of 
fields, stratigraphic features such as reefs—to geophysics, 
1 i tr] y i Vi ; i ae Associate Professor of Geology 
lithologic, electrical and radioactivity well logs—and Meuiediliniady Gonds Medan 


methods of locating wells. Refer to this book for accurate, College of Texas 
dependable facts, methods, and al a RO 


By WILLIAM L. RUSSELL 











data on petroleum geology—the 

kind of information that will aid For sale by 

in more successful, economical SOCIETY OF EXPLORATION GEOPHYSICISTS 
discovery and production of oil Box 7248 Tulsa 18, Oklahoma 
and gas. 
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5300 BROWNWAY ROAD 
HOUSTON, TEXAS 


MOhawk 4531 
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217 ALBERTA CORNER BLDG. 


CALGARY, ALBERTA 


2-5038 
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AMERICAN PAULIN SYSTEM 








MODEL M-1 . . . Range 
6,000 feet (—1000’ to-+- 
5000’) in intervals of 1- 


MODEL M-2 . . . Range 
10,000 feet (0’ to 10,- 
000’) in intervals of 2’ 


MODEL M-5 . . . Range 
15,000 feet (0’ to 15,- 
000’) in intervals of 5’ 


MODEL MM-1. . . Range 
5,000 meters (0 to 5,000 
meters) in intervals of 1 
fe meter 


47 $300 EACH with leather 
case, Thermometer, Mag- 
nifier, and Operational 
Procedures. 








Also 
available 
in Metric 
System 












a Don’t guess at a 
altitude readings — 
Only American Paulin 
System Altimeters are 
graduated in easily 
read 1 foot divisions. 







ACCURATE 
DEPENDABLE 
LIGHTWEIGHT 













The World’s Standard... 


TERRA SURVEYING ALTIMETER 


$200 EACH with leather 
case, Thermometer, Mag- 
nifier, and Operational 
Procedures, 


Literature and Technical Publications Available on Request 


nic”o AMERICAN PAULIN SYSTEM TERR 


1847 S. FLOWER © LOS ANGELES 15 





MODEL SA-2 . . . Range | MODEL SA-5 .. . Range 
10,600 feet (—900’ to-++ | 15,000 feet (—500’ to-+ 
9,700’) in intervals of 5’ | 14,500’)in intervals of 10’ 


MODEL SA-1 . . . Range 
4,360 feet (—760’ to-+ 
3,600’) in intervals of 2’ 
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DIFFICULT 
AREA 
PROBLEMS? 


LET US HELP YOU 














Multiple plant arrangements are 
made practical by Electro Tech’s 
light but highly efficient and rugged 
Type EVS seismic detector. 












ACTUAL SIZE 


Present day competition demands the best in seismic 
exploration tools and teehniques. The lower operating cost 
and improved performance of the Type EVS seismic 
detector over older, eum bersome equipment soon pays for 
the relatively low Jou are invited to contact 
us for further infow s and other equipment 
manufactured by .E ing complete seis- 
mograph systems. 
We will gladly arrange a 


stration on your prospect. — 


































vee plant demon- 









ELECTRO TECHNI 


504 WAUGH DRIVE 


HOUSTON, TEXAS 
JU-1597 
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GEOGRAPHY IN THE 20TH CENTURY 
Edited by Griffith Taylor, University of Toronto 


Twenty authors, each a specialist in his field, treat the growth, fields, techniques, 


aims and trends of geography. 


From The Contents 


Part I-EVOLUTION OF GEOGRAPHY 
AND ITS PHILOSOPHICAL BASIS 


The French and German School of Geog- 
raphy 

The West Slav Geographies 

Environmentalism and Possibilism 


Part II—THE ENVIRONMENT AS A 
FACTOR 


The Progress of Geomorphology 
Geographical Aspects of Meteorology 
Climatic Influences 

Soils and Their Geographical Significance 
Settlement of the Modern Pioneer 
Geography and Arctic Lands 


Exploration of Antarctica 
Geography and The Tropics 
Geography and Regionalism 


Part III—SPECIAL FIELDS OF 
GEOGRAPHY 


Geography is a Practical Subject 

Geography and Empire 

Racial Geography 

The Sociological Aspects of Geography 

Urban Geography 

Geography and Aviation 

The Field of the Geographical Society 

Geography in Practice in The Federal 
Government, Washington 

Geopolitics and Geopacifics 


56 charts, maps and illustrations with a glossary of geographical terms. Over 600 pages. 


PHILOSOPHICAL LIBRARY 


$8.75 


Publishers 
15 East 40th Street, Desk 595, New York 16, N.Y. 
Expedite Shipment by Prepayment 
Special student bulk rate on 10 or more 


























GE.0O PHYS ECs 
Now Available On 
MICROCARDS 


Volume 3, No. 1 2 cards .50 
Volume 3, No. 3 4 cards 1.00 
Volume 3, No. 4 2 cards -50 
Volume 6, No. 1 3 cards oy 
Volume 7, No. 1 4 cards 1.00 
Volume 8, No. 1 3 cards ey, 
Volume 8, No. 2 4 cards 1.00 
Volume 8, No. 3 3 cards YP. 
Volume 8, No. 4 4 cards 1.00 
Volume 9, No. 1 4 cards 1.00 
Volume 9, No. 2 5 cards 1,25 
Volume 9, No. 3 4 cards 1.00 
Volume 9, No. 4 6 cards 1.50 
Volume 10, No. 1 5 cards 1.25 
Volume 10, No. 2 6 cards 1.50 
Volume 12, No. 1 5 cards 1.25 





SOCIETY OF EXPLORATION 
GEOPHYSICISTS 
Box 7248 Tulsa 18, Oklahoma 
Back Issues of Periodicals 

Journal of the Society of Petroleum Geophysi- 
cists 

1935 Volume VI Number | 

- Geophysics 

1936 Volume | Number |, 2, 3, (complete) 

1937 i I, 2, 3, 4 

1938 W 2 

1939 IV i, 2304 

1940 v 2, 3, 4 

1941 vi a3 

1942 Vil 2, 3, 4 

1943 Vill (out of print) 

1944 Ix (out of print) 

1945 x 3, 4 

1946 Xl 1; 2:3, 4 

1947 Xi 2, 3, 4 
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The Client’s subsurface expectations were “amd 
confirmed by SSC’s reflection survey. run surveys: poperato 
THE RESULT: A new Pennsylvanian reef Ol and Gas 
igi in southeastern Borden County, July 12, en 
exas. 











Summarizing SSC’s Approach to SUCCESSFUL Reef Exploration, the 
following fundamental procedures are recommended: 


1. Complete Geological and Geophysical Coordination. 
2. Detailed Surface and Sub-surface Coverage. 

3. Highest Record Quolity. 

4. Detailed and Thorough Interpretive Technique. 

5. Isopach Analyses. 

6. Velocity Analyses. 
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G.S.I. 
SALUTES 


Field Party 


312 


for outstanding 
Safety Record 
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n recognition of five continuous years of seismic operation with- 










out a single lost-time accident... G.S.I. congratulates Party 312! 
G.S.I.’s seismic parties contributed to a distinguished safety record 
in 1950, with a total of 359 crew-months of accident-free opera- 





tion—30 crew-years. The members of the field crews made this 
outstanding record possible. 


A safety program, of which the above record is an example, is 
a portion of G.S.I.’s long-range planning for efficient operations. 








G.S.1."s Nationol Sofety Film Award Winner—'‘'Moke 
No Mistoke’’—Iis available for free showing. Write for 
porticulors on this outstanding color motion picture. 





P7 Geopnysicat Service Inc. 
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